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We review the basic aspects of electrons in graphene (two-dimensional graphite) exposed to a 
strong perpendicular magnetic field. One of its most salient features is the relativistic quan- 
tum Hall effect the observation of which has been the experimental breakthrough in identifying 
pseudo-relativistic massless charge carriers as the low-energy excitations in graphene. The effect 
may be understood in terms of Landau quantisation for massless Dirac fermions, which is also the 
theoretical basis for the understanding of more involved phenomena due to electronic interactions. 
We present the role of electron-electron interactions both in the weak-coupling limit, where the 
electron-hole excitations are determined by collective modes, and in the strong-coupling regime 
of partially filled relativistic Landau levels. In the latter limit, exotic ferromagnetic phases and 
incompressible quantum liquids are expected to be at the origin of recently observed (fractional) 
quantum Hall states. Furthermore, we discuss briefiy the electron-phonon coupling in a strong 
magnetic field. Although the present review has a dominating theoretical character, a close con- 
nection with available experimental observation is intended. 
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I. INTRODUCTION TO GRAPHENE 

The experimental and theoretical study of graphene, 
two-dimensional (2D) graphite, has become a major is- 
sue of modern condensed matter research. A milestone 
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was the experimental evidence of an unusual quantum 
Hall effect reported in September 2005 by two different 
groups, the Manchester group led by Andre Geim and 
a Columbia-Prince ton collaboration led by Philip Kim 
and H orst Stormer ( Novoselov et "aZl. l2005aHzhang et all 
120051 ). 

The reasons for this enormous scientific interest are 
manyfold, but one may highlight some major motiva- 
tions. First, one may underline its possible technological 
potential. One of the first publications on graphene in 
2004 by the Geim group reported indeed an electric field 
effect in graphene, i.e. the possibility to control the car- 
rier density in the graphene sheet by sim ple application 
of a gate voltage I Novoselov et a/.l . |2004[ ). This effect is 



ground state 



excited state 
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a fundamental element for the design of electronic de- 
vices. In a contemporary publication Berger et al. re- 
ported on the fabrication and the electrical contacting 
of monol ayer graphene samp les on epitaxially grown SiC 
crystals ( Berger et al\ . \206^ . Today's silicon-based elec- 
tronics reaches its limits in miniaturization, which is on 
the order of 50 nm for an electric channel, whereas it has 
been shown that a narrow graphene strip with a width 
of only a few nanometers may be used as a transistor 
( Ponomarenko et all . 120081 ) . i.e. as the basic electronics 
component. 

Apart from these promising technological applications, 
two major motivations for fundamental research may be 
emphasized. Graphene is the first truely 2D crystal ever 
observed in nature and possess remarkable mechanical 
properties. Furthermore, electrons in graphene show rel- 
ativistic behavior, and the system is therefore an ideal 
candidate for the test of quantum-field theoretical models 
that have been developed in high-energy physics. Most 
promenently, electrons in graphene may be viewed as 
massless charged fermions living in 2D space, particles 
one usually does not encounter in our three-dimensional 
world. Indeed, all massless elementary particles hap- 
pen to be electrically neutral, such as photons or neutri- 
nos.^ Graphene is therefore an exciting bridge between 
condensed-matter and high-energy physics, and the re- 
search on its electronic properties unites scientists with 
various thematic backgrounds. 

Several excellent reviews witness the enormous re- 
search achievements i n graphene. In a first step t hose by 
Geim an d Novoselov dCeim and Novoselovi l2007f ) and by 
de Heer ( de Heer et all 2007 ) aimed at a rather global 
experimental review of exfoliated and epitaxial graphene, 
respectively. Furthermor e, the review by Castro Neto 
( Castro Neto et al I l2009^ was concerned with general 
theoretical issues of electron s in graphene. Apar t from 
the review by Abergel et al. ( Abergel et a/.l . [2010t ). more 
recent reviews concentrate on the subfields of graphene 



^ The neutrino example is only partially correct. The 
observed oscillation between different neutrino fla- 
vors (ug o 1/t) requires indeed a tiny non-zero m ass 
llFukuda. Y. et al. (Super-Kamiokande Collaboration)! . 1 199 j) . 
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Figure 1 Electronic configurations for carbon in the ground state 
(left) and in the excited state (right). 



research, which have themselves grown to a consider- 
able size and that require reviews on their ow n. As 
an ex ample one may cite the review by Peres ( Peresl 
l20fflf ). which is concerned with transport properties of 
graphene, o r that by Kotov an d co-workers on interac- 
tion effects ( Kotov et all I20f0h . The present theoretical 
review deals with electronic properties of graphene in a 
strong magnetic field, and its scope is delimited to mono- 
layer graphene. The vast amount of knowledge on bilayer 
graphene certainly merits a review on its own. 



A. The Carbon Atom and its Hybridizations 

In order to understand the crystallographic structure 
of graphene and carbon-based materials in general, it is 
useful to review the basic chemical bonding properties of 
carbon atoms. The carbon atom possesses 6 electrons, 
which, in the atomic ground state, are in the configu- 
ration Is^2s^2p^, i.e. 2 electrons fill the inner shell Is, 
which is close to the nucleus and which is irrelevant for 
chemical reactions, whereas 4 electrons occupy the outer 
shell of 2s and 2p orbitals. Because the 2p orbitals (2p^, 
2py, and 2p2) are roughly 4 eV higher in energy than the 
2s orbital, it is energetically favorable to put 2 electrons 
in the 2s orbital and only 2 of them in the 2p orbitals 
(Figll]). It turns out, however, that in the presence of 
other atoms, such as e.g. H, O, or other C atoms, it is 
favorable to excite one electron from the 2s to the third 
2p orbital, in order to form covalent bonds with the other 
atoms. 

In the excited state, we therefore have four equivalent 
quantum-mechanical states, |2s), \2px), \2py), and |2p^). 
A quantum-mechanical superposition of the state |2s) 
with n \2pj) states is called sp" hybridization. The sp^ 
hybridization plays, e.g., an important role in the context 
of organic chemistry (such as the formation of acetylene) 
and the sp'^ hybridization gives rise to the formation of 
diamonds, a particular 3D form of carbon. Here, how- 
ever, we are interested in the planar sp^ hybridization, 
which is the basic ingredient for the graphitic allotropes. 

As shown in Fig. ^ the three sp^-hybridized orbitals 
arc oriented in the xy-plane and have mutual 120° angles. 




Figure 2 ('aj Schematic view of the sp'^ hybridization. The orbitals 
form angles of 120°. (b) Benzene molecule (CgHe). The 6 carbon 
atoms are situated at the corners of a hexagon and form covalent 
bonds with the H atoms, (c) The quantum- mechanical ground 
state of the benzene ring is a superposition of the two configurations 
which differ by the position of the tt bonds, (d) Graphene may be 
viewed as a tiling of benzene hexagons, where the H atoms are 
replaced by C atoms of neighboring hexagons and where the tt 
electrons are delocalized over the whole structure. 
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Figure 3 (a) Honeycomb lattice. The vectors Si, S2, and 53 
connect nn carbon atoms, separated by a distance a = 0.142 nm. 
The vectors ai and 3.2 are basis vectors of the triangular Bravais 
lattice, (b) Reciprocal lattice of the triangular lattice. Its primitive 
lattice vectors are and aj . The shaded region represents the first 
Brillouin zone (BZ), with its center F and the two inequivalent 
corners K (black squares) and K' (white squares). The thick part 
of the border of the first BZ represents those points which are 
counted in its definition such that no points are doubly counted. 
The first BZ, defined in a strict manner, is, thus, the shaded region 
plus the thick part of the border. For completeness, we have also 
shown the three inequivalent cristallographic points M, M' , and 
M" (white triangles). 



The remaining unhybridizcd orbital is perpendicular 
to the plane. 

A prominent chemical example for such hybridization 
is the benzene molecule the chemical structure of which 
has been analyzed by August Kekule in 1865 ( Kekulel 
[1865. 1866[). The molecule consists of a hexagon with 
carbon atoms at the corners linked by a bonds [Fig. [2] 
(b)]. Each carbon atom has, furthermore, a covalent 
bond with one of the hydrogen atoms which stick out 
from the hexagon in a star-like manner. In addition to 
the six a bonds, the remaining 2p^ orbitals form three 
TT bonds, and the resulting double bonds alternate with 
single a bonds around the hexagon. Because a double 
bond is stronger than a single a bond, one may ex- 
pect that the hexagon is not perfect. A double bond 
(C=C) yields indeed a carbon-carbon distance of 0.135 
nm, whereas it is 0.147 nm for a single a bond (C-C). 
However, the measured carbon-carbon distance in ben- 
zene is 0.142 nm for all bonds, which is roughly the av- 
erage length of a single and a double bond. This equiv- 
alence of all bonds in benzene was explained by Linus 
Pauling in 1931 withi n a quantum-m echanical treatment 
of the benzene ring (jPaulinel ll960D . The ground state 
is indeed a quantum-mechanical superposition of the two 
possible configurations for the double bonds, as shown 
schematically in Fig. [5] (c). 

These chemical considerations indicate the way to- 
wards carbon-based condensed matter physics - any 
graphitic compound has indeed a sheet of graphene as its 
basic constituent. Such a graphene sheet may be viewed 
simply as a tiling of benzene hexagons, where the hydro- 
gen are replaced by carbon atoms to form a neighboring 
carbon hexagon [Fig. [2](d)]. However, graphene has re- 
mained the basic constituent of graphitic systems during 
a long time only on the theoretical level. From an exper- 
imental point of view, graphene is the youngest allotrope 



and accessible to electronic-transport measurements only 
since 2004. 

For a detailed discussion of the different fabrication 
techniques, the most popular of wh ich are the exfoliation 
technique ( Novoselov et al\ . l2005bl ) and th ermal graphi- 
tization of epitaxially-grown SiC crystals ( Berger et al\ . 
|2004( ). we refer the reader t o exi s ting experimental re - 
views ( Geim and Novoselovl 120071 : Ide Heer et al. |,[2003). 
Notice that, more recently, large-scale grap hene has been 
fabri cated by chemical vapor deposition ( Reina et all 
I2OO9D that seems a promising technique not only for fun- 
damental research but also for technological applications. 



B. Crystal Structure of Graphene 

As already mentioned in the last section, the carbon 
atoms in graphene condense in a honeycomb lattice due 
to their sp'^ hybridization. The honeycomb lattice is not 
a Bravais lattice because two neighboring sites are in- 
equivalent from a crystallographic point of view.^ Fig. 
[3] (a) illustrates indeed that a site on the A sublattice 
has nearest neighbors (nn) in the directions north-east, 
north-west, and south, whereas a site on the B sublattice 
has nns in the directions north, south-west, and south- 
east. Both A and B sublattices, however, are triangular 
Bravais lattices, and one may view the honeycomb lat- 
tice as a triangular Bravais lattice with a two- atom basis 
(A and B). The distance between nn carbon atoms is 
a = 0.142 nm, which is the average of the single (C-C) 
and double (C=C) covalent a bonds, as in the case of 



^ This needs to be clearly distinguished from a chemical point of 
view according to which they may be equivalent as in the case 
of graphene where both types of sites consist of carbon atoms. 
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benzene. 

The three vectors which connect a site on the A sub- 
lattice with a nn on the B sublattice are given by 

a 

- \^-y/'6e.j, + Gj^j , da = -aty, 

(1) 

and the triangular Bravais lattice is spanned by the basis 
vectors 



ai 



and 



a2 



\/3e 



(2) 



The modulus of the basis vectors yields the lattice spac- 
ing, a — "v/Sa = 0.24 nm, and the area of the unit cell 
is Auc ~ \/ia? /2 = 0.051 nm^. The density of carbon 
atoms is, therefore, nc = 2/A„c = 39 nm~^ = 3.9 x 10^^ 
cm~^. Because there is one tt electron per carbon atom 
that is not involved in a covalent a bond, there are as 
many valence electrons as carbon atoms, and their den- 
sity is, thus, = nc = 3.9x 10^^ cm~^. As discussed in 
detail below, this density is not equal to the carrier den- 
sity in graphene, which one measures in electric transport 
measurements. 

The reciprocal lattice, which is defined with respect to 
the triangular Bravais lattice, is depicted in Fig. [3] (b). 
It spanned by the vectors 



a, = 



2tt 
\/3a 



and 



3a " 



(3) 



Physically, all sites of the reciprocal lattice represent 
equivalent wave vectors. The first Brillouin zone [BZ, 
shaded region and thick part of the border of the hexagon 
in Fig. [3] (b)] is defined as the set of inequivalcnt points 
in reciprocal space, i.e. of points which may not be con- 
nected to one another by a reciprocal lattice vector. The 
long wavelength excitations are situated in the vicinity 
of the F point, in the center of the first BZ. Furthermore, 
one distinguishes the six corners of the first BZ, which 
consist of the inequivalcnt points K and K' represented 
by the vectors 



±K = ± 



(4) 



The four remaining corners [shown in gray in Fig. [3] (b)] 
may indeed be connected to one of these points via a 
translation by a reciprocal lattice vector. These cristal- 
lographic points play an essential role in the electronic 
properties of graphene because their low-energy excita- 
tions are centered around the two points K and K\ as is 
discussed in detail in the following section. We empha- 
sise, because of some confusion in the literature on this 
point, that the inequivalence of the two BZ corners, K 
and K' , has nothing to do with the presence of two sub- 
lattices, A and B, in the honeycomb lattice. The form of 
the BZ is an intrinsic property of the Bravais lattice, in- 
dependent of the possible presence of more than one atom 
in the unit cell. For completeness, we have also shown, 
in Fig. [3](b), the three crystallographically inequivalcnt 
M points in the middle of the BZ edges. 



C. Electronic Band Structure of Graphene 

As we have discussed in the previous section, three 
electrons per carbon atom in graphene are involved in the 
formation of strong covalent a bonds, and one electron 
per atom yields the tt bonds. The tt electrons happen to 
be those responsible for the electronic properties at low 
energies, whereas the a electr ons form energy ba nds far 
away from the Fermi energy ( Saito et alV Il998() . This 
section of the introduction is, thus, devoted to a brief 
discussion of the energy bands of tt electrons within the 
tight-binding approximation, which was originally calcu- 
lated for the ho neycomb lattice by P. R. Wallace in 1947 
(IWallacel . [l947l) . 



1. Tight-binding model for electrons on the honeycomb lattice 

In the case of two atoms per unit cell, we may write 
down a trial wave function 



(5) 



momentum k. Both ijj'^\v) and ip'^' [r) arc Bloch func 



^k(r) =ak^(-^)(r)+M^f)(r), 

where Ck and 6k are complex functions of the quasi- 
momcntur 
tions with 

V-i'^ W = E e'^''^'''0^''(r + (6) 

where j = A/B labels the atoms on the two sublatticcs 
A and B, and 5j is the vector which connects the sites 
of the underlying Bravais lattice with the site of the j 
atom within the unit cell. The 0(^^(r -I- 5j — R;) are 
atomic orbital wave functions for electrons that are in the 
vicinity of the j atom situated at the position R; — 5j 
at the (Bravais) lattice site R;. Typically one chooses 
the sites of one of the sublatticcs, e.g. the A sublattice, 
to coincide with the sites of the Bravais lattice. Notice 
furthermore that there is some arbitrariness in the choice 
of the phase in Eq. ([6]) - instead of choosing exp(ik • R; ) , 
one may also have chosen exp [ik- (R; — Sj)], for the atomic 
wave functions. The choice, however, does not affect the 
physical properties of the system because it simply leads 
to a redefinition of the weights Ck and 6k which aquire a 
different relative phase ( Bena and Montambauxi |2009|) . 

With the help of these wave functions, we may now 
search the solutions of the Schrodinger equation ifV'k = 
Ek^k- Multiplication of the Schrodinger equation by ipl^ 
from the left yields the equation tp^H^k = ek'0k''/'kj 
which may be rewritten in matrix form with the help 

ofEq. dm 



K,&k)^k 



Ok 



ek(ak:^k)'5k 



flk 



6k y ^ \ 6k 

Here, the Hamiltonian matrix is defined as 



(7) 



(8) 
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and the overlap matrix 



/ I (A) 



,{A) 




(9) 



accounts for the non-orthogonahty of the trial wave func- 
tions. The eigenvalues ek of the Schrodinger equation 
yield the electronic bands, and they may be obtained 
from the secular equation 





\a3 32/ 











^ A" 


ai 
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Figure 4 Tight-binding model for the honeycomb lattice. 



det [Hk - e^5k] = 0, (10) 

which needs to be satisfied for a non-zero solution of the 
wave functions, i.e. for ak 7^ and 6k 7^ 0. The label A 
denotes the energy bands, and it is clear that there are as 
many energy bands as solutions of the secular equation 
([TU)) . i.e. two bands for the case of two atoms per unit 
ceU. 



a. Formal solution. Before turning to the specific case of 
graphene and its energy bands, we solve formally the sec- 
ular equation for an arbitrary lattice with several atoms 
per unit cell. The Hamiltonian matrix ([8|) may be writ- 
ten, with the help of Eq. ([6]), as 

n^ = N{e<~^^s^ + q) (11) 

where (5y = 5^ - S,), 

^ ^ e'*--^' f d\ (r)</)(^) (r + % - R/) = ^ 

(12) 

and wc have defined the hopping matrix 
4^ =^e''"^' f d^rcl)^'>{r)AV(P^^\r+6,,-Ri) . (13) 

Here, N is the number of unit cells, and we have sep- 
arated the Hamiltonian H into an atomic orbital part 
H" = -{h'^/2m)A + V{r - Ri + Sj), which satisfies the 
eigenvalue equation (r-f-^j—R/) = e^^'>(l)^^\r+5j — 

R;) and a "perturbative part" AV which takes into ac- 
count the potential term that arises from all other atoms 
different from that in the atomic orbital Hamiltonian. 
The last line in Eq. pT|) has been obtained from the fact 
that the atomic wave functions 0'*^ (r) are eigcnstates of 
the atomic Hamiltonian with the atomic energy e*-'^ 
for an orbital of type i. This atomic energy plays the 
role of an onsite energy. The secular equation now reads 
det[4^ - (ck - e^''^)Sk ] = 0- Notice that, if the the atoms 
on the different sublattices are all of the same electronic 
configuration, one has e^'^ = eo for all i, and one may 
omit this on-site energy, which yields only a constant 
and physically irrelevant shift of the energy bands. 



b. Solution for graphene with nearest-neighbor and next- 
nearest-neighour hopping. After these formal considera- 
tions, we now study the particular case of the tight- 
binding model on the honeycomb lattice, which yields, 
to great accuracy, the tt energy bands of graphene. Be- 
cause all atomic orbitals are Pz orbitals of carbon atoms, 
we may omit the onsite energy cq, as discussed in the last 
paragraph. We choose the Bravais lattice vectors to be 
those of the A sublattice, i.e. Sa = 0, and the equivalent 
site on the B sublattice is obtained by the displacement 
= Sab = ^3 (sec Fig. The nn hopping amplitude 
is given by the expression 

t = J d\c^'^*{r)AV^''{r + S3), (14) 

and we also take into account next-nearest neighbor 
(nnn) hopping which connects neighboring sites on the 
same sublattice 

t„„„ = yciV(/.^*(r)AT/0^(r + ai) (15) 

Notice that one may have chosen any other vector 5j or 
a2, respectively, in the calculation of the hopping ampli- 
tudes. Because of the normalization of the atomic wave 
functions, we have / d'^rcf)^^^* {r)(j)^^'^ {r) = 1, and we con- 
sider furthermore the overlap correction between orbitals 
on nn sites, 

s = J d\cj)'^*{r)(t>^ir + 63). (16) 

We neglect overlap corrections between all other orbitals 
which are not nn, as well as hopping amplitudes for larger 
distances than nnn. 

If we now consider an arbitrary site A on the A sub- 
lattice (Fig. HI), we may sec that the off-diagonal terms 
of the hopping matrix (|13p consist of three terms corre- 
sponding to the nn Bi, B2, and , all of which have the 
same hopping amplitude t. However, only the site B3 is 
described by the same lattice vector (shifted by ^3) as the 
site A and thus yields a zero phase to the hopping matrix. 
The sites Bi and B2 correspond to lattice vectors shifted 
by a2 and a.^ = a2 — ai, respectively. Therefore, they 
contribute a phase factor exp(«k • a2) and exp(ik ■ aa), 
respectively. The off-diagonal elements of the hopping 
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matrix may then be written as = tj^ 
well as those of the overlap matrix s^^ 



it 



is 



AA 



^BB 



57k 



BAy 
„BA 



as 

r, 

1 , due to the above-mentioned normaliza- 
tion of the atomic wave functions), where we have defined 
the sum of the nn phase factors 



7k = 1 



jka2 



(17) 



The nnn hopping amplitudes yield the diagonal elements 
of the hopping matrix, 



^AA _ ,BB _ r,. 



3 

nnn ^ ^ COs(k • a^) 

i=l 



(l7k| 



and one obtains, thus, the secular equation 



det 



^-AA 
^k 



Ck it - sek)7k 



it - sek)7k t 



AA 



Ek 



with the two solutions (A ~ ±) 



j-AA 
''k 



At|7k| 



1 + As|7k| 



(18) 
(19) 

(20) 



This expression may be expanded under the reasonable 
assumptions s <C 1 and i„„„ ^ t, which we further jus- 
tify at the end of the paragraph. 



tt'' + 



^ t' 



At|7k| - st\"fk\ 

3 

3 + 2 ^ cos(k • Sii) 



4nnl7kr +At|7k| 



-At 



3 -I- 2^ cos(k • a^), 



(21) 



where we have defined the effective nnn hopping ampli- 
tude t'„„„ = tnnn — st, and we have omitted the unim- 
portant constant — 3i„„„ in the second step. Therefore, 
the overlap corrections simply yield a renormalization 
of the nnn hopping amplitudes. The hopping ampli- 
tudes may be determined by fitting the energy dispersion 
([21]l obtained within the tight-binding approximation to 
those calculated nume rically in more sophisticated band- 
structure calculations ( Partoens and Peeterd. 2006) or to 
spect roscopic measurements ( Mucha-Kruczvhski et all . 
120081 ). These yield a value of t ~ —3 eV for the nn 
hopping amplitude and fj^^n — ^-U, which justifies the 
above-mentioned expansion for t'^^^jt ^ 1. Notice that 
this fitting procedure does not allow for a distinction be- 
tween the "true" nnn hopping amplitude tnnn and the 
contribution from the overlap correction ~st. We, there- 
fore, omit this distinction in the following discussion and 
drop the prime at the effective nnn hopping amplitude, 
but one should keep in mind that it is an effective pa- 
rameter with a contribution from nn overlap corrections. 



^ The hopping matrix element corresponds to a hopping from 
the B to the A sublattice. 




Figure 5 Energy dispersion as a function of the wave- vector com- 
ponents fc^ and fcy, obtained within the tight- binding approxima- 
tion, for t„nn/t = 0.1. One distinguishes the valence (tt) band from 
the conduction (tt*) band. The Fermi level is situated at the points 
where the tt band touches the n* band. The energy is measured in 
units of t and the wave vector in units of I /a. 



c. Energy dispersion of tt electrons in graphene. The en- 
ergy dispersion pip is plotted in Fig. [5]for tnnn/t = 0.1. 
It consists of two bands, labeled by the index A = ±, 
each of which contains the same number of states. Be- 
cause each carbon atom contributes one tt electron and 
each electron may occupy either a spin-up or a spin- 
down state, the lower band with A = — (the tt or valence 
band) is completely filled and that with A = -I- (the tt* 
or conduction band) completely empty. The Fermi level 
is, therefore, situated at the points, called Dirac points, 
where the tt band touches the tt* band. Notice that only 
if ^nnn = the energy dispersion ([2T|) is electron-hole 
symmetric, i.e. = — ej^^. This means that nnn hop- 
ping and nn overlap corrections break the electron-hole 
symmetry. The Dirac points are situated at the points 
k''^ where the energy dispersion (j2ip is zero. 



<=ko 



0. 



(22) 



Eq. ([22]) is satisfied when J]^d = 0, i.e. when 



Rc7k 



1 + cos 



-|- COS 



\/3a 
V3a 



(-fcf + Vsk^ 



(23) 



and, equally, 
Im7j.r> 



(fcf + Vsk^) 



^i-k^ + V3k^) 



0. (24) 
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The last equation may be satisfied by the choice 
and Eq. ([23|) . thus, when 



0, 



1 + 2 cos 



= 



47r 
3 a/3 a 



(25) 

Comparison with Eq. Q shows that there are, thus, two 
inequivalent Dirac points D and D' , which are situated 
at the points K and K' , respectively. 



An 

±K = ± — 

3\/3a 



(26) 



Although situated at the same position in the first BZ, it 
is useful to make a clear conceptual distinction between 
the Dirac points D and D' , which are defined as the 
contact points between the two bands tt and tt* , and the 
crystallographic points K and K' , which are defined as 
the corners of the first BZ. There are, indeed, situations 
where the Dirac points move away from the points K and 
K' , as we will discuss in Sec. II.DI 

Notice that the band Hamiltonian ^ respects time- 
reversal symmetry, = ^Ik' "which implies e_k = ek 
for the dispersion relation. Therefore, if is a solution 
of ek = 0, so is — k"^, and Dirac points thus necessarily 
occur in pairs. In graphene, there is one pair of Dirac 
points, and the zero-energy states are, therefore, doubly 
degenerate. One speaks of a twofold valley degeneracy, 
which survives when we consider low-energy electronic 
excitations that are restricted to the vicinity of the Dirac 
points, as is discussed in Sec. II.C.2I 



d. Effective tight-binding Hamiltonian. Before considering 
the low-energy excitations and the continuum limit, it is 
useful to define an effective tight-binding Hamiltonian, 



Here, 1 represents the 2x2 one-matrix 

1 



1 = 



1 



(27) 



(28) 



This Hamiltonian effectively omits the problem of non- 
orthogonality of the wave functions by a simple renor- 
malization of the nnn hopping amplitude, as alluded to 
above. It is therefore simpler to treat than the origi- 
nal one ([8|) the eigenvalue equation of which involves the 
overlap matrix iSk, while it yields the same dispersion re- 
lation (|21l) . The eigenstates of the effective Hamiltonian 
([27]) are the spinors 



(29) 



the components of which arc the probability amplitudes 
of the Bloch wave function ([5]) on the two different sublat- 
ticcs A and B. They may be determined by considering 



Ai|7k|*^ 



the eigenvalue equation 'Hkitnnn — "7^k ~ ^i-r/ki^ki 
which does not take into account the nnn hopping cor- 
rection. Indeed, these eigenstates are also those of the 
Hamiltonian with 7^ because the nnn term is pro- 
portional to the one- matrix 1. The solution of the eigen- 
value equation (|29p yields 



I7k| 



and, thus, the eigenstates 



(30) 



1 

V2 



1 

Xe'f^ 



(31) 



where <y9k ~ arctan(Im7k/Rc7k). 

As one may have expected, the spinor represents an 
equal probability to find an electron in the state VP^ on 
the A as on the B sublattice because both sublattices are 
built from carbon atoms with the same onsite energy e*^*). 



2. Continuum limit 

In order to describe the low-energy excitations, i.e. 
electronic excitations with an energy that is much smaller 
than the band width ~ \t\, one may restrict the exci- 
tations to quantum states in the vicinity of the Dirac 
points and expand the energy dispersion around ±K. 
The wave vector is, thus, decomposed as k = ±K -t- q, 
where |q| <C |K| ~ 1/a. The small parameter, which gov- 
erns the expansion of the energy dispersion, is therefore 
|q|a < 1. 

It is evident from the form of the energy dispersion 
([21]) and the effective Hamiltonian that the basic entity 
to be expanded is the sum of the phase factors 7k. As 
we have already mentioned, there is some arbitrariness 
in the definition of 7k, as a consequence of the arbitrary 
choice of the relative phase between the two sublattice 
components - indeed, a change 7k — > 7kexp(j/k) in Eq. 
(|17p for a real and non-singular function /k does not af- 
fect the dispersion relation ([^T|) . which only depends on 
the modulus of the phase-factor sum. For the series ex- 
pansion, it turns out to be more convenient not to use the 
expression ([T7|) . but one with /k = k ■ 8?,, which renders 
the e xpression more symmetric (|Bena and Montambauxl . 
[200i), 



e"^-^^7k 



ik-^i I jk-^2 



(32) 



In the series expansion, we need to distinguish further- 
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more the sum at the K point from that at the K' point, 

3 



7." 



i=i 



„±i27r/3 



1 + iq- ^1 - -(q- ^1 



1 + iq • ^2 - 7T(q • ^2 



l + ,q.^3__(q.^3)2 



(33) 



By definition of the Dirac points and their position at the 
BZ corners K and K' , we have 7q^°'' = 7±k = 0. We 
Umit the expansion to second order in |q|a. 



a. First order in |q|a. The first-order term is given by 



'q 



a 
'2 



^±i27r/3 



qy)e 



=F427r/3 



3a , 



-iqya = T — {qx ±iqy), 



(34) 



which is obtained with the help of sin(±27r/3) = ±\/3/2 
and cos(±27r/3) = —1/2. This yields the effective low- 
energy Hamiltonian 



nf-^ =^hvF{q,a- +^qyay), 

where we have defined the Fermi velocity"* 

Ma 3\t\a 



Vf 



2ft 



and used the Pauli matrices 
^ Q 1 and 



2ft 



t'V - 



~i 

1 



(35) 



(36) 



(37) 



Furthermore, we have introduced the valley pseudospin 
^ = ±, where ^ = + denotes the K point at +K and 
^ — ~ the K' point at — K modulo a reciprocal lattice 
vector. The low-energy Hamiltonian ([35|) does not take 
into account nnn-hopping corrections, which are propor- 
tional to |7kP and, thus, occur only in the second-order 
expansion of the energy dispersion [at order Odqja)^]. 
The energy dispersion ([2T|) therefore reads 



(38) 



independent of the valley pseudospin ^. We have already 
alluded to this twofold valley degeneracy in Sec. ILC-H in 



* The minus sign in the definition is added to render the Fermi 
velocity positive because the hopping parameter t ~ — 3 eV hap- 
pens to be negative, as mentioned in the last section. 



the framework of the discussion of the zero-energy states 
at the BZ corners. From Eq. (|38p it is apparent that 
the continuum limit |q|a <C 1 coincides with the limit 
|e| <C \t\, as described above, because |eq| = 3ia|q|/2 <^ 
\t\ then. 

It is convenient to swap the spinor components at the 
K' point (for ^ = -), 



(39) 



i.e. to invert the role of the two sublattices. In this case, 
the effective low-energy Hamiltonian may be represented 
as 



off ,4 



: £,hvF{qxO-'' + qy<^") = fti'FT^ ® q • <T, (40) 



i.e. as two copies of the 2D Dirac Hamiltonian H]j — 
vfp ■ (T (with the momentum p = ftq), where we have 
introduced the four-spinor representation 



(41) 



in the last line via the 4x4 matrices 



(T 

-cr 



(42) 



and cr = (cr^, cr^). In this four-spinor representation, the 
first two components represent the lattice components at 
the K point and the last two components those at the 
K' point. We emphasise that one must clearly distin- 
guish both types of pseudospin: (a) the sublattice pseu- 
dospin is represented by the Pauli matrices , where 
"spin up" corresponds to the component on one sublat- 
tice and "spin down" to that on the other one. A rotation 
within the SU(2) sublattice-pseudospin space yields the 
band indices A = ±, and the band index is. thus, in- 
timitely related to the sublattice pseudospin. (b) The 
valley pseudospin, which is described by a second set of 
Pauli matrices , the z-component of which appears in 
the Hamiltonian ([40|) . is due to the twofold valley degen- 
eracy and is only indirectly related to the presence of two 
sublattices. 

The eigenstates of the Hamiltonian (|40p are the four- 
spinors 



/ 1 




1 

71 



/ 



1 



where we have, now, 



(fid = arctan 



Qy 



(43) 



(44) 
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valence 
band 
(1 = -) 




momentum 

Figure 6 Relation between band index A, valley pscudospin ^, and 
ciiirality rj in grapiiene. 



b. Chirality. In high-energy physics, one defines the he- 
licity of a particle as the projection of its spin onto the 
direction of propagation ( Weinb er^ 1 1 9 9 5() . 



q - cr 



(45) 



which is a Hermitian and unitary operator with the eigen- 
values T] = ±, riq\ri = ±) = zL\ri = ±). Notice that 
cr describes, in this case, the true physical spin of the 
particle. In the absence of a mass term, the helicity op- 
erator commutes with the Dirac Hamiltonian, and the 
helicity is, therefore, a good quantum number, e.g. in 
the description of neutrinos, which have approximately 
zero mass. One finds indeed, in nature, that all neutri- 
nos are "left-handed" (77 = — ), i.e. their spin is antipar- 
allel to their momentum, whereas all anti-neutrinos are 
"right-handed" (?/ = +). 

For massive Dirac particles, the helicity operator ((45|) 
no longer commutes with the Hamiltonian. One may, 
however, decompose a quantum state j^*) describing a 
massive Dirac particle into its chiral components, with 
the help of the projectors 



I*/ 



1 



and 



1 



(46) 



In the case of massless Dirac particles, with a well-defined 
helicity \^!) = |r/ = ±), one simply finds 



I* 



and 



1 - ^qi 



0, 



^q. 



'^q 



-) =0, 



such that one may then identify helicity and chirality. 
Because we are concerned with massless particles in the 
context of graphene, we make this identification in the 
remainder of this review and use the term chirality. 

For the case of graphene, one may use the same def- 
inition (|45p . but the Pauli matrices define now the sub- 
lattice pscudospin instead of the true spin. The operator 



Tjq clearly commutes with the massless 2D Dirac Hamil- 
tonian (|40|) . and one may even express the latter as 



"-q 



(49) 



which takes into account the two-fold valley degeneracy, 
in terms of the valley pscudospin ^ = ± . The band index 
A, which describes the valence and the conduction band, 
is therefore entirely determined by the chirality and the 
valley pscudospin, and one finds 



(50) 



which is depicted in Fig. [SI 

We notice finally that the chirality is a preserved quan- 
tum number in elastic scattering processes induced by 
impurity potentials Vimp = V^(r)l that vary smoothly on 
the lattice scale. In this case, inter- valley scattering is 
suppressed, and the chirality thus conserved, as a conse- 
quence of Eq. ([50]) . This effect gives rise to the absen ce of 
backscattering in graphene (iShon and Andol . Il998l ) and 
is at the origin of Klein tunneling according to which a 
massless Dirac particle is fully transmitted, under nor- 
mal incidence, thro ugh a high electrostatic barrier with- 
out being reflected (jKatsnelson et al This rather 
counter-intuitive result was first considered as a paradox 
and led to the fo rmula t ion o f a charged vacuum in the 
potential barrier ( Kleinl . [19291 ) . which may be indentified 
in the framework of band theory with a Fermi level in 
the valence band. 



c. Higher orders in \q\a. Although most of the fundamen- 
tal properties of graphene are captured within the effec- 
tive model obtained at first order in the expansion of the 
energy dispersion, it is useful to take into account second- 
order terms. These corrections include nnn hopping cor- 
rections and off-diagonal second-order contributions from 
the expansion of 7k. The latter yield the so-called trigo- 
nal warping^ which consists of an anisotropy in the energy 
dispersion around the Dirac points. 

The diagonal second-order term, which stems from the 
nnn hopping, is readily obtained from Eq. ([M)) . 



9a' 



(47) 

independent of the valley index ^. 



(51) 



The off-diagonal second-order terms are <7q 
— hvpa(c[x — iS,qy)' /A. Notice that there is a natural 
energy hierarchy between the diagonal and off-diagonal 
second-order terms when compared to the leading lin- 
ear term; whereas the off-diagonal terms are on the or- 
der 0(|q|a) as compared to the energy scale /I'yj-lql, the 
diagonal term is on the order 0((i„„„/i)|q|a) and thus 
roughly an order of magnitude smaller. We therefore 
take into account also the off-diagonal third order term 
t7q*''^' = — ^/iWFa^(g^-f i^(7j^)|qp/8, which also needs to be 
considered when calculating the high-energy corrections 
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of the energy levels in a magnetic field (see Sec. III.B[) . 
Up to third order, the off-diagonal terms therefore read 



(52) 



where one may omit the valley-dependent sign before the 
y-components of the wave vector by sweeping the sublat- 
tice components in the spinors when changing the valley. 

In order to appreciate the influence of the second-order 
off-diagonal terms on the energy bands, we need to cal- 
culate the modulus of 7^, 



3a , 



l-e^cos(3^q) 



(53) 



where we have used the parametrization = |q| cosLp^ 
and qy = |q|sin(/3q, and where we have restricted the 
expansion to second order. Finally, the energy dispersion 
([21]) expanded to second order in |q|a reads 



9a' 



Afi.Wi?|q| 



1 



^— cos(3cpq) 



"(54) 

As mentioned in Sec. ll.C.li it is apparent from Eq. 
([M]) that the nnn correction breaks the electron-hole 

This is, however, a rather 
|q|at„„„/t, to the first-order 



-A 



symmetry ^ 
small correction, of order 

effective Hamiltonian (|40|) . The second-order expansion 
of the phase factor sum 7q yields a more relevant cor- 
rection ~ the third term in Eq. ([M]). that is of order 
|q|a ^ \q\atnnn/t ~ to the linear theory. It depends ex- 
plicitly on the valley pseudospin ^ and renders the energy 
dispersion anisotropic in q around the K and K' point. 
The tripling of the period, due to the term cos(3iy9q), is 
a consequence of the symmetry of the underlying lattice 
and is precisely the origin of trigonal warping. 

The trigonal warping of the dispersion relation is visu- 
alized in Fig. [71 where we have plotted the contours of 
constant (positive) energy in Fourier space. The closed 
energy contours around the K and K' points at low en- 
ergy are separated by the high-energy contoms around 
the r point by the dashed lines in Fig. [71 (a) at en- 
ergy \t + tnnnl the crossing points of which correspond 
to the M points. As mentioned above, the dispersion 
relation has saddle points at these points at the border 
of the first BZ, which yield van Hove singularities in the 
density of states. In Fig. [71 (b), we compare constant- 
energy contours of the full dispersion relation to those 
obtained from Eq. ([54l) calculated within a second-order 
expansion. The contours are indistinguishable for an en- 
ergy of e = 1^1/3 ~ 1 eV, and the continuum limit yields 
rather accurate results up to energies as large as 2 eV. 
Notice that, in today's exfoliated graphene samples on 
Si02 substrates, one may probe, by field-effect doping 
of the graphene sheet, energies which are on the order 
of 100 meV. Above these energies the capacitor breaks 



(a) ^ 

^^^^^ 
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1 
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-1 

A -1 
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Figure 7 Contours of constant (positive) energy in reciprocal 
space, (a) Contours obtained from the full dispersion relation II21I I. 
The dashed line corresponds to the energy t -\- tn„„, which sepa- 
rates closed orbits around the K and K' points (black lines, with 
energy e < t + tnnn) from those around the F point (gray line, with 
energy e > t + t„„n)- (b) Comparison of the contours at energy 
e = 1 eV, 1.5 eV, and 2 eV around the K' point. The black lines 
correspond to the energies calculated from the full dispersion re- 
lation 1 I2III and the gray ones to those calculated to second order 
within the continuum limit l|54|l . 




axis of 
deformation 



Figure 8 Quinoid-type deformation of the honeycomb lattice - 
the bonds parallel to the deformation axis (double arrow) are mod- 
ified. The shaded region indicates the unit cell of the oblique 
lattice, spanned by the lattice vectors ai and a2. Dashed and 
dashed-dotted lines indicate next-nearest neighbors, with charac- 
teristic hopping integrals tnnn and t'„„„, respectively, which are 
different due to the lattice deformation. 



down, and Fig. [71 (a) indicates that the continuum limit 
([5^ yields extremely accurate results at these energies. 

We finally mention that, when higher-order terms in 
|q|a are taken into account, the chirality operator ([45l) 
no longer commutes with the Hamiltonian. Chirality is 
therefore only a good quantum number in the vicinity of 
the Dirac points. 



D. Deformed Graphene 

In the previous section, we have considered a per- 
fect honeycomb lattice which is invariant under a 27r/3 
rotation. As a consequence, all hopping parameters 
along the nn bonds 5j were equal. An interesting 
situation arises when the graphene sheet is deformed, 
such that rotation symmetry is broken. In order to 
illustrate the consequences, we may apply a uniax- 
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Figure 9 Band dispersion of the quinoid-type deformed the hon- 
eycomb lattice, for a lattice distortion of Sa/a = —0.4, with t = 3 
eV, t„„n/t = 0.1, dt/da = -5 eV/A, and dtnnn/da = -0.7 eV/A. 
The inset shows a zoom on one of the Dirac points, D' . 



ial strain in the y-direction,^ 



in which case one obtains a quinoid-type deformation 
(Fig. El). The hoppina; t' along 63 is then differ- 



20071) . 
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, dt , 

t^t=t + —Sa. 

da 



(55) 



Furthermore, also four of six nnn hopping integrals are 
affected by the strain (see Fig. [5]), 



tn 



-^t' -t 



da 



■Sa. 



(56) 



If one considers a moderate deformation e = 5a/a <C 1, 
the effect on the hopping ampl itudes may be e stimated 
with the help of Harrison's law ( Harrisonl Il98l[ ). accord- 
ing to which t ~ Ch^ /ma^ , where C is a numerical pref- 
actor of order 1. One therefore finds a value 



dt_ 

da 



2t 

a 



-4.3eV/A and i' = <(1 - 2e) (57) 



which coincides well with the value dt/da c± 5 eV/A , 
which may be fou nd in the literature ( Dillon et aZ.l . ll977t 
ISaito etall . Il998h . The estimation of the modified nnn 
hopping integral t'„„„ is slightly more involved. One may 
use a law tnnn{b,a) sa t(a)exp[— (& — a)/d( a)] farniliar i n 
the context of the extended Hiickel model (ISaleml . [l966l) . 
where b is the nnn distance, and d ~ a/3.5 « 0.4 A is 
a caracteristic distance related to the overlap of atomic 
orbitals. In undeformed graphene one has b = a-\/3, 



^ In our simplified model, we only consider one bond length 
changed by the strain. T he more general cas e has been con- 
sidered by Peirera et al, l|Pereira et al However, the 
main effects are fully visible in the simplified model. 



whereas in quinoid-type graphene b' = b{l + e/2), which 
gives 



t'nnn = *nnn(l - 2e + be/2d). 



(58) 



The electronic properties of quinoid-type graphene 
may then be described in terms of an effective Hamil- 
tonian of the type (jST]) . 



"Hk — tnnnhk'^ + t 



with (|Goerbig et all l2008t) 



z cos a + Z 



f' 



7;: 

7k 



(59) 



+ kya[- + e 



■ cos 



^/Sk^a 



kya + £ 



(60) 



and the off-diagonal elements 



7k = 2e''=«"(3/2+e) I y^kA + (1 _ 2e). (61) 



The resulting energy dispersion 



■ Ai|7k| 



(62) 



is plotted in Fig. [9] for an unphysically large defor- 
mation, e = 0.4, for illustration reasons. Notice that 
the reversible deformations are limited by a value of 
e ^ 0.1 ■..0.2 bey ond which the graphene sheet cracks 
( Lee et al I I2OO8I) . One notices, in Fig. [9l two effects 
of the deformation: i) the Dirac points no longer coin- 
cide with the corners of the first BZ, the form of which 
is naturally also modified by the deformation; and ii) the 
cones in the vicinity of the Dirac points are tilted, i.e. the 
nnn hopping term ([60|l breaks the electron-hole symme- 
try already at linear order in |q|a. These two points are 
discussed in more detail in the following two subsections. 



1. Dirac point motion 

In order to evaluate quantitatively the position of the 
Dirac points, which are defined as the contact points be- 
tween the valence (A = — ) and the conduction (A = -I-) 
bands, one needs to solve the equation 7kD = 0, in anal- 
ogy with the case of undeformed graphene discussed in 
Sec. II.C.ll One then finds 



and 



,D .2 ft' 
k„ a = £ —pz arccos — — 



(63) 



where the valley index ^ = ± denotes again the two in- 
equivalent Dirac points D and D' , respectively. As al- 
ready mentioned, the Dirac points D and D' coincide, for 
undistortcd graphene, with the crystallographic points K 
and K', respectively, at the corners of the first BZ. The 
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A<0 



2A A>0 



Figure 10 Topological semi-metal insulator transition in the model 
I I64II driven by the gap parameter A. (a) Two well-separated Dirac 
cones for A <^ 0, as for graphene. (b) When lowering the modulus 
of the (negative) gap parameter, the Dirac points move towards a 
single point, (c) The two Dirac points merge into a single point 
at the transition (A = 0). The band dispersion remains linear in 
the gy-direction while it becomes parabolic in the g^i-direction. (d) 
Beyond the transition (A > 0), the (parabolic) bands are separated 
by a b and gap A (insulating phase) . From iMontambaux et al\ , 
l2009al . 



gj,-direction, as one would expec t for Dirac points , but it 
is quadratic in the g^i-direction ( Dietl "aZl l2008[ ). This 



is a general feature of merging points, which may only 
occur at the F point or else at half a reciprocal lattice 
vector G/2, i.e. in the center of a BZ borde r line (such as 
the M points) (jMontambaux et all . l2009a[ ) . Indeed, one 
may show that in the case of a time-reversal symmetric 
Hamiltonian, the Fermi velocity in the a;-dircction then 
vanishes such that one must take into account the 
quadratic order in in the energy band. Notice that 
such hybrid semi-Dirac points, with a linear-parabolic 
dispersion relation, are unaccessible in graphene be- 
cause unphysica l ly large stra i ns w ould be required 



cause unpnysica i iy large stra i ns w ouia De required 
(|Lee et all l2008l : iPereira et ail , [looi). However, such 



points may exist in other physical sys t ems such as cold 
atom s in optical latt ic es (iHou et all 20091: Lee et al 



20091: IWunsch et all \20m IZhao and Paramekanti . 



20061 : iZhu et all |2003), the quasi-2D organ ic ma- 
terial a;-(BEpT-TTF )2l.-i (|Katavama et all 120061: 
iKobavashi et all l2007l) or V02/Ti02 hcterostructures 



(iBaneriee et aLl . 120091 ) 



2. Tilted Dirac cones 



distortion makes both pairs of points move in the same 
direction due to the negative value of dt/da. However, 
unless the parameters are fine-tuned, this motion is dif- 
ferent, and the two pairs of points no longer coincide. 

One further notices that Eq. (|63p has (two) solutions 
only for t' < 2t. Indeed, the two Dirac points merge 
at the characteristic point M" at the border of the first 
BZ (see Fig. [3]). The point t' = 2t is special insofar 
as it characterizes a topological phase transition between 



cones and a 


band insulator (for t' > 2t) (iDietl et al. 


2008: Esaki et all 120091: IMontambaux et al. l2009allb 


Pereira et al. 


. \20m IWunsch et all 12008). In the vicin- 



ity of the transition, one may expa nd the Hamiltonian 
(1591) a round the merging point Al" I Montambaux et all 
l2009al lbl). and one finds^ 



njM _ 







2m* 



ihCQy 







(64) 



in terms o f the mass m* = 2h^/3ta^ and the velocity 
c = 3ta/n (jMontambaux et a^.l . [2009bl) . The gap param- 
eter A = t' — 2t changes its sign at the transition - it is 
negative in the semi-metallic and positive in the insulat- 
ing phase, where it describes a true gap (Fig. [TO)) . 

The Hamiltonian ([64]) has quite a particular form in 
the vicinity of the merging points: it is linear in the 



We do not consider the diagonal part of the Hamiltonian, here, 
i.e. we choose tnnn = 0, because it does not affect the position 
of the Dirac points. 



Another aspect of quinoid-type deformed graphene and 
a consequence of the fact that the Dirac points no longer 
coincide with the BZ corners K and K' of high crystallo- 
graphic symmetry is the tilt of the Dirac cones. This 
may be appreciated when expanding the Hamiltonian 
((59)) to linear order around the Dirac points ^k^, in- 
stead of an expansion around the point M" as in the 
last subsection. In contrast to the undeformed case (|5ip. 
the d iagonal components fek now yield a linear contribu- 
tion ( Goerbig et a/.l . [20081 ) . tnrm/i^ko+ql — C^wg-ql, in 
terms of the tilt velocity 

2\/3 

wqx = —^{tnnnasm29 + t'„,^,^asm6) and woy = 0, 

(65) 

where we have defined 9 = arccos(— t'/2i). The linear 
model is therefore described by the Hamiltonian,^ 

+ WxQxCr'^ + WyqyCjy), (66) 

with the renormalized anisotropic velocities 

y/?>ta 



■ sin. 9 



and 



Wy 



3t'a / 2 

2TV+r 



^ This model may be viewed as the minimal form of the generalized 
Weyl Hamiltonian (with cr" = 1) 

M=0,...,3 

which is the most general 2x2 matrix Hamiltonian that yields 
a linear dispersion relation. 
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Diagonalizing the Hamiltonian 
sion relation 



66|) yields the disper- 



observe the effect, e.g. in angle-resolved photoemis- 
sion spectroscopy (ARPES) measurements ( Damascellil . 
|2004[ ) tha t have been successfu lly applied to epitaxial 



(q) = Hwq ■ q 4- Xh^w^q^ + WyQy, (67) graphene (iBostwick et all l2007t ) and graphitic samples 

( Zhou et ai . 20061 ). Notice that the Dirac cones are 



and one notices that the first term {hwQ-q) breaks indeed 
the symmetry (q) = e^(— q) in each valley, i.e. it tilts 
the Dirac cones in the direction opposite to wq , as well as 
the electron- hole symmetry eA(q) = — e_A(q) at the same 
wave vector.^ Indeed, the linearity in q of the generalized 
Weyl Hamiltonian (|66)) satisfies only the symmetry = 

— "Hlq inside each valley. 

Furthermore, one notices that the chiral symme- 
try is preserved even in the presence of the tilt term 
if one redefines the chirality operator (|45p as rjq ~ 

{w^qxcr^ + Wyqyay)/ ^w'lql + w^ql, which naturally 

commutes with the Hamiltonian (|66p . The eigenstates 
of the chirality operator are still given by 



(68) 



with tan(^k = Wyqy/wxqx, and one notices that these 
states are also the natural eigenstates of the Hamiltonian 

One finally notices that not all values of the tilt pa- 
rameter Wq are indeed physical. In order to be able to 
associate A = -|- to a positive and A = — to a negative 
energy state, one must fulfill the condition 



Wo < 1, 
in terms of the tilt parameter 



Wo 



'wcteV / wpy \ ^ 



Wy J 



(69) 



(70) 



In the particular case of the deformation in the y-axis, 
which is discussed here and in which case Woy = [see 
Eq. (|55|) ]. the general form of the tilt parameter reduces 
to Wo = wox/wx- Unless this condition is fulfilled, the 
iso-energctic lines are no longer ellipses but hyperbolas. 
In quinoid-type defo rmed graphene , the tilt parameter 
may be evaluated as ( Goerbig et al\ . |2008| ) 



(71) 

where we have used Eqs. ((57|) and ((58|) . Even at mod- 
erate deformations (e < 0.1), the tilt of the Dirac cones 
is on the order of 5%, and one may therefore hope to 



In the absence of the tilt term Swq • ql, this is a consequence 
of the symmetry u^'H.fJ^ = — "H, which is satisfied both by the 
effective Hamiltonian l|27p for tnn-n = and the linearised version 
II40I I in each valley for undcformed graphene. 



naturally tilted in a-(BEDT -TTF)2l.s (jKatavama et all 
l2006l : lKobavashi et a/.l . [2007[ l ■ where the Dirac points oc- 
cur at positions of low crystallographic symmetry within 
the first BZ. 



II. DIRAC EQUATION IN A MAGNETIC FIELD AND 
THE RELATIVISTIC QUANTUM HALL EFFECT 

As already mentioned in the introduction, a key ex- 
periment in graphene research was the discovery of a 
particular quantum Hall effect (jNovoselov et al I l2005al: 
IZhang et al. which unveiled the relativistic nature 

of low-energy electrons in graphene. For a deeper under- 
standing of this effect and as a basis for the following 
parts, we discuss here relativistic massless 2D fermions 
in a strong quantizing magnetic field (Sec. III.Ap . The 
limits of the Dirac equation in the treatment of the high- 
field properties of graphene are discussed in Sec. III.Bl 
and we terminate this section with a discussion of the 
relativistic Landau level spectrum in the presence of an 
in-plane electric field (Sec. III.C|) and that of deformed 
graphene (Sec. III.D[) . 



A. Massless 2D Fermions in a Strong Magnetic Field 

In order to describe free electrons in a magnetic field, 
one needs to replace the canonical momentum p b y the 
gauge-invariant kinetic momentum (| JacksonL [l999t ) 



p ^ n = p + eA(r) 



(72) 



where A(r) is the vector potential that generates the 
magnetic field B = V x A(r). The kinetic momentum is 
proportional to the electron velocity v, which must natu- 
rally be gauge-invariant because it is a physical quantity. 

In the case of electrons on a lattice, the substitution 
(|72p . which is then called Peierls substitution, remains 
correct as long as the lattice spacing a is much smaller 
than the magnetic length 



(73) 



which is the fundamental length scale in the presence 
of a magnetic field. Because a = 0.24 nm and Ib — 
26 nm/ ^/Bff], this condition is fulfilled in graphene for 
the magnetic fields, which may be achieved in today's 
high-field laboratories (~ 45 T in the continuous regime 
and ~ 80 T in the pulsed regime). 

With the help of the (Peierls) substitution ([7^ . one 
may thus immediately write down the Hamiltonian for 
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charged particles in a magnetic field if one knows the 
Hamiltonian in the absence of the field, 



■H(p) ^ H{U) = H{p + eA) = i7^(p, r). 



(74) 



Notice that because of the spatial dependence of the 
vector potential, the resulting Hamiltonian is no longer 
translation invariant, and the (canonical) momentum 
p = fi,q is no longer a conserved quantity. For the Dirac 
Hamiltonian (|40p . which we have derived in the preced- 
ing section to lowest order in |q|a, the Peierls substitution 
yields 



off, 5 



(75) 

We further notice that, because electrons do not only 
possess a charge but also a spin, each energy level re- 
sulting from the diagonalization of the Hamiltonian (j75p 
is split into two spin branches separated by the Zeeman 
effect Az = gjisB, where g is the q-factor of the host 
material [g ~ 2 for graphene ( Zhang et al\ . |2006() ] and 
/is = eh/2mQ is the Bohr magneton, in terms of the 
bare electron mass rriQ. In the remainder of this sec- 
tion, we concentrate on the orbital degrees of freedom 
which yield the characteristic level structure of electrons 
in a magnetic field and therefore neglect the spin degree 
of freedom, i.e. we consider spinless fermions. Effects 
related to the internal degrees of freedom are discussed 
in a separate section (Sec. |V| in the framework of the 
quantum-Hall ferromagnet. 



1. Quantum-mechanical treatment 

One may easily treat the Hamiltonian (|75p quantum- 
mechanically with the help of the s t andar d canonical 
quantization (jCohen-Tannoudii et all Il973[) . according 
to which the components of the position r = (x, y) and 
the associated canonical momentum p = (j)x,Py) sat- 
isfy the commutation relations [a;,p^] = [y,Py] = ih and 
[x,y\ = [pxjPy] = [x.Py] = [v.Px] =0. As a conse- 
quence of these relations, the components of the kinetic 
momentum no longer commute, and, with the help o f the 
commutator relation ( Cohen- Tannoudii et a?.l . ll973l ) 



PlJ{02)] 



(76) 



between two arbitrary operators, the commutator of 
which is an operator that commutes itself with both Oi 
and ©2, one finds 



[n.,H,] 



ox 



dy 



in terms of the magnetic length ([75]) . 

For the quantum-mechanical solution of the Hamilto- 
nian (j75p , it is convenient to use the pair of conjugate op- 
erators Ha; and Ily to introduce ladder operators in the 
same manner as in the quantum-mechanical treatment 



of the one-dimensional harmonic oscillator. These lad- 
der operators play the role of a complex gauge-invariant 
momentum (or velocity), and they read 

a - (H, - tlly) and «^ - ^ (H. + *H,) , (78) 

where we have chosen the appropriate normalization such 
as to obtain the usual commutation relation 



[a,a^l = 1. 



(79) 



It turns out to be helpful for practical calculations to 
invert the expression for the ladder operators ([75]) . 



H, 



V2l 



^ 'a^ + d) and H.^ = (a^ - d) . (80) 



2. Relativistic Landau levels 



In terms of the ladder operators (|78|) . the Hamiltonian 
([75|) becomes 



Lb 



d 
at 



(81) 



One remarks the occurence of a characteristic frequency 
Lu' ~ y/2vp/lB, which plays the role of the cyclotron 
frequency in the relativistic case. Notice, however, that 
this frequency may not be written in the form eB/mi, 
because the band mass is strictly zero in graphene, such 
that the frequency would diverge.^ 

The eigenvalues and the eigenstates of the Hamilto- 
nian (jST]) are readily obtained by solving the eigenvalue 
equation Hgipn ~ e-n'tjjn, in terms of the 2-spinors, 



(82) 



We thus need to solve the system of equations 

S^hhj' dVn — en Un and S^hu)' d^ Un — Sn Vn, (83) 

which yields the equation 



(84) 



for the second spinor component. One may therefore 
identify, up to a numerical factor, the second spinor com- 
ponent Vn with the eigenstate |n) of the usual number 
operator d^d, with a^ajn) — n\n) in terms of the inte- 
ger n > 0. Furthermore, one observes that the square 
of the energy is proportional to this quantum number. 



® Sometimes, a density-dependent cyclotron mass mc is formally 
introduced via the equality u)' = eB/mQ. 



15 



(a) 



(b) 




-,n=4 
-,n=3 
-,n=2 



magnetic field B 





Figure 11 (a) Relativistic Landau levels as a function of the mag- 
netic field, (b) Semi-classical picture of cyclotron motion described 
by the cyclotron coordinate tj, where the charged particle turns 
around the guiding center R. The gray region depicts the uncer- 
tainty on the guiding center, as indicated by Eq. I I98II . 



= {huj')^n. This equation has two solutions, a posi- 
tive and a negative one, and one needs to introduce an- 
other quantum number A = ±, which labels the states of 
positive and negative energy, respectively. This quantum 
number plays the same role as the band index (A = + for 
the conduction and A = — for the valence band) in the 
zero-i?-ficld case discussed in the preceding section. One 
thus obtains the spectrum ( McCluri . Il95(: 



A 



1^ 



'2n 



(85) 



of relativistic Landau levels (LLs) that disperse as XyBn 
as a function of the magnetic field [see Fig. [TlT a)]. Notice 
that, as in the B = Q case, the level spectrum is two-fold 
valley-degenerate. 

Once we know the second spinor component, the first 
component is obtained from Eq. (|83p . which reads m„ oc 
avn ^ a-\n) ^ I't- — 1) because of the usual equations 



a' \n 



) = Vn+l\n + 1) and a\n) = ^/^\n - 1) (86) 



for the ladder operators, where the last equation is valid 
for n > 0. One then needs to distinguish the zero-energy 
LL {n = 0) from all other levels. Indeed, for n = 0, the 
first component is zero because 



a\n = 0) = 
In this case one obtains the spinor 





|n = 0) 



(87) 



(88) 



In all other cases {n ^ 0), one has positive and nega- 
tive energy solutions, which differ among each other by 
a relative sign in one of the components. A convenient 
representation of the associated spinors is given by 



V2 \ eA|n> 



(89) 



The particular form of the n = spinor ([88)) associated 
with zero-energy states merits a more detailed comment. 
One notices that only the second spinor component is 
non-zero. Remember that this component corresponds 
to the B sublattice in the if-valley (^ = -I-) and to the 
A sublattice in the i^'-valley {£, — —)- the valley pseu- 
dospin therefore coincides with the sublattice pseudospin, 
and the two sublattices are decoupled at zero energy. No- 
tice that this is also the case in the absence of a magnetic 
field, where the relation ([50]) between the chirality, the 
band index and the valley pseudospin is only valid at non- 
zero values of the wave vector, i.e. not exactly at zero 
energy. Indeed, the chirality can no longer be defined as 
the projection of the sublattice pseudospin on the direc- 
tion of propagation q/|q|, which is singular at q = 0. At 
zero energy, it is therefore useful to identify the chiral- 
ity with the valley pseudospin. Notice, however, that this 
particularity concerns, in the absence of a magnetic field, 
only a non-extensive number of states (only two) because 
of the vanishing density of states at zero energy, whereas 
the zero-energy LL n = is macroscopically degenerate, 
as discussed in the following paragraphs. 



a. LL degeneracy. A particular feature of LLs, both rel- 
ativistic and non-rclativistic ones, consists of their large 
degeneracy, which equals the number of flux quanta 
Nb — Ax B/{h/e) threading the 2D surface A occupied 
by the electron gas. From the classical point of view, 
this degeneracy is related to the existence of a constant 
of motion, namely the position of the guiding center, i.e. 
the center of the classical cyclotron motion. Indeed, due 
to translation invariance in a uniform magnetic field, the 
energy of an electron does not depend on the position of 
this guiding center. Translated to quantum mechanics, 
this means that the operator corresponding to this guid- 
ing center R = {X, Y) commutes with the Hamiltonian 
•H(p + eA). 

In order to understand how the LL degeneracy is re- 
lated to the guiding-center operator, we formally decom- 
pose the position operator 



r = R + 77 



(90) 



into its guiding center R and the cyclotron variable 
V ~ {Vx,Vv): depicted in Fig. [TlTb). Whereas the 
guiding center is a constant of motion, as mentioned 
above, the cyclotron variable describes the dynamics of 
the electron in a magnetic field and is, classically, the 
time-dependent component of the position. Indeed, the 
cyclotron variable is perpendicular to the electron's ve- 
locity and thus related to the kinetic momentum 11 by 



(91) 



— - and 77,, = , 

eB eB ' 

which, as a consequence of the commutation relations 
([TT]) . satisfy 



[^x,Uy] 
{eBf 



iP 



(92) 
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whereas they commute naturaUy with the guiding-center 
components X and Y . Equation (|92p thus induces the 
commutation relation 



(93) 



in order to satisfy [x, y\ = 0. 

These commutation relations indicate that the compo- 
nents of the guiding-center operator form a pair of conju- 
gate variables, and one may introduce, in the same man- 
ner as for the kinetic momentum operator 11, the ladder 
operators 



1 



-{X + iY) and P 



1 



-{X-iY), (94) 



which again satisfy the usual commutation relations 
[6, 6^] = 1 and which naturally commute with the Hamil- 
tonian. One may then introduce a number operator wb 
associated with these ladder operators, the eigcnstates of 
which satisfy the eigenvalue equation 



(95) 



One thus obtains a second quantum number, an integer 
TO > 0, which is necessary to describe the full quantum 
states in addition to the LL quantum number n, and the 
completed quantum states (|88p and ([89)1 then read 



n— O.m 



and 



r Xn.m 





= 0, to) 



1 I \n — 1, m) 
71 V ^Mn,m) 



(96) 



(97) 



respectively. 

One may furthermore use the commutation relation 
(|93)) for counting the number of states, i.e. the degener- 
acy, in each LL. Indeed, this relation indicates that one 
may not measure both components of the guiding cen- 
ter simultaneously, which is therefore smeared out over a 
surface 



AX AY = 2ttII, 



(98) 



as it is depicted in Fig. [TlT b). The result ((98|) for the 
surface occupied by a quantum state may be calculated 
rather simply if one chooses a particular gauge, such as 
the Landau or the symmetric gauge for the vector poten- 
tial, but its general derivation is rather involved ( Imrvl . 
Il997t ). This minimal surface plays the same role as the 
surface (action) h in phase space and therefore allows 
us to count the number of possible quantum states of a 
given (macroscopic) surface A, 



Nf 



A 



AXAY 



A 

2nl% 



riB X A, 



(99) 



where we have introduced the flux density 

1 



riB 



2TTli 



B 



(100) 



which is nothing other than the magnetic field measured 
in units of the flux quantum /i/e, as already mentioned 
above. The ratio between the electronic density rig; and 
this flux density then deflnes the filling factor 



nei 
riB 



~eB 



(101) 



which characterizes the filling of the different LLs. 



b. The relativistic quantum Hall effect. The integer 
quantum Hall effect (I QHE) in 2D electron systems 



Klitzing etUl . ll980D is a manifestation of the LL 



quantization and the macroscopic degeneracy (jlOO[) of 
each level, as well as of semi-classical electron localiza- 
tion due to the sample impurities. In a nutshell, this 
energy quantization yields a quantization of the Hall re- 
sistance 



Rh 



(102) 



where N = [i/] is the integer part of the filling factor 
pOip . while the longitudinal resistance vanishes. The 
resistance quantization reflects the presence of an incom- 
pressible quantum liquid with gapped single-particle and 
density excitations. In the case of the IQHE, at integer 
filling factors, the gap is simply given by the energy dif- 
ference between adjacent LLs which must be overcome 
by an electron that one adds to the system. Notice that 
if one takes into account the electron spin and a vanish- 
ing Zceman effect, the condition for the occurence of the 
IQHE is satisfied when both spin branches of the last LL 
n are completely filled, and one thus obtains the Hall- 
resistance quantization at the filling factors 



IQHE 



2n, 



(103) 



i.e. for even integers. Odd integers may principally be 
observed at higher magnetic fields when the Zeeman ef- 
fect becomes prominent, and the energy gap is then no 
longer given by the inter-LL spacing but by the Zeeman 
gap. This picture is naturally simplistic and needs to be 
modified if one takes into account electronic interactions 



Strictly speaking, the IQHE requires only the breaking of trans- 
lation invariance, which in a diffusive sample is due to impuri- 
ties. In a ballisti c sample, trans lation invariance is broken via 
the sample edges l lBiittikeil . Il992l ) . 

A simultaneous measurement of the Hall and the longitudinal re- 
sistance requires a particular geometry with at least four electric 
contacts [ for a r ecent review on the quantum Hall effect, see Ref. 
l lGoerbiel . 120091 ')]. 
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- their consequences, such as the fractional quantum Hall 
effect or ferromagnetic states arc discussed, in the con- 
text of graphene, in Sec. El 

The phenomenology of the relativistic quantum Hall 
effect (RQHE) in graphene is quite similar to that of 
the IQHE. Notice, however, that one is confronted not 
only with the two-fold spin degeneracy of electrons in 
graphene (in the absence of a strong Zccman effect), 
but also with the two-fold valley degeneracy due to the 
presence of the K and K' points in the first BZ, which 
govern the low-energy electronic properties. The filling 
factor therefore changes by steps of 4 between adjacent 
plateaus in the Hall resistance. Furthermore, the fill- 
ing factor (|10ip is defined in terms of the carrier density 
which vanishes at the Dirac point. This particle-hole 
symmetric situation naturally corresponds to a half-filled 
zero-energy LL n = 0, whereas all levels with A = — arc 
completely filled and all A = -I- levels are unoccupied. 
In the absence of a Zeeman effect and electronic inter- 
actions, there is thus no quantum Hall effect at ly ~ 0, 
and the condition of a completely filled (or empty) n = 
LL is found ior v = 2 {v ^ —2). As a consequence, the 
signature of the RQ HE is a Hall-resistance quantization 
at the filling facto rs ( Gusvnin and Sharapovll2005i 120061 : 
iPeres et ad I2006D 



,i?QH£ ^2(2n+l), 



(104) 



which needs to be contrasted to the series (|103p of the 
IQHE in non-relativistic 2D electron systems. The se- 
ries (|104p has indeed been obs erved in 2005 wi t hin the 
quantum Hall me asurements ( Novoselov et al\ . l2005at 



IZhang et al. which thus revealed the relativistic 

character of electrons in exfoliated graphene. More re- 
cently, the RQHE has been observe d also in ep i taxial 
graphene with moderate mobilit ies ( Jobst et al\ . l2010l : 
IShen et all l2009l : IWu et al\ . ^m^ . 



c. Experimental observation of relativistic Landau levels. 
The %/ Bn dispersion of relativistic LLs has been observed 
experimentally in transmission spectroscopy, where one 
shines monochromatic light on the sample and measures 
the intensity of the transmitted light. S uch experiments 
have been performed both on e pitaxial (ISadowski et all . 
I2006D and exfoliated graphene (j jiang et all \2007^^i ^. 

When the monochromatic light is in resonance with a 
dipole-allowed transition from the (partially) filled (A, n) 
to the (partially) unoccupied LL (A',n ± 1), it is ab- 
sorbed due to an electronic excitation between the two 
levels. Notice that, in a non-relativistic 2D electron gas, 
the only allowed dipolar transition is that from the last 
occupied LL n to the first unoccupied one n -I- 1. The 
transition energy is hiuc, independently of n, and one 
therefore observes a single absorption line (cyclotron res- 
onance) that is robust to electron-el ectron interactions, 
as a consequence of Kohn's theorem ( Kohnl . Il96ll) . 

In graphene, however, there are many more allowed 
transitions due to the presence of two electronic bands. 



ma \- 

E 
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Figure 1 2 Transmission spe ctroscopy on epitaxial multilayer 
graphene l lPlochocka et aEl . l20081 . The inset shows a representative 
transmission spectrum. The main figure represents the positions of 
the absorption lines as a function of the square-root of the mag- 
netic field. The dashed lines correspond to transitions calculated 
at linear order, in agreement with the Dirac equation, whereas one 
notices downward deviations in the high-energy limit. 



the conduction and the valence band, and the transitions 
have the energies 



1 



CV2^ 



(105) 



where A = + denotes an intraband and A = 

— an inter band transition ("Abergcl and FaPkoL l2007t 
llvengar et al.. . ,2007. : .Sadowski et al.. . 2006[k One there- 
fore obtains families of resonances the energy of which 
disperses as A„_a oc -s/B, as it has been observed in the 
experi ments [see Fig. I12[ whe re we show the results from 
Ref. (|Plochocka et all l2008t )]. Notice that the dashed 
lines in Fig. [12] are fits with a single fitting parameter 
(the Fermi velocity wf), which matches well all experi- 
mental points for different values of n in the low-energy 
regime. 

Moreover, the relativistic LLs have later been directly 
observed in scanning-tunn eling sp ectroscopy in graphene 
on a graphite subst rate^^ (iLi et al... 20 09al) as well as on 
epitaxial graphene ( Song et all 201(^ . 



B. Limits of the Dirac Equation in the Description of 
Graphene Landau Levels 

Transmission spectroscopy is an ideal tool for the study 
of the high-energy part of the LL spectrum when consid- 
ering the transitions (A = — , n) (A = 4-, n ± 1), for 
n ^ 1. As discussed in Sec. 11.0.2) one expects devi- 
ations [of order C(|qpa^)] from the linear dispersion in 
this limit. These deviations renormalize the energy of 
the LLs and thus the transition energies. 



With the help of the same techni que, relativistic LLs had before 
been identified even in graphite l lLi and Andreil . l2007t) . 
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In order to quantify the effect ( Ploefiocka et all . \200^ . 
we may use the Peierls substitution ([72|) and the expres- 
sions (f50| in the terms ([?T|) and ([5^ corresponding to 
the higher-order diagonal and off-diagonal band terms, 
respectively. This yields the Hamiltonian 



H 



h' h* 



h' 

where the diagonal elements read 



h' = hu 



I 3^n 



and the off-diagonal ones are 



. awi 



(106) 



(107) 



(108) 



Naturally, to lowest order in a/Is, one obtains again the 
Hamiltonian (jSip . The dimcnsionless parameters wi and 
W2 are artificially added to the expressions and play the 
role of fitting parameters in the comparison with experi- 
mental measurements, as will be discussed below. They 
measure the deviation from the tight-binding-model ex- 
pectation, wi = 11)2 = Notice furthermore that, since 
we arc interested in the n ^ 1 limit, we do not care 
about corrections related to the ordering of the ladder 
operators, such that we identify a'^^a^ ~ a^a'^'^ ~ (a'^a)'^ 
in the following parts. 

In the calculation of the LL spectrum, one may pro- 
ceed in the same manner as in Sec. III. A. 21 - the eigen- 
value equation ((M|) for the second spinor component now 
becomes 



(109) 



which is asymptotically correct in the large-n limit, where 
we may neglect the commutator [/i^, h'] on the right-hand 
side of the equation. The combination = h^h}^ is 
now interpreted as some fake Hamiltonian which needs 
to be diagonalizcd in order to obtain the modified LLs. 
Notice that n remains a good quantum number if one 
considers h' on the right-hand side of the eigenvalue equa- 
tion. The left-hand side consists of a term 



l\2 



a' a 



-\ 2 



(110) 



that contains powers of a) a and thus respects the quan- 
tum number n, but in addition it contains the trigonal- 
warping term 



Wl(fujj') 

2V2h 



-{a 



ot3 



(111) 



The commutator would yield relative corrections that are on 
the order of 1/n and a/Is as compared to the energy scale 
{tnnn/t){a/lg)n that dominates h' . 



which does not commute with a) a and which needs to be 
treated apart. If we neglect this trigonal- warping term 
for a moment, the LL energies are obtained from the 
quadratic equation 



{huj') 



Awl 



3- / 3^nnn^ 



^/2^/i 



(112) 

In order to account for the trigonal-warping term in 
the eigenvalue equation (jl09p . we may use a perturbative 
treatment, which is justified because of the small param- 
eter a/Is- There is no contribution at first order since 
{n\a^^^'^\n) = due to the orthogonality of the eigenstates 
(n|n') — 5n,n' ■ At second order, one obtains 



[huj' 



3n [1 + 0(1/7 



(113) 



which needs to be added to the right-hand side in Eq. 
(|112p . Interestingly, trigonal warping thus yields the 
same correction to the energies of the relativistic LLs as 
the third-order term in the expansion of the band disper- 
sion, although trigonal warping occurs at second order in 
the absence of a magnetic field, as we have discussed in 
Sec. II.C.2I This effect is due to the anisotropy of the 
band correction; in the presence of a magnetic field, the 
cos(3(pq) term in Eq. (I54p is averaged over the angle (/3q, 
and therefore only contributes at second order in the per- 
turbation theory described above. This eventually yields 
a correction of order (q/Ib)'^ to the LL energy, as docs the 
third-order term in the correction of the band dispersion. 

One finally obtains, in the large-n limit, where these 
correcti ons become relevant, th e energies of the relativis- 
tic LLs (jPlochocka et all |2008| ) 



Vp 3tr, 



t h 



(114) 



independent of the valley index ^, where ©(l/n) stands 
for corrections of order 1/n. Notice that the fitting pa- 
rameters wi and W2 cannot be determined independently 
from a fit to the spectroscopic measurement, but only the 
combination uP' = {wf + 2i(j|)/3. Equation (|114|1 gener- 
alizes a calculation for the relat ivistic LLs when o nly nnn 
hopping is taken into account (jPeres et al liooi). 

In Fig. (fT2l) . we show experimental results obtained 
from high-field trans mission spectroscopy on multi-layer 
epitaxial graphene ( Plochocka et al\ . 120081 ). Qualita- 
tively, one notices a downward renormalization of the 
transition energies 



(115) 



in the interband regime for large values of n, in agree- 
ment with Eq. 11141 Notice that because transmission 
spectroscopy is sensitive to energy-level differences, the 
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nnn correction in Eq. pi4p yields only a correction on 
the order of {tnnn/t){a/lB)/n < 1%/n at B - 25 T as 
compared to the energy scale t{a/lB)n of the transition, 
whereas the other term yields a correction on the order 
of {u/Ib)'^ xn ^ 0.5% X n. The latter corrections thus be- 
come more relevant in the large-n limit than the nnn cor - 
rection. Indeed, the experiment ( Plochocka et al\ . 120081 ) 
was not capable of probing the electron-hole symmetry 
breaking associated with the nnn term, whereas a quan- 
titative study of the high-energy transitions revealed a 
good semi-quantitative agreement with the calculated LL 
spectrum pi4p . However, it has been shown that the 
simple-minded tight-binding approach (with u; = 1) un- 
derestimates the higher-order band corrections and that 
the best fit to Eq. (|114p is obtained for a value of w = 2.8. 
The origin of this discrepancy is yet unexplained, and 
it may be interesting to perform high-field transmission 
spectroscopy measurement also on single-layer exfoliated 
graphene in order to understand whether the stronger 
downward renormalization of the LLs is due to interlayer 
couplings in the epitaxial multi-layer sample. 



such that the LLs (|85p . which scale as V B' 
{E /veBYY^^ with the magnetic field, read 



'2n, 



'B\l 



(117) 



where the primes indicate the physical quantities in the 
moving frame of reference. When measuring the energy 
in the original (lab) frame of reference, the above energy 
spectrum also needs to be transformed into this frame 
of reference, which amounts to being multiplied by an- 
other factor ()116p . such that the spectrum of relativistic 
LLs in the presence of an in-plane electric field becomes 
(iLukose et all |2007| ) 



e^,n;k = \^[l-{ElvpBff/S 



'2n 



h^k. (118) 



The quantum-mechanical derivation of this result will be 
discussed in more detail in Sec. IlI.D.ll in the context of 
the generalized Weyl Hamiltonian in a magnetic field. 



C. Landau Level Spectrum in the Presence of an Inplane 
Electric Field 

A remarkable consequence of the relativistic character 
of electrons in graphene and the Lorentz invariance of 
the Dirac equation is their behavior in crossed magnetic 
and electric fields, where the magnetic field remains per- 
pendicular to the graphene sheet and the electric field is 
applied in the plane. Remember that in a non-relativistic 
2D electron systems, the electric field E = Eey (in the 
y-direction) simply lifts the LL degeneracy and adds a 
term h{E/B)k to the LL energies, where k is the wave 
vector in the x-dircction. At a fixed wave vector fc, the 
LL spacing is unaffected by the in-plane field. 

The situation is different for relativistic electrons in 
graphene, as a consequence of the Lorentz invariance of 
the Dirac equation. One may choose a reference frame 
in which the electric field vanishes as long as the drift 
velocity vd = E/B is smaller than the Fermi velocity, 
which plays the role of an upper bound for the physically 
significant velocities in the same manner as the speed 
of light in relativity (magnetic regime).^* In addition 
to the lifted LL degeneracy, the LL spacing is reduced 
( Lukose et "^■ 120071: If eres and Castroll20Q7^ . as may be 
seen from a Lorentz boost into the reference frame which 
moves at the drift velocity and in which the electric field 
vanishes. In this reference frame, the magnetic field is 
reduced by the factor 



D. Landau Levels in Deformed Graphene 



As we have discussed in Sec. II. Dl a uniaxial strain 
deforms the graphene sheet and modifies the electronic 
structure. The induced anisotropy of the Fermi velocity 
Wx 7^ Wy is essentially washed out by the magnetic field, 
which yields an effective averaging over the Fermi surface, 
I'F v'p = ^WxWy. More spectacular are the two fur- 
ther consequences of the deformation; (a) the tilt of the 
Dirac cones accounted for in the generalized Weyl Hamil- 
tonian ([66|) and (b) the topological phase transition due 
to the Dirac point motion. The implication for the LL 
spectrum are briefly reviewed in the following sections. 



1. The generalized Weyl Hamiltonian in a magnetic field 



With the help of the Peierls substitution (|72|1 and the 
expression of the kinetic momentum in terms of ladder 
operators ((80|) . the generalized Weyl Hamiltonian ((66l) 
may be cast into the form 



where 



h^2wxWy f 2|i(ae*¥' + H.c.) 



fi(ae*'^ -KH.c.) 



Woe 



_ WQx ^ .WQy 

Wx Wy ' 



(119) 

(120) 



v/1 - {E/vFBf 



(116) 



In the opposite case, vu > vp, one may choose a referen ce system 
in wh ich the magnetic field vanishes (electric regime) l|Jacksonl . 
I1999D . 



in terms of the effective tilt parameter ([70|) and the angle 
<y5 between the a;-axis and the direction of the effective tilt 
{wQxlwx,WQy/wy), renormalized by the Fermi velocities 
Wx and Wy in the a;- and y-direction, respectively. 

The Hamiltonian (jll9p may be solved quantum- 
mechanically in a straight-forward, but length y manner 
(|Morinari et ad 120091: IPeres and Castroll2007l ). Instead, 
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one may al so obtain the re s ult in a simpler semi-classical 
treatment (IGoerbig et ai . 2008), with the help of the 



Onsager relation (jLifshitz and KosevichL Il956l : lOnsageil 
I1952I ) according to which the surface S{e) enclosed by 
a trajectory of constant energy e in reciprocal space is 
quantized as 



S{e)ll 



(121) 



where n is an integer denoting the energy level which co- 
incides with the Landau level in the full quantum treat- 
ment. The additional contribution 7 is related to a Berry 
phase acquired by an electron during its cyclotron orbit. 
Usually, one has 7 = 1/2 except if there is an extra Berry 
phase of tt, which in our case yields 7 = 0, as for graphene 
with no tilt ( Mikitik and Sharla 1. 119991) . If one considers 
a density of states which scales as p{e) oc e", the energy 
levels thus scale as 



[i3(n + 7)]i/(i+"), 



(122) 



in the large-n limit. 

Because the density of states vanishes linearly at the 
Dirac point, as in the case of no tilt, i.e. a = 1, the 
scaling argument (|122p yields the energy levels. 



eA,n - Av2- — y/n, 

IB 



(123) 



as for unconstrained graphene, apart from a renormaliza- 
tion of the Fermi velocity. The latter is readily obtained 
from the calculation of the total number of states below 
the energy e within the positive energy cone. 



N+ie) 



1 



dqxdqy 



1 



2T:tfvf 2 ' 



(124) 

where wc have defined qx/y = Wx/yqx/y, and the renor- 
malized Fermi velocity is 



(125) 



in terms of the effective tilt parameter ([70]) . This yields 
the result 



,„ = A^^(l-^^)3/4 



2n (126) 



for the LL spectrum in the presence of a tilt, which 
coincides with the one obtained from the full quan- 
tum treatment ([Morinari et m. 120091 : IP^es and Castrol 
I2OO7I ). One finally notices that the LL spacing becomes 
zero for wq = 1, which corresponds to the condition 
of maximal tilt for the Dirac cones, as discussed in Sec. 
II. Dl indeed for values of Wq larger than 1, the isoener- 
getic lines are no longer closed elliptic orbits but open 
hyperbolas, for which the energy is not quantised. 



2. Tilted Dirac cones in a crossed magnetic and electric field 

One notices that the form (|126[) of LLs for tilted Dirac 
cones is the same as that of the LL spectrum (|118p if 
one interprets the drift velocity vd = E/vpB as an ef- 
fective electric-field induced tilt. The magnetic regime 
E/B < vp corresponds then to the regime of closed or- 
bits (wq < 1) and the open hyperbolic orbits may be 
identified with the electric regime E/B > vp- Mathe- 
matically, the generalized Weyl Hamiltonian with an in- 
plane electric field may still be cast into the form (|119[) 



E 

h-ki, 



(127) 



where is the same as that of E g. (I119p if one replac es 
the tilt parameter wq exp{i(p) by (jGoerbig et all l2009l ) 



wdE)e'^^^''^ ^ ^ 



.W^y 



Here, the renormalizcd tilt velocity is given by 

E X B 



= iw^x,w^y) = Wo - C- 
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(128) 



(129) 



and the angle ip^ is the angle between this velocity and 
the X-axis. 

The resulting energy spectrum is given by 



hy/WxWy 



[1 ~ w^{E)^f^ V2^ 



h-k. 
B 



(130) 

Naturally one obtains the result (jllSp for undeformed 
graphene in an in-plane electric field, for Wx = Wy = vp 
and Wo = 0, as well as the LL spectrum (|126p for the 
generalized Weyl Hamiltonian with tilted Dirac cones for 
zero in-plane field (£' = 0). However, the most interest- 
ing situation arises when both the tilt and an in-plane 
field are present, in which case one observes a lifting of 
the valley degeneracy that is maximal when the electric 
field is applied perpen dicular to the tilt velocity, E _L wo 
(IGoerbig et aLl . l2009D . 

Notice that, in order to obtain an effect on the or- 
der of '--^ 1%, extremely large electric fields would be 
required (on the o rder of 10^ V/m) for a 10% deforma- 
tion of the lattice ( Goerbig et aLl . [2009( ). It seems there- 
fore difficult to observe the effect in graphene, e.g. in 
high-field transmission spectroscopy or transport mea- 
surments, whereas the effect may be more visible in 
Q— (BE pT-TTF)2l.ci, wh e re the Dirac cones are na t urally 
tilted (|Katavama et all 120061: iKobavashi et all 120071 ) 
and where lower electric fields would be required for a 
comparable effect due to a roughly ten times smaller ef- 
fective Fermi velocity. 



III. ELECTRONIC INTERACTIONS IN GRAPHENE - 
INTEGER QUANTUM HALL REGIME 

In the preceding sections, we have discussed the elec- 
tronic properties of graphene within a one-particle model, 
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i.e. we have neglected the Coulomb interaction between 
the electrons. In many materials, the one-particle picture 
yields the correct qualitative description of the electronic 
properties and is modified only quantitatively if one in- 
cludes the el ectron-electron interactions within perturba - 
tion theory (jGiuliani and Vignald . 120051: iMahanl . Il993( ). 
Notice, however, that there exists a class of materials 
- strongly correlated electron systems - the electronic 
properties of which may not be described correctly, not 
even on the qualitative level, within a one-particle pic- 
ture. 

In order to quantify the role of the electronic inter- 
actions, i.e. the correlations, in graphene, one needs to 



compare the characteristic Coulomb energy Ei, 



at the average inter-electronic distance (e is the dielectric 
constant describing the environment the 2D electron gas 
is embedded in) to the kinetic one Ekinikp) at the same 
length scale, given in terms of the Fermi wave vector kp, 



F ' 



Ei, 



Eki 



(131) 



If this dimensionless interaction parameter becomes very 
large, » 1, the electrons are strongly correlated. In 
non-relativistic 2D metals with a parabolic band dis- 
persion, Ekin ^ fi^kp/mb, the dimensionless parameter 
reads 



mf,e 



1 



a^kp ' 



(132) 



in terms of the effective Bohr radius Oq = aosmQ/mi,, 
where ag = 0.5 A is the Bohr radius in vacuum and 
mb /too the ratio between the band and the bare electron 
mass. The relevance of electronic correlations therefore 
increases in the dilute limit when ^ ^ Cg . Notice that the 
parameter rs, which is also called Wigner-Seitz radius, 
plays the role of a length measured in units of the effective 
Bohr radius Og. 

The same argument applied to graphene yields a com- 
pletely different result. Whereas the scaling of the 
Coulomb energy remains the same, that of the kinetic 
energy is changed due to the linearity of the band dis- 
persion. As a consequence the dimensionless interaction 
parameter in graphene reads 



Ejn 



Eki' 



hevf 



2.2 

e 



(133) 



independent of the carrier density.^^ The correlations 
are therefore in an intermediate regime but may be de- 
creased if the graphene sheet is embedded in an environ- 
ment with a large dielectric constant. Notice that the 



In contrast to an electron system with a parabolic band disper- 
sion, this parameter can no longer be interpreted as a dimension- 
less radius, and we therefore use the notation olq rather than r^. 



expression (|133p is the same as that of the fine structure 
con stant a = lh£c = 1/137 in quantum electrodynam- 



ics (IWeinberd . [l995l) if one replaces the Fermi velocity by 
the velocity of light, which is roughly 300 times larger. 
One therefore calls aa alternatively the graphene fine 
structure constant. 



Long-range versus short-range interactions. Another im- 
portant aspect of interacting electrons is the range of the 
interaction potential. Whereas the underlying Coulomb 
potential e^/er is long-range, short-range interaction 
models such as the Hubbard model are often - and suc- 
cessfully - used in the description of correlated met- 
als. The use of such a short-range interaction may be 
justified by the screening properties of interacting elec- 
trons, whi ch are correctly captured in a Thomas-Fe rmi 
approach I Giuliani and Vignald . 120051 : iMahaiil . Il993l ) ac- 
cording to which the Coulomb interaction potential is 
screened above a characteristic screening length Xtp ^ 
l/kxp.^^ In 2D, the Thomas-Fermi wave vector 



(134) 



is given in terms of the dimensionless interaction param- 
eter (|13ip and the Fermi wave vector kp.^"^ 

One notices that, for metals with a parabolic dis- 
persion relation, the Thomas-Fermi wave vector is sim- 
ply given in terms of the inverse effective Bohr radius, 
kpp ^ l/ogi independent of the electronic density. Un- 
less the band mass is very small as compared to the bare 
electron mass or the dielectric constant of the host ma- 
terial very large, the Coulomb interaction is therefore 
screened on the atomic scale. A description of such sys- 
tems in the framework of short-range interaction mod- 
els, such as the Hubbard model, then becomes better 
justified than in systems with a small band mass or a 
prominent dielectric constant (such as in 2D electron 
systems in GaAs heterostructures). Typical examples, 
where short-range interaction model yields valuable phys- 
ical insig ht, are heavy- fe rmion compounds [for a review 
see Ref. (IColemanl . [2003l) ]. 

The situation is again drastically different in graphene 
where the Thomas-Fermi wave vector (jl34|) becomes 



2.2 

aakp ~ — kp 



(135) 



i.e. it vanishes at the Dirac points where the carrier 
density goes to zero, and the screening length then di- 
verges. Notice that even for doped graphene, where 



^® Notice that the Thomas-Fermi approach is restricted to static 
screening effects, whereas dynamic screening require a more com- 
plex treatment, e.g. in the framework of the random-phase ap- 
proximation. 

In three space dimensions, the relation reads k^p ~ rgkp. 
Due to this divergence of the screening length, one principally 
needs t o describe screeni ng beyond the level of linear-response 
theory l lKatsnelsonl . 120061 ') . 
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one may typically induce carrier densities on the order 
of 10^^ cm"^, the screening length is Xtf ^ 10 nm. i.e. 
much larger than the lattice scale. 

One thus comes to the conclusion that the rele- 
vant electronic interactions in graphene are long-range 
Coulomb interactions that may not be captured, in con- 
trast to other materials with a parabolic ban d di spe rsion , 
within models such as the Hubbard model (iHerbuti . i2006i 
l2007aD . We therefore investigate, in this section, the fate 
of the long-range Coulomb interaction in a strong mag- 
netic field. In Sec. IIII.AI we decompose the Coulomb in- 
teraction Hamiltonian in the two-spinor basis of the low- 
energy electronic wave functions in graphene and com- 
ment on its symmetry with respect to the valley pseu- 
dospin. The role of these interactions in the particle-hole 
excitation spectrum is studied in Sec. IIII.Bl where we 
discuss the resulting collective excitations in the IQHE 
regime, which allows for a perturbative treatment. The 
strong-correlation regime of partially filled LLs (regime 
of the fractional quantum Hall effect) is presented sepa- 
rately in Sec. |Vl 

A. Decomposition of the Coulomb interaction in the 
Two-Spinor Basis 

Generally, the Coulomb interaction for 2D electrons 
may be accounted for by the Hamiltonian 

Htnt = \Y1 (136) 
q 

in terms of the Fourier components p(q) = 
/ c?^rexp(— iq • r)^^ {r)'ip{r) of the electronic den- 
sity 4'''' {r)i(j{r) and the 2D Fourier transformed l/r 
Coulomb potential, v{q) — 27re^/e|q|. If one takes 
into account the electronic spin a =t, i, the Coulomb 
interaction respects the associated SU(2) symmetry, 
and the Fourier components are then simply the 
sum of the densities Pcr(q) in both spin orientations. 



p(q) = p^(q) -I- /0^(q)- For notational convenience, we 
neglect the spin index in the following discussion keeping 
in mind that the spin SU(2) symmetry is respected. The 
density operators may be decomposed in the basis of 
the spinor wave functions ([96)1 and (|97)) for relativistic 
electrons in graphene, 

X'n' ,ra' 

(137) 

where c^^J ^^.^ are fcrmion operators in second quantiza- 
tion that annihilate (create) an electron in the quantum 
state 

and = ( i^^f^^i;^^^ (138) 

In order to avoid confusion in the case of inter-valley 
coupling, we use now a representation in which the first 
spinor component represents the amplitude on the A 
sublattice and the second on the B sublattice for both 
valleys. Contrary to the expressions and (|97p . the 
state ()138|) is valid for both n = and n 7^ by us- 
ing the short-hand notation 1* = ■\/(l — (5„^o)/2 and 
2* = 1^(1-1- (5„,o)/2. Furthermore, we have explicitly 
taken into account the rapidly oscillating part exp(i^K-r) 
due to the two different valleys, whereas the expressions 
(|96)l and (jOT)) are only concerned with the slowly varying 
envelope function. Explicitly, the Fourier components of 
the density operator (|137p then read 

p(q)= E •^ln,v„'(q)pi;CA'„'(q): (i39) 

An, A n 

in terms of the reduced density operators 



pi'^,'A'„'(q) = E 



-»[q+(«-«')K]-R 



/ "^An,™;?"^^'"',™';?'' 



(140) 



which may also be interpreted as magneto- exciton operators associated with a particular inter-LL transition (see Sec. 
lIV.Ap , and the /orm factors 

J-,„,v„'(q) ^ -^L«v„'(q) = 1:1:' {n - 1 le-^'i-'^l n' - l) + AA'2:2:, {n |e-'^i-^| n') (141) 
for intra-valley and 



•^A„,A'„'(q)=Ai;'2; n 



-i(q+2K)-?7 



n'-l)-A'i:,2!, (n-1 



-i(q+2K)-?7 



^y-Xxni-<l) (142) 



for inter-valley processes. Here, we have used the decomposition r = R -f r; of the position operator into its guiding 
center and cyclotron coordinate (see Sec. III.A.ip and the fact that /i(r7)/2(R)|n, m) = fi{ri)\n) (E) /2(R)|to), for two 
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arbitrary functions /i and f2- The full expressions for the matrix elements in Eqs. (|140p . (|141|) . and p42|) may be 
found in Appendix 1X1 

In terms of the reduced density operators ()140p . the interaction Hamiltonian (|136p reads 

H^nt = i E E ^'M;f.tA.„.(q)pt;tA3n3(-q)pti:,A.„.(q). (i43) 

q Ai?ii . ..X4n4 



v,.-!;, (b) . v,,-!; 



(C)\v,,-1; 




Vj.l; (d)\V5,l; 




V4.-!; 




Figure 13 Diagrammatic representation of the interaction ver- 
tex (we use the short-hand notation Ui = {\ini,m,i) for the 
quantum numbers); (a) vertex associated with terms of the form 



^4(1) °^ ■"Aim' ■■^4714(1)' ^''^ vertex of Umklapp type, 



ni . . . A, 

"i'fn'J.'.xUA^'^' °^ backscattering type, t^^^^^^'l'^^ Jq) 

and ('dj vertex respecting the SU(2) valley-pseudospin symmetry 
5, -5, 5', -5' c \ 

Aini...A4n4VM/* 



where the interaction vertex is defined as 
27re2 



V 



Aini...A4n4 



(q) 



£|q| 



-^Aim.Aanal q^l-^; 



A2?l2,A4n4 



(q)- 

(144) 



1. SU(2) valley symmetry 

One notices that, in contrast to the SU(2) symmetry 
associated with the physical spin, the Hamiltonian (jl43|) 
does not respect a similar valley-pseudospin symmetry 
due to possible inter-valley couplings. An SU(2) valley- 
pseudospin symmetry would be respected for the case 

= and ^2 = ^4; i-G- if the interaction vertex (jl44|) 



V 



«1...?4 

Ai?ii . . . A4n4 



(q) « 3^1. .^3^1 



(145) 



One may show, however, that the SU(2) valley- 
pseudospin symmetry is approximately respected when 
considering the different classes of interaction vertices de- 
picted in Fig. [13l 

• Consider the diagram in Fig. llSf a). which rep- 



resents a vertex of the type 

Aini...A4n 



Ai Til . . . A4n4 



(q) or 



V 



^(q). In this case, the particle on the 
left remains in the same valley whereas that on the 



right changes the valley. Such a process would re- 
quire a momentum transfer of ±K, i.e. of the wave 
vector connecting the two valleys, and therefore 
does not respect momentum conservation, in the 
absence of a magnetic field. Naturally, momentum 
is not a good quantum number here due to to the 
magnetic field, but momentum conservation man- 
ifests itself by an exponential suppression of such 
processes. In order to appreciate this point, we 
need to consider the Gaussian in the form factors 
(fMB and (fT42D . 



-^An.A'n' >>q/ ^ ^ ' 



(146) 



as discussed in Appendix \K\ [see Eq. (jA2p ]. One 
therefore sees that the interaction vertex contains 
a Gaussian term 



V 



(q) « 



^ ^-(|qr + |K|V4)/l,/2 



(147) 



where # represents an unimportant numerical fac- 
tor and where we have shifted the momentum 
q' = q ± K/2 in the second step. The processes 
associated with the diagram in Fig. fTST a) are thus 
exponentially suppressed in ~ 10"'/i?[T] and 

may safely be neglected in the range of physically 
accessible magnetic fields. 

• The same fate is reserved for the diagram in Fig. 
I13r b). which represents a process of Umklapp type. 
In this case, the vertex reads 



|q+K|2 + |q-K|^);l,/4 



(148) 



which is again exponentially small in Ig/a . 

• The situation is different for backscattering-type di- 
agrams [Fig. [T3Fc)]. in which case the interaction 
vertex is 



-(|q±K|^-Hq±K|^);|,/4 



(149) 



One may then redefine the wave vector q' = q±K, 
which is eventually an integration variable in the 
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interaction Hamiltonian (|143p . and the interaction 
vertex becomes 



, 27re^ , 



","11/2 (150) 



27re '2,2 



elKI 



9 



As an order of magnitude, with |K| 1/a, one then 
notices that the backscattering interaction vertex is 
suppressed by a factor oi a/ls ^ 0.005 x B\T] as 
compared to the leading energy scale /eIb- 

• The leading interaction vertex is therefore the 
SU(2) valley-pseudospin symmetric one depicted in 
Fig- ITST d). for which the rapidly oscillating contri- 
bution at K vanishes, as may be seen directly from 
the form factors (|142p . 



The above argument, which generalizes symmetry con- 



(lAlicea and Fishei 


, 2ooe 


; Doretto and Morals Smithl. 


7l: GoerbiE et all 
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Herbutl. 2007bl). shows that 



although the valley SU(2) symmetry is not an exact sym- 
metry, such as the SU(2) symmetry associated with the 
physical spin, it is approximately respected by the long- 
range Coulomb interaction. Valley-symmetry breaking 
terms are due to lattice effects beyond the continuum 
limit and therefore suppressed by the small factor a/Z^, 
which quantifies precisely corrections due to effects on the 
lattice scale. If one takes into account the additional spin 
degree of freedom, the resulting four-fold spin- valley de- 
generacy may then be described within the larger SU(4) 
symmetry group, which turns out to be relevant in the 
description of strong-correlation effects in partially filled 
LLs (Sec. El. 



2. SU(4) spin-valley symmetric Hamiltonian 

The SU(4)-symnictric part of the interaction Hamiltonian (|143|) finally reads 

HtnT ^\Y^ vZn^...XinS°^)P^^^iM^i'^-^)P^2n^MnM). (151) 

q Aini...A4n4 

where the symmetric interaction vertex is 
in terms of the reduced density operators 

PAn,A'n'(q) = Pxix-n'i^) = E E E (™ l^"''"'''! 4„,rn;« ..^A'™' ,m' , (153) 



where we have explicitly taken into account the spin in- 
dex a =t, i in the last line. 

We finally notice that the graphene form factors (|14ip 
may also be rewritten in terms of the LL form factors 

Fn,A^) = {n\e~'^-'^\n') , (154) 

which arise in a similar decomposition of the Coulomb 
interaction in Landau states in the non-relativistic 2D 
electron gas, as 

J^A«,A'«'(q) = i;in'^n-i,n'-i(q) + AA'2;2,*,F„,„-(q). 

(155) 

To summarize the differences and the similarities between 
the interaction Hamiltonians in graphene and the non- 
relativistic 2D electron system, one first realizes that its 
structure is the same if one replaces the LL form factor 
(|154p by the graphene form factors (|14ip and if one takes 
into account the larger (approximate) internal symmetry 
SU(4); due to the spin- valley degeneracy, instead of the 
spin SU(2) symmetry. 



In the remainder of this section, we neglect the 
symmetry-breaking part of the Hamiltonian and consider 
the Coulomb interaction to respect the SU(2) valley sym- 
metry. 

B. Particle-Hole Excitation Spectrum 

The considerations of the previous subsection allow us 
to discuss the role of the Coulomb interaction within 
a perturbative approach in the IQHE regime for v = 
±2(2n -I- 1), where the (non-interacting) ground state 
is non-degenerate and separated by the cyclotron gap 
y/2(hvF /lB)Wn + \ — \/n) from its excited states. Quite 
generally, the inter-LL transitions evolve into coher- 
ent collective excitations, as a consequence of these 
Coulomb interactions. Prominent examples in the 
non-relativistic 2D electron gas are the upper-hybrid 
mode (sometimes also called magneto-plasmon) , which 
is the magnetic-field counterpart of the usual 2D plas- 
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Figure 14 Zero-field particle-hole exctiation spectrum for doped 
graphene. (a) Possible intraband (I) and interband (II) single-pair 
excitations in doped graphene. The exctiations close to the Fermi 
energy may have a wave-vector transfer comprised between q = 
(la) and q = 2qp (lb), in terms of the Fermi wave vector qp. 
(b) Spectral function Imno(q,i.i;) in the wave- vector /energy plane. 
The regions corresponding to intra- and interband excitations are 
denoted by (I) and (II), respectively. 



mon (iGiuliani and Vignale . l2005l) . and magneto- excitons 
( Kallin and Halperinl . 1984 ). In the present subsection, 
we discuss how these modes manifest themselves in 
graphene in comparison with the non-relativistic 2D elec- 
tron gas. 



1. Graphene particle-hole excitation spectrum at B = 

Before discussing the particle-hole excitation spectrum 
(PHES) for graphene in the IQHE regime, we briefly re- 
view the on e for i? = as well as its associated collec- 
tive modes (Andol. 2006at Hwang and Das Sarmal 120071 : 



IShund . 119861: 1 Wunsch et al\ . \2QQ^ . Quite generallv. the 



PHES is determined by the spectral function 

1. 



5(q,c.) 



-Imn(q,cj), 



(156) 



which may be viewed as the spectral weight of the al- 
lowed particle-hole excitations, in terms of the polariz- 
ability H(q, cj), whi ch plays the role of a densit y -density 
respo nse function ( Giuliani and Vignalel 120051 : iMahanl . 
Il993l) . 

The particle-hole excitations for non-interacting elec- 
trons in doped graphene are depicted in Fig. 1141 ^^ In 
contrast to the PHES of electrons in a single parabolic 
band (the non-relativistic 2D electron gas) , there are two 
different types of excitations: intraband excitations [la- 
beled by I in Fig. fT^a)]. where both the electron and the 
hole reside in the conduction band (CB), and interband 
excitations [labeled by II in Fig. [Mi a)] , where an electron 
is promoted from the valence band (VB) to the CB. In un- 
doped graphene, there exist naturally only interband ex- 
citations (II). If the electron and the hole have an energy 



We consider here only the case of a Fermi energy tp in the con- 
duction band, for simplicity. 



close to the Fermi energy, the allowed excitations imply 
a wave- vector transfer that lies in between q = Q (la) and 
q = 2qp (lb). At non-zero values e of the transfered en- 
ergy, one needs to search for available quantum states at 
larger wave vectors, and the particle-hole pair wave vec- 
tor is then restricted to e/hvp < q < 2qF + e/hvp, as a 
consequence of the linear dispersion relation in graphene. 
This gives rise to the region I, which describes the intra- 
band particle-hole continuum, and its linear boundaries 
in the PHES described by the spectral function in Fig. 
Mh)- 

In addition to intraband excitations, one notices that 
interband excitations become possible above a threshold 
energy oi ep, where an electron at the top of the VB (at 
q = 0) may be promoted to an empty state slightly above 
the Fermi energy. The associated wave- vector transfer is 
naturally q = qp- The point {qp,ep) marks the bottom 
of the region II in Fig. [Ml b). which determines the region 
of allowed interband excitations (interband particle-hole 
continuum). Direct interband excitations with zero wave- 
vector transfer are possible above an energy of 2ep. 

Another aspect of the PHES in Fig. [M] is the strong 
concentration of spectral weight around the central diag- 
onal w = ?ii;F|q|- This concentration is a particularity 
of gr aphene due to the electrons' chirality (|Polini et al\ . 
I2008D . Indeed, if one considers a 2qp backscattering 
process in the vicinity of the Fermi energy in the CB, 
Eq. (jSO]) indicates that the chirality, i.e. the projec- 
tion of the sublatticc-pseudospin on the direction of the 
wave vector, is preserved. The inversion of the direc- 
tion of propagation in the 2qp process would therefore 
require an inversion of the A and B sublattices that is 
not supported by most of the scattering or interaction 
processes. This effect is reflected by a strong suppres- 
sion of the spectral weight when approaching the right 
boundary of the region I in the PHES associated with 
processes of the type lb in Fig. [T4l(a). Similarly, the 
conservation of the electrons' chirality ([5D|) favors 2qp 
processes in the interband region (II) and the suppres- 
sion of direct q = interband excitations of the type Ha 
in Fig. [MT a). Notice that, although the direction of the 
wave vector is inverted in a 2qF process, this indicates 
still the absence of backscattering because the group ve- 
locity V — Vqe^/S = Aii^q/jql remains unchanged - the 
change in the sign due to the inversion of the wave vector 
is indeed canceled by the one associated with the change 
of the band index. 



a. Formal calculation of the spectral function. In order 
to obtain the spectral function, it is apparent from 
Eq. (|156|) that one needs to calculate the polarizabil- 
ity H(q, w) of the 2D system, which may be found with 
the help of the Green's functions G(q, w), 

H(q,..) = -zTr J ^ ^ G(q', a;')G(q + q', c + c.'), 

(157) 
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Figure 15 (a) Particle- hole bubble diagram (polarizability), in 
terms of Green's functions G(q,i.i;) (lines). (b) Spectral func- 
tion Im ni{pyi(q, for doped graphene in the wave- vector /energy 
plane. The electron-electron interactions are taken into account 
within the RPA. We have chosen oq = 1 here. 



where Tr means the trace since the Greens functions 
are 2x2 matrices as a consequence of the matrix char- 
acter of the kinetic Hamiltonian. Diagrammatically, 
the polarizabihty may be represented by the so-cahed 
bubble diagram shown in Fig. [TST a). and one finds 
for non- interacting electrons i n graphene (Ando', '2006a!; 
Hwang and Das Sarmal. 120071: iPolini et all \2008i iShunal . 



19861 : IWunsch al\ . \20od) 



n°(q,-) = -| 



(4: 



K ] 



CAA'(q',q'+q) 

^ , — f. . H- nil) H- 7.n 

(158) 

where eq = A?iu_F|q| ~ ep is the energy of the quantum 
state V'A(q) measured from the Fermi energy e_F, 5 = 4 
takes into account the four-fold spin- valley degeneracy, 
and n{e^) is the Fermi-Dirac distribution function that 
reduces to a Heavyside step function n{e^) = 6(— Cq) 
at zero temperature. Equa tion (11581) is nothing other 
than the Lind hard function ( Giuliani and Vignalel 120051 : 
iMahanl . Il993| ). apart from the factor 



<^AA'(q',q' + q) = 



_ 1 + AA' COS0q',q'+q 



(159) 



in terms of the angle 6'q'.q'+q between q' and q' + q, 
which takes into account the particular chirality prop- 
erties of graphene ~ as already mentioned above, this 
chirality factor vanishes for backscattcring processes, i.e. 
for intraband (A = A') 2qp processes with 6'q',q/-|_q = tt 
as well as for interband (A = —A') q ~ processes with 
V.q'+q = or 2tt. 

Notice finally that the quantity 6 in Eq. p58p is an 
infinitesimal energy in the case of pure graphene and may 
be used (for finite values) as a phenomenological measure 
of the impurity broadening 6 ~ h/r, in terms of a life time 
T of the excitations. 



level of the random-phase approximation (RPA), which 
amounts to calculating a geometric series of bubble 
diagrams and which has shown to yield reliable re- 
sults for do ped graphen e (iHwang and Das Sarmal . 120071 : 



ISabio et all [2008. : .Wunsch et all |2006|). The RPA has 



also b een ap plied to undoped graphene ( Gonzalez et all 

f1, 'l99 9|), but its validity has been qu estioned 
mgad haraiah et all . [2008: .Kotov et al. [. I2OO7D because 
of the vanishing density of states, which wou ld require to 
take into account diagrams beyond the RPA ( Katsnelsonl . 
|2006[ ). The RPA polarizability then becomes 



q,wj 



(160) 



in terms of the polarizability (jl58|) for non-interacting 



electrons and the dielectric function 

2^e2 



-RPA 



(q,^) = 1 



^|q| 



n°(q,c.) 



(161) 



The spectral function associated with the RPA polar- 
izability (|160p . which is shown in Fig. [THb). reveals 
the characteristic coherent 2D plasmon mode, which 
corresponds to the solution of the implicit equation 
e^^^{q,ujpi) = and the dispersion relation of which 
reads 



2e'eF 



(162) 



in the small-g limit (jShungj . Il986l: IWunsch et ^ ■ |2006[) . 

Interestingly, this equation is valid also for non- 
relativ istic e l ectron s in conventional 2D electron sys- 
tems (ISternl . 119671) if one takes into account the dif- 
ference in the density dependence of the Fermi energy 
{tp = Tmei/mb for non-relativistic 2D electrons and 
ep = fivpy/nnei in graphene) as well as that in the Fermi 
velocity {vp = \f2tpjrn\t for non-relativistic electrons, 
as compared to a constant vp in graphene). Notice that 
the dispersion relation is restricted to small values of q 
(as compared to the Fermi wave vector lip\ whereas the 
numerical solution presented in Fig. [15] indicates that 
the asymptotic dependence of the plasmon mod e is in- 
deed given by the central d iagonal a;„/i(q) > vpq ( Shund . 
Il986l : IWunsch et all l2006( ). Therefore, contrary to the 
plasmon in 2D metals with a parabolic dispersion rela- 
tion, the plasmon in graphene does not enter region I, 
but only the interband particle-hole continuum (region 
II). In this region, the Landau damping is less efficient, 
and the coherence of the mode thus survives to a certain 
extent without decaying into incoherent particle-hole ex- 
citations. 



2. Polarizability for B / 

b. Polarizability in the random-phase approximation. The 

diagrammatic approach is particularly adapted for tak- In the case of a strong magnetic field applied pcrpen- 
ing into account the electronic interactions on the dicular to the graphene sheet, one needs to take into ac- 
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count the quantization of the kinetic energy into rela- larizabihty. One finds a similar expression for the zero- 
tivistic LLs described in Sec. III.A.H as well as the spine- temperature polarizability of non-interacting electrons as 
rial eigenfunctions ,„ in the calculation of the po- in Eq. (|158p . 



Q{eF - Xhuj'^/n) 
nB(q,w)=.9 2^ 

\n,X'n' 



I 



\ri.\ 



n'(q)l' 



(163) 



in terms of the graphene form factors (|14ip and the char- 
acteristic frequency = splvpjlB introduced in Sec. 
III.A.ll One notices that th e first part is nothing other 
than a Lindhard function (|Giuliani and Vignald . [200l 
iMahanl . [l993h for relativistic LLs filled up to the Fermi 
energy tp = ^(wf/^b)-\/2]Vf, which is chosen to be situ- 
ated between a completely filled (Np) and a completely 
empty + LL in the CB (IQHE regime). The second 
factor is the modulus square of the graphene form factors 
which plays the role of the chirality f actor C\,\i{al ^ q'-|-q) 
in the absence of a magnetic field (jRoldan et aly l2009l 

[2oi3;[shiin3l2oo3).^° 

As for the zero-field case, one may distinguish two con- 
tributions to the polarizability, one that may be viewed as 
a vacuum polarizability n'"'^(q, w) and that stems from 
interband excitations when the Fermi level is at the Dirac 
point, and a second one that comes from intraband exci- 
tations in the case of doped graphene. Because undoped 
graphene with zero carrier density does not correspond 
to an IQHE situation - as we have already discussed in 
Sec. IILAI the zero-energy LL n = is only half-filled 
then -, we define, here, the vacuum polarizability with 
respect to the completely filled zero-energy level. 

In order to describe more explicitly the different con- 
tributions to the polarizability , we define the auxiliary 
quantities ( Roldan et al\ . [20091 ) 



n 



Xn.X'n 



Xn.Xri 



(q)r 



\hvF\ 



-UJ-) (164) 



where w+ — > < 
—hjj — id and 



indicates the replacement hu + i5 



nA„(q,a;) = ^ ^ nA„,A'«'(q,a;) (165) 

A' n'=0 



A similar expression for t he po l arizability has also been obtained 
in Refs. l|Berman et all \200A iTahir and Sabeehl . l200d ) though 
with approximate form factors. 
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Figure 16 Particlc-holc excitation spectrum for graphene in a per- 
pendicular magnetic field. We have chosen Np = 3 in the CB and 
a LL broadening of 5 = O.Obhvp /Is (a,) and <5 = 0.2hvp /Ib (b). 
The ultaviolet cutoff is chosen such that Nc = 70. 



which verify nAn(q, i^) = -n_Ari(q,a;) 
larization may then be defined as 



The vacuum po- 



n™^(q,c.) 



(166) 



where Nc is a cutoff that delimits the validity of the 
continuum approximation. Notice that, already in the 
absence of a magnetic field, the validity of the contin- 
uum approximation is delimited by a maximal energy 
^ t. One may then introduce an upper level cutoff with 
the help of cn^ = fi{vp /Ib)\/'2Nc ^ t, that leads to 
Nc ^ 10^/i3[T], which is a rather high value even for 
strong magnetic fields. However, due to the fact that the 
separation between LLs in graphene decreases with n, 
it is always possible to obtain reliable semi-quantitative 
results from smaller values of Nc- 

The spectral function ^(q,^;) = — ImH^(q, a;)/7r is 
shown in Fig. [16] for Np = 3 and for two different values 
of the phenomenologically introduced LL broadening S. 
One notices first that the spectral weight is restricted to 
the two regions I and II corresponding to the intraband 
and interband particle-hole continuum, respectively, in 
the zero-field limit. This is not astonishing because the 
electron-hole pair wave vector remains a good quantum 
number also in the presence of a magnetic field and be- 
cause the overlap between the electron and hole wave 
functions is largest in these regions; if one considers the 
pair with its overall charge neutrality, its motion is unaf- 
fected by the magnetic field. Indeed, the pair momentum 
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may be viewed as the sum of the pseudomomenta asso- 
ciated with the guiding-center variable for the electron, 
R X ez/l%, and the hole — R' x e^/lg, respectively. Each 
of the pseudomomenta is naturally a constant of motion 
because so is the guiding center, as we have discussed in 
Sec. III. A. II One therefore obtains the relation 



q=ARxe2/?| or Ai? = |q|; 



2 



(167) 



where AR = R — R' is the distance between the guiding 
center of the electron and that of the hole. The bound- 
aries of the PHES in Fig. [12] may then be obtained from 
the decomposition ((90)) . which yields ?/ — ?7 < Ai? < 
■q' + T], with the help of the average values r] = {\t]\) = 
IbV^u + 1 and rj' = IsV'^^n' + 1, 



V2n' + 1- V2n +l<qlB< V2n' + 1 + V2n + 1. (168) 

Because the energy scales also with ^/n, one obtains the 
linear boundaries of the particle-hole continua as in the 
zero-field case mentioned above. 

In contrast to these similarities with the zero-field 
PHES, one notices a structure in the spectral weight 
that is due to the strong magnetic field. As a con- 
sequence of the relativistic LL quantization, the spec- 
tral weight corresponds to inter-LL transitions at ener- 
gies oj ~ ^/2h{vF /lB){\/n — \\fnl), where n > Np and 
n' < Np (for A = +) or n' > (for A = -). For larger 
values of Np, or quite generally when increasing the en- 
ergy, the level density increases due to the ^/n scaling of 
the LLs and the transitions. The LL structure is there- 
fore only visible in the lower part of PHES, in the clean 
limit S = O.O^Hvf/Ib [Fig. [TCTa)]. whereas the inter-LL 
transitions are blurred at larger energies or even for the 
lower transitions in the case of less clean samples [Fig. 
M.h), foT 5 = 0.2hvF/lB].^^ 

In addition to the (blurred) LL structure in the PHES, 
one notices another structure of the spectral weight, 
which is organised in lines parallel to the central diago- 
nal uj = hvF\<i\- This weight is again decreased when ap- 
proaching the right boundary of the intraband continuum 
(region I) and the left one of the intcrband continuum 
(region II) , due to the above-mentioned chirality proper- 
ties of electrons in graphene. The emergence of diagonal 
lines is a consequence of the graphene form factors (|14ip 
the modulus square of which intervenes in the polariza- 
tion function. Indeed, these form factors ^A(n+m),A'n(q) 
are (associated) Laguerre po lynomials with n + 1 zeros 
( Gradshtevn and Rvzhikl . [2000i ) due to the overlap be- 
tween the wave function of the hole in the level X'n and 
that of the electron in the LL A(n + m) (|Roldan et all 
These zeros in the inter-LL transitions are organ- 
ised in lines that disperse parallel to the central diagonal 
and thus give rise to zones of vanishing spectral weight. 
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Figure 17 Particle- hole exeitation spectrum for graphene in a per- 
pendicular magnetic field. The Coulomb interaction is taken into 
account within the RPA. We have chosen Np = 3 in the CB and a 
LL broadening of 5 = O.Obhvp/lg (a) and <5 = 0.2hvF /Ig (b). The 
ultaviolett cutoff is chosen such that Nc = 70. 



Interestingly, it is this structure of diagonal lines that 
survives in more disordered samples in which the hori- 
zontal lines associated with inter-LL transitions start to 
overlap, i.e. once the LL spacing is smaller than the level 
broadening (S.^^ 



3. Electron-electron interactions in the random-phase 
approximation: upper-hybrid mode and linear magnetoplasmons 

The PHES of non-interacting electrons in graphene 
gives already insight into the collective modes which 
one may expect once electron-electron interactions are 
taken into account. Indeed the regions of large spec- 
tral weight evolve into coherent collective excitations as 
a consequence of these interactions. Because the re- 
gions of large spectral weight are organised in lines par- 
allel to the central diagonal uj — ?iVF|q|, as mentioned 
above, one may expect that the dominant collective ex- 
citations are roughly linearly dispersing modes instead 
of the more conventional weakly dispersing magneto- 
excitons, which emerge fr om the inter-LL transitions 
(iKallin and Halperinl Il984l) . It has though been argued 
that such magneto-excitons may play a role at low en- 
ergie s in clean samples with low doping ( Ivengar et all . 
I2OO7I ) and that they may renormalize the cyclotron en - 
ergy at zero wave vector (jBvchkov and Martinet . |2008[) . 

As in the zero-field case, we take into account the 
Coulomb interaction within the RPA [see Eq. (|160[) ]. 
The resulting spectral function is shown in Fig. [T7]for the 
same choice of parameters as in the non-interacting case 
(Fig. [T5|) . Furthermore, we have chosen a dimensionless 
interaction parameter ac = 1, here, which corresponds 



The value 5 = 0.2hvp /Ib is a reasonable estim ate for today 's 
exfoliated graphene samples on an Si02 substrate l lAndol . r2007ar) . 



This behavior is in stark contrast to that of non-relativistic 2D 
electrons, where the LL spacing is constant and given by the 
cyclotron energy heB/rrif,. If this quantity is larger than the 
level broadening 5, there is no qualitative difference between low 
and high energies, and the horizontal lines associated with the 
inter-LL excitations (multiples of the cyclotron energy) remain 
well separated. 
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to a dielectric constant of e ~ 2. 

One notices the prominent mode that evolves in the 
originally forbidden region for particle-hole excitations. 
This mode, which is called upper-hybrid mode, is the 
magnetic-field descendent of the 2D plasmon mode (|162p 
and acquires a density-dependent cyclotron gap ujc ~ 
eBvp/ep = eBvp /hvp^nTLei. Its dispersion relation in 
the small-g limit is then given by 



(169) 



' 2he'^vpy^7Tnei 



eBvj, 
hvpy/nnei 



as may be shown easily within a hydrodynamic ap- 
proach that has been successfully applied to the upper- 
hybrid mode in non-r elativistic 2D electron systems 
( Chin and Quinnl . Il974l) . It is apparent from Fig. [17] 
that this mode remains coherent also within the region 
II, which corresponds to the B = interband particle- 
hole continuum. 

In addition to the upper-hybrid mode, one notices lin- 
early dispersing coherent modes in the regions I and II 
that emerge from the lines of large spectral weight in 
the non-interacting PHES, as expected from the quali- 
tative discussion above . We may ca l l these mod es lin- 
ear magneto-plasmons ( Roldan et aZ.I . [2009ll2010l ) in or- 
der to distinguish them clearly from the upper-hybrid 
mode (|169p and the weakly dispersin g magneto-excitons 
at low doping ( Ivengar et all |2007| ). These modes are 
more prominent in the interband than in the intraband 
region although they are also visible there. They arc gen- 
uine to graphene with its characteristic V Bn LLs that 
inevitably overlap in energy, above a critical LL n, if 
level broadening is taken into account, and they may not 
be captured in the usual magneto-exciton approximation 
where the collective mode s are adiabatically connected 
to the inter-LL excitations tevchkov and Ma rtinet 120081 : 
llvengar et oil . [2007; .Kallin and Halpcrin. . .1983 ). 



4. Dielectric function and static screening 

We terminate this section with a brief discussion of the 
dielectric function (|161|1 in the static limit 



-IiPA/„\ _ ^RPA 



(q)=e^*^^(q,cj = 0) = l- 



2TTe' 



n"(q,c^ = 0), 
(170) 

by comparing the B 7^ to the zero-field case, as shown 
in Fig. [THl As mentioned above, one may distinguish 
two separate contributions to the static polarizability, 
the vacuum polarizability n^"'^(q) = n"°^(q, w = 0) due 
to interband excitations and the intraband contribution 
n''°^'(q) = n''°P(q, uj — 0), which is only present in doped 
graphene. 



n"(q,a; = 0) = n™"(q) + U'^°P{q). 



(171) 
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Figure 18 Static polarization function n''(q, t»; = 0) for non- 
interacting electrons in graphene without (a) and with (b) a, mag- 
netic field. To compare both cases, we have chosen a Fermi wave 
vector qp = \/2Np + 1/Ib = V7/Ib that corresponds to Np = 3. 
The full black line represents the total polarizability, whereas the 
red dotted and the blue dashed lines show the intraband and the 
interband contributions, respectively. From (.Roldan et ai . 2010.) . 



One notices that up to 2qp, the zero-field static po- 
larizability [Fig. [TST a)] remains constant. Indeed, the 
interband contribution (bl ue dashed line) increases lin- 
early with the wave vector (|Andcil . [2006al : lGonzalez et all 
Il999ll 



n^'''=(q) = 



(172) 



and thus cancels the linear decrease of the intraband con- 
tribution (red dotted line). 



n'^°P(|q|<2g^):^p(6^^)fl-^ 



where 



(173) 



(174) 



is the density of states per unit area at the Fermi energy. 
At wave vectors larger than 2qp, the intraband contri- 
bution dies out, and the polarizability is dominated by 
interband excitations. 

With the help of these two contributions, we may 
rewrite the static dielectric function p70[) as 



.RPA 



(q) 



1 



n'iop(q) 



where we have defined 

^ — ^RPA 



(|q| C5o) 



(175) 



(176) 



i.e. the value that the static dielectric function ap- 
proaches at large wave vectors. Notice that this is pre- 
cisely the RPA result for the intrinsic dielect ric constant 
for undoped graphene ( Gonzalez et aLl . ll999( ). The above 
form of the static dielectric function may be cast into a 
Thomas-Fermi form. 



e 



TP 



(q) 



So 



1 



1L 
q 



in terms of the effective coupling constant 
£00 H-7raG/2' 



(177) 



(178) 
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Figure 19 Effective coupling constant Oq as a function of tlie bare 
coupling aQ. The dashed line indicates the asymptotic value 2/tt 
obtained for large values of the bare coupling (r^ 3> !)• 



which is plotted in Fig. [T9l 

One notices that interband excitations yield a contri- 
bution to the dielectric constant that originally takes into 
account the dielectric environment in which the graphene 
sheet is embedded, £—)•£*= eSao- This is a direct conse- 
quence of the linear behavior of the vacuum polarization 
(|172p as a function of the wave vector and thus specific 
to graphene. Furthermore, one may also define an effec- 
tive Thomas-Fermi wave vector q^p = qxp/sao = ctQlF^ 
which describes the screening length in the presence of 
the vacuum polarization. As a consequence of the sat- 
uration of the effective coupling constant p78p at large 
values of ac, the effective Thomas-Fermi vector is thus 
always on the order of the Fermi wave vector unless 
ac ^ 1, where Uq ~ aa- The relevant effective param- 
eters are summarized in the table below for freestanding 
graphene and graphene on mainly used substrates. 



graphene 


e 


aa 


£ oo 




in vacuum 


1 


2.2 


4.5 


0.5 


on Si02 


2.5 


0.9 


2.4 


0.38 


on h-BN 


2.3 


1 


2.4 


0.39 


on SiC 


5.5 


0.4 


1.6 


0.25 



Table I Dielectric constant e, Sao, bare coupling aa, and effective 
coupling a*Q for graphene in vacuum and popular substrates. 



Finally, the screened Coulomb interaction potential 
may be approximated as 



Vscr{q) 



eeoo(l + a*(jqF/q)q' 



(179) 



One notices from this expression that processes at wave 
vectors q ^ qp, where the interband excitations play a 



The dielectric constant e is then the average of the dielectric 
constant of the substrate material and that of the vacuum. 




Figure 20 Static dielectric function for graphene in the IQHE 
regime for Np = 1, 2 and 3 (blue, red and green curves, 
ing order). From l|Roldan et al l \201(t ). 



m mcreas- 



negligiblc role [see Fig. [TST a)]. are still governed by the 
bare coupling constant ao ^ Vgcril ^ <lF)l^vpq. How- 
ever, processes at or above the Fermi wave vector, such 
as those that are relevant in the electronic transport, are 
characterized by the effective coupling constant a*Q ^ 
Vscril ^ If) l^vpq, which saturates at a value of 2/7r as 
mentioned above. If we consider the value (jlSSp . ac — 
2.2, for the bare coupling constant of graphene in vac- 
uum, the effective coupling is roughly four times smaller 
a*Q ~ 0.5, such that the electrons in doped graphene ap- 
proach the weak-coupling limit. The situation is differ- 
ent in undoped graphene, wh er e recent renormaliza tion- 
group (|Herbut et al\ . l2009ai rbl: lJuricic et aR. l2009l) and 



lattice g auge theoretical calculations ( Prut and Lahdel . 



lattice ga 

l2009allbf ) indicate a flow towards strong coupling at mod- 
erate values of ac- 

In Fig. fTSl b). we have plotted the static polarizability 
for graphene in the IQHE regime. Qualitatively, the re- 
sult agrees with the zero-field behavior, with a (roughly) 
linearly increasing vacuum polarizability and a decreas- 
ing intraband contribution, apart from some superim- 
posed oscillations due to the overlap functions that are 
refiected by the form factors (|14ip . An important differ- 
ence is manifest in the small wave- vector limit of the po- 
larizability. In contrast to the zero-field case, where the 
polarizability saturates at a non-zero density of states, 
the system is gapped in the IQHE regime with a result- 
ing vanishing density of states at the Fermi energy. This 
gives rise to a q^ behavior of the polarizability at small 
wave vectors. Furthermore, the static die l ectric function , 
which is shown in Fig. [20] (|Roldan et all l2010l : IShizuval . 
|2007[ ). does no longer diverge in this limit, contrary to 
the zero-field Thomas- Fermi result (|177p . Indeed, the 
small-g behavior may be approximated as 



.RPA 



(9) 



1 oc acNp^^qlB, 



(180) 



which is the same as for non-relativistic 2D elec- 
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trons ( Aleiner and Glazmanl . Il995l ).^^ The maximum 



of the static dielectric function is obtained at qIb ^ 
l/qpl B ^ 1/V 2Nf + 1, at the value £niax ^ e^^'^iq ^ 
1/IbV^^f + 1) ~ aoNp. It therefore scales as Smax oc 
Np in co ntrast to a y/Np scali n g in n on-relativistic 2D 
systems ( Aleiner and Glazmanl . 19951) . At large wave 



vectors, the static dielectric function saturates at the 
same value Eoo as in zero magnetic field. 



IV. MAGNETO-PHONON RESONANCE IN GRAPHENE 

In the previous section, we have discussed the role of 
electron-electron interactions in the IQHE regime, where 
a perturbative (diagrammatic) approach may be applied. 
Similarly, one may treat the electron-phonon interaction 
in a perturbative manner in this regime. This is the 
subject of the present section, before discussing again 
electron-electron interactions in the strong coupling limit 
of partially filled LLs (Sec. |V|. 

As a consequence of the honeycomb-lattice structure 
of graphcne, with two incquivalent sublattices, there are 
four in-plane phonons, two acoustic and two optical ones. 
Each phonon type occurs in a longitudinal (longitudinal 
acoustic, LA, and longitudinal optical, LO) and a trans- 
verse (transverse acoustic, TA, and transverse optical, 
TO) mode. In addition to these four phonons, one finds 
two out-of-plane phonons, one acoustic (ZA) and one op- 
tical (ZO) [for a re view of phonons in graphene, s ee Refs. 
(ISaito et all Il998l ) and (IWirtz and Rubiol . |2004[ )]. Here, 
we concentrate on the in-plane optical phonons LO and 
TO, which couple to the electronic degrees of freedom. 
More specifically, we discuss these phonons at the F point 
(E2g modes) in the center of the first BZ, which are at- 
tributed to the G-band at hiu^h — 0.2 eV in the Raman 



spectra [see e.g. Refs. (Ferrari et aL . 12006 : iGraf et~al. 



2007D ] 



20071: iGupta et all 120061: iPisana et all 120071: lYan et at 



One of the most prominent effects of electron-phonon 
coupling in met als and semi conductors is the so-called 
Kohn anomaly ( Kohnl . Il959l ). which consists of a sin- 
gularity in the phonon dispersion due to a singular- 
ity in the electronic density-density response function. 
The analog of the Kohn anomaly in graphene yields a 
logarithmic divergence of the phonon frequency when 
the bare phon on frequency coincides with twice the 
Fermi energv (lAndol. 2006b • ICastro Neto and Guineal . 



l2007l : iLazzeri and Mauril |2006|) . We investigate how this 
renormalization manifests its e lf in graphene i n a s trong 
mag netic fi eld (|Andol . l2007bl: iGoerbig et all . |2007[ ). In 
Sec. IIV.AI we consider the specific form of the electron- 
phonon coupling and discuss its consequences for the 



■^'^ Notice, however, that the expression becomes exact only in the 
large- A'^i^ limit and that in non-relativistic 2D electron systems, 
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phonon 

electron , 
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electron-phonon 
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Figure 21 Electron-phonon coupling in graphene. (a) Optical 
phonons in graphene, with a wave vector q in the vicinity of the F 
point at the center of the first BZ [see part (c)]. The amplitude of 
the LO phonon is in the direction of propagation (black arrows), 
that of the TO phonon perpendicular (red arrows). The optical 
phonons modify the bond lengths of the honeycomb lattice, (h) As 
a consequence of the modified bond lengths, the electronic hopping 
is varied, and the electron-phonon coupling is off-diagonal in the 
sublattice index. 



renormalization of the optical phonons at the F point 
in Sec. IIV.BI More specifically, we consider both 
non-resonant (Sec. lIV.B.ip and resonant coupling (Sec. 
IIV.B.2p . the latter being specific to graphene in a mag- 
netic field when the phonon is in resonance with an al- 
l owed inter-LL trans ition (magneto-phonon resonance) 
(IGoerbig et aLl . l2007D . 



A. Electron-Phonon Coupling 

The LO and TO phonons in graphene are schematically 
represented in Fig. l2lT a). As already mentioned above, 
we concentrate on phonons at small wave vectors, in the 
vicinity of the F point. The origin of the electron-phonon 
coupling may easily be understood from the variation of 
the bond length caused by the phonon, which affects the 
electronic hopping amplitude between nn carbon atoms. 
As we have discussed in Sec. II. D[ the efi' ect may be quan - 
tified with the help of Harrison's law ( Harrison! . Il98ll ). 
which yields dt/da ~ -4.3 eV/A [see Eq. The 
order of magnitude for the bare electron-phonon energy 
is then obtained by multiplying this variation with the 
typical length scale y/h/Muph, which characterizes the 
amplitude of a lattice vibration of frequency ujph within 
the harmonic approximation. The intervening mass M is 
that of the carbon a t om. Indeed, a tight - bindin g calcu- 
lation (jAndol . l2006bl : llshikawa and Andol . l2006l ) corrob- 
orates this argument, apart from a numerical prefactor 
3/2, and yields a bare electron-phonon coupling 



5= -oTT- 




0.26 eV. 



(181) 



This value ag r ees w ell with ab-initio calculations 
( Piscanec et although it is though slightly 

lower than the value determined experimentally , which 
ranges froni q ~ 0.3 eV ( Faugeras et all 20091 : 



aLl. 12007!) 



the coupling constant rs also depends on A^^, 



-1/2 



iPisana et all . I2007D to g ~ 0.36 eV (|Yan et 

Furthermore, one notices that, because the electron- 
phonon coupling is mediated by a bond variation between 
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sites that belong to two different sublattices, the coupling 
constant is off-diagonal in the sublattice basis. This is 
diagrammatically depicted in Fig. I^TTb). 



where cr A u(r) = [<t x u(r)]^ = a^Uy{r) — (T''u^(r) is 
the 2D cross product between the Pauli matrices and the 
displacement field. 



1. Coupling Hamiltonian 



2. Hamiltonian in terms of magneto-exciton operators 



The above considerations help us to understand the 
different terms in the Hamiltonian, 



H = H, 



ph 



H. 



coupl 



which serves as the basis for the analysis of the electron- 
phonon coupling. The Hamiltonian for 2D electrons in a 
magnetic field, 



off,? 
B 



V'?(r) 



(182) 



may be written, in second quantization, in terms of the 
one-particle Hamiltonian (|75p and the fermionic fields 



Xn.m 



V'An,m;?(r)cA 



where '0An,m;5(r) is the wave function in position space 
associated with the spinor (|138p . 

The lattice vibration is characterized by the relative 
displacement u(r) between the two sublattices, which 
may be decomposed in terms of the bosonic operators 
6^,q and 6t q, 



u(r) = 




2NucMLUi,{q) 



(183) 

where e^^q denotes the two possible linear polariza- 
tions {fi = LO,TO) at the wave vector q and Nuc = 
A/{3y/3a'^ /2) is the number of unit cells in the system. 
The phonon Hamiltonian then reads, on the level of the 
harmonic approximation. 



(184) 



in terms of the phonon dispersion Ld^{q). Notice that, at 
the r point, the frequencies of the LO and TO phonons 
coincide, and one has ujph = w^(q = 0). It is then con- 
venient to describe the phonon modes in terms of circu- 
larly polarized modes, UQ{r) = [^^.(r) -I- iUy{r)]/\/2 and 

Mo(r) = WoW- 

Finally, taking into account the above considerations 

on the elec tron-p honon co upling, the coupling Hamilto- 
nian reads (lAndo . 2006b: Castro Neto and Guineal . [2007t 
llshikawa and Andol . 1200611 



H, 



coupl 



2Mu 



ph 



J d'r^ir) [<TAu(r)]V'c(r), 

(185) 



As a consequence of the off-diagonal character of the 
electron-phonon coupling pSSp . one notices that the in- 
tervening matrix elements are proportional to (5„.„±i, and 
one thus obtains the selection rules 



An -> A'(n± 1), 



(186) 



in analogy with the magneto-optical selection rules dis- 
cussed in Sec. III. A. II Furthermore, if we fix the energy 
of the dipole transition (|105p to be^^ A„ = A„ ,\=- = 
\/2h{vF IIb)W''^ -1-14- v^)j there are two possible tran- 
sitions, which may be distinguished by the circular po- 
larization of the phonon they are coupled to. Indeed, 
the form of the electron-phonon coupling (|185p indicates 
that the O-polarized phonon is coupled to the transi- 
tion — {n + l) — > +n, whereas the 0-polarized phonon 
couples to the —n — » +{n + 1) interband transition 
(jGoerbig et aLl . l2007D . 

It is then convenient to introduce magneto-exciton op- 
erators, associated with the above-mentioned intcr-LL 
transitions 



<P0\^^^) - / . '^+Ti.m:f'^-("+l).™;^' 



^ ra 

where the index A =0, O characterizes the angular mo- 
mentum of the excitation and where the normalization 
factors 

= ^2(1 + 5,,^o)Nb [v-(n+l) - V+n] 

and TVp = ^2{l + 5n.o)NB [v-n - v+^n+i)] 

are used to ensure the bosonic commutation rela- 
tions of the exciton operators, [(f>A{n,S,),(f)\,{n' = 
5A,A'5^,^'^n,n' ^ including the two-fold spin degeneracy. 
These commutation relations are obtained within the 



mean-field approximation with (c^ 



An,m;^ 

(i5a,- -I- S\^+i^\n), where < i^An < 1 
is the partial filling factor of the n-th LL, normalized to 
1. 

One notices that the magneto-exciton operators are, 
apart from a normalization constant, nothing other 



We consider mainly interband transitions here, which may have a 
chance to be in resonance with the phonon of energy huiph ~ 0.2 
eV. 
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^l-(n+l)(*l 



than the reduced density operators p40p . (/(^(n,^) cx 
0) and (/>^(n,^) cx p^^„+i) _„(q = 
0), respectively, at zero wave vector. Notice fur- 
thermore that, because of the relative sign ^ between 
the kinetic part p82p and the clcctron-phonon cou- 
pling Hamiltonian (|185p . the optical phonons couple 
to the valley-anti-symmetric magneto-exciton combina- 
tion (l)A,as{n) = [0^(71, C = +) - <^^(n,C = -)]/V2- 
This nee ds to be contrasted to the magneto-optical 
coupling (lAbergel and Fal'kol . 120071: llvengar aZI . 120071: 
ISadowski et ad . l2006t) . where the photon couples to the 



valley-symmetric mode 0^,s(n) — [^yil^T-, C = +) + 
q>A{n,^ = -)]/V2. 

The magneto-exciton operators ([T87| allow one to 
rewrite the electron-phonon Hamiltonian at the T point 
(q = 0) in a bosonic form as ( Goerbig et aLl . [2007() 



T=s,as A,n 



Y.9A{n) [b\(i^A;as{n) +bA<i>\,M] > (188) 



A.n 



in terms of the effective coupling constants 



90{n) = g\J (1 + 5nfi)lyjv_(n + i) - J/+„, 

and go(n) = (1 + Sn.ohy/v^n ~ j^+(„+i), (189) 

with the constant 7 = 3-\/3a^/27r/^. One therefore re- 
marks that, although the bare coupling constant g is 
rather large [see Eq. (|18ip ]. the effective coupling is re- 
duced by a factor of a/ls, 



9A(n)^g—^l...2mcV./BfT]. (190) 



B. Phonon Renormalization and Raman Spectroscopy 



The Hamiltonian (|188p shows that a phonon may be 
destroyed by exciting a magneto-exciton, and the associ- 
ated Dyson equation for the dressed phonon propagator 
D(cj) reads 

DA{io) = Do{lo) + Do{uj)xa{<^)Da{<^), (191) 
in terms of the bare bosonic phonon propagator 

1 2w 



■'ph 



and 



XA 



(-)- E 



7i=Nf + 1 J^f 



(192) 



(193) 



The form of the last expression is transparent; the 
magneto-exciton is a boson, and its propagator is there- 
fore of the same form as that of the bare phonon. It 



is equivalent to a particle-hole propagation associated 
with a polarization bubble [see Fig. fTSl a)]. but the ex- 
pression ()193p also takes into account the square of the 
effective coupling constant which is due to the double 
occurence of the electron-phonon coupling - first when 
the phonon is converted into a magneto-exciton and the 
second time when the magneto-exciton is destroyed by 
creating a phonon. The last term in Eq. (|193p takes 
into account the (only) possible intra-band magneto- 
exciton from the last filled LL Np to Np -\- 1 with energy 
Ajv^ = V2{hvp/lB)iVW 



1 — y/Np), which we have 
omitted in the Hamiltonian (|188p because it is irrelevant 
for resonant coupling. The parameter Nc is the same 
high-energy cutoff, defined by ejv^ = h{vp /Ib)\/'2.Nc ~ t, 
as in Sec. lIH.BI of the preceding chapter. 

The renormalizcd phonon frequencies loa may be ob- 
tained from the Dyson equation p9ip . by searching the 
poles of the dressed phonon propagator 

DAiio)-^ = = i5^(cD^)-i - Xa{C:^a). (194) 
and one finds (|Andol . l2007bl : lGoerbig et aLl . 120071 ) 



ph 



E ■ 



Ang\{n) , KM,g\{Np) 



Al 



A2 



(195) 



1. Non-resonant coupling and Kohn anomaly 



Before discussing resonant coupling, i.e. when the 
phonon frequency is in resonance with a possible inter- 
LL excitation, in a strong magnetic field, we comment on 
the relation between Eq. p95p and the (non-resonant) 
renormalization of the phonon frequency in zero mag- 
netic field. The zero-field limit may indeed be obtained 
from Eq. (|195p if one replaces the sum by an in- 
tegral j dn, i.e. if the spacing between the LLs van- 
ishes, Ajvf — > 0. Linearizing Eq. p95p . if one replaces 
LJA ^ i-^ph in the denominators, and using the approxi- 
mation ^/n + vn+1 ~ 2-y/n yields 



y^^-^ln 
4 



LOph , f ujph + 2^/2Npvp/li 



UJph - 2y/2NpVp/lB 

(196) 

where A = {2/V^TT){g/t)^ ~ 3.3 x 10"^ is the dimen- 
sionles s electron-phonon co upling constant introduced in 
Refs. (|Andol . l2006bll2007bD . and 



Wo — LJph + 



ph 



A2 



(197) 



is the physical phonon frequency at zero doping. In- 
deed, the frequency ujph is not relevant in a physi- 
cal measurement in graphene even if it occurs in the 
Hamiltonian, but one measures ^-t zero doping and 
B = 0. Equat ion (I196D coincides precisely with the 
zero-field result ( Andol . 2006bl : ICastro Neto and Guineal . 
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boson modes read ijAndcl . l2007bl: iGoerbig et ^ . l2007t) 
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Figure 22 ('aj From Ref. l|Goerbig et a/.l . [2007h . Anticrossing of 
coupled phonon-magneto-exciton mo des as a func tion of the mag- 
netic field, (b) From Ref. (Faugeras et g/.l . 120091 '). Colour plot of 
the Raman spectra as a function of the magnetic field. The contin- 
uous white lines indicate the magnetic field for which the phonon is 
in resonance with an inter-LL excitation of energy A„. Top: data 
of the full S-field range. Bottom: zoom on the range from to 10 
T. 



l2007l : lLazzeri and Mauri |2006| ) if one identifies the chem- 
ical potential with the energy of th e last filled LL, /.i = 
V^jT^hvp/lB (|Goerbig et aJ . 120071 ). 



2. Resonant coupling 

Apart from the non-resonant coupling discussed in the 
preceding section, the high-field elcctron-phonon cou- 
pling reveals a linear effect when the phonon is in res- 
onance with a particular magneto-exciton, huph ~ A„. 
In this case, the sum on the right-hand side in Eq. 
(|195p is dominated by a single term and may be ap- 
proximated by 2{u!pii/h)g'^{n)/ {hujj^ — A„). This re- 
sults in a fine structure of mixed phonon-magneto- 
exciton modes, (ftAjasin) cos 6 + 6^sin0 with frequency 
a)J^ and (/)yi as(n) sin 0—6^ cos with frequency (D^ [where 
cot 29 — (A„— ?kZ)o)/25_4]. The frequencies of these mixed 



,~,± 



in) 



h 



- ^0 T 




(198) 

and the resulting phonon-magneto-exciton anticrossing 
is depicted in Fig. I^ a). 

The above-mentioned anticrossing of the coupled 
phonon-magneto-exciton modes has been observed in 
recent Raman experiments on epitaxial graphene. Re- 
member that Raman spectroscopy is sensitive to the 
phon on component of the mixed modes ( Faugeras et all . 
l2009ll . The results are shown in Fig. [SSf b) and corrob- 
orate the theoretical ly expected behavior ( Andol . [2007bt 
iGoerbig et all I2007D . Indeed, one may obtain the oscil- 
lating behavior from a numerical solution of Eq. (|195p if 
one expresses the equation in terms of wq instead of ujph 
and if one takes into account a finite broadening of the 
levels. If the phonon is out of resonance with an inter-LL 
transition, its frequency is essentially field-independent 
and coincides with the energy of the F,2g line at 1586 
cm~^ ~ 0.2 eV. When it approaches the resonance (by 
increasing the magnetic field), its energy is shifted up- 
wards as a consequence of the anticrossing but rapidly 
dies out in intensity once the magneto-exciton component 
becomes dominant in the ui^ mode. Upon further in- 
crease of the magnetic field, the mode becomes more 
phonon-like and therefore visible in the Raman spectra. 

The fine structure of the high-field resonant electron- 
phonon coupling may furthermore be investigated by 
sweeping the chemical potential when the magnetic field 
is held fixed at resonance. The effect is most pronounced 
for the resonance hujph — A„=Oi which is expected at 
S ~ 30 T [see Fig. ESJa)]. In this case, the mode 
consists of an equal-weight superposition of the phonon 
and the magneto-exciton {cos 9 = sin 9 = 1/^/2), and 
the Fi2g band would appear as two lines, at the energies 
hui^ = hCjQ ± for the case of undoped graphene. 
With the above estimation (|18ip for the bare electron- 
phonon coupling constant, one obtains for the line split- 
ting 2g_4 ~ 16 meV 130 cm~^), w hich largely exceeds 



the G-band width observed in Refs. (Ferrari et al. 



Graf et al. 1 . 120071 : iGuDta" et al. 1. 120061 : 
Yan et ad [20071 ). 



Pisana et al. 



2006; 



2007 



It is apparent from the expressions (|190|) for the ef- 
fective coupling constants and g^ that the splitting 
may be controled by the LL filling factor. Exactly at 
zero doping, both coupling constants coincide, (?o = 50i 
but upon electron doping the transition — 1 — > associ- 
ated with the 0-polarization becomes weaker due to the 
reduced number of final states in n = 0, whereas the 
— > +1 transition (with polarization O) is strengthened. 



^® Notice, however, that only an oscillation of the phonon mode, 
and not a split ting, has been observ ed in the experiment by 



Faugeras et al. l lFaugeras et aZ.Ll2009l ). 
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As a consequence, the associated coupling constants are 
increased and decreased, respectively, until the coupling 
constant vanishes at u — 2. 

The above-mentioned filling-fac tor dependence has a 
direct impact on the Raman lines ( Goerbig et aZj . [2007t ). 
Whereas at = 0, one expects two lines separated by 
the energy 2gQ = 2gQ , the degeneracy in the circular po- 
larization is lifted between < v < 2?'' One therefore 
expects to observe four lines instead of two, where the in- 
ner ones are associated with the polarization O, whereas 
the outer ones with increased splitting correspond to the 
opposite polarization O. The separation between the in- 
ner lines vanishes then at ^ = 2, where the splitting of the 
outer lines is maximal and where one expects to observe 
three lines. 



in a partially filled relativistic LL (Sec. IV. A|) . This 
model yields a qualitative understanding of the above- 
mentioned correlated electronic phases in the context of 
graphene, as compared to non-relativistic electrons. In 
Sec. IV. Bl we apply this model to the quantum Hall fcrro- 
magnetism with an internal SU(4) symmetry that is the 
relevant symmetry in graphene LLs and discuss its rela- 
tion with the experimentall y observed degenera cy lifting 
of the zero-energy LL n = (jZhang et a/.l . l2006l ). We ter- 
minate this section with a review of the specific FQHE in 
graphene (Sec. IV.CI). which has recently been observed in 
the two-terminal ( Bolotin et a/.l . [20091 : IDu et ai . 20091) as 



well as in the four-te rminal geometry (jPean et all 20111: 
iGhahari et ^ . l201lh . 



V. ELECTRONIC CORRELATIONS IN PARTIALLY 
FILLED LANDAU LEVELS 

This last section is devoted to the physics of in- 
teracting electrons in the strong-correlation limit of a 
partially filled LL. The motivation stems from non- 
relativistic quantum Hall systems in GaAs heterostruc- 
tures, where these correlations lead to the formation 
of incompressible quantum-liquid phases, w hich display 
the f ractional quantum Hall effect (FQHE) (jTsui et all . 
Il982t ). as well as of exotic electron- solid phases , 
such as the high-fie l d Wig ner crystal (jAndrei et all 
Il988l : IWilliams et all Il99ll) or th e theoretical l y pre- 
dicted bubble a nd st ri pe phases (iFogler et all 11996 
iKoulakov et al ], 119961 : iMoessner and Chalkeil 
The latter 
aniso tropic 



are likely 
tra nspor t 



to be at the origin of highly 

. . propertie s at half-filled higher 

LLs (|Du et all Il999l : iLillv et all Il999[ ). pa rticular elec- 
tron t ransport un der microwave irradiation ( Lewis et all 
l2005l . |2004| . |2002|) . and an intrigui ng reentrance of the 
IQHE in n = 1 an d n = 2 (jCooper et all Il999l: 
lEisenstein et adl2002| ). 

It is therefore natural to ask whether such strongly- 
correlated phases exist also in graphene and if so what 
the differences are with respect to non-relativistic 2D 
electrons. Moreover, the fact that the electrons reside 
at the surface opens up the possibility to probe these 
phases by spectroscopic means, such as scanning tunnel- 
ing spectroscopy, which has already been applied success- 
fully in the analysis of the electr on density distribution 
of exfoliated ( Martin et all l2008l ) and epitaxial graphene 
(iMallet adl2007l). as well as graphene on graphite sub- 



strates (|Li et allmom . 



After a short discussion of the Coulomb interaction in 
graphene as compared to non-relativistic 2D electrons, 
we introduce the basic model of interacting electrons 



We present the argument for a Fermi energy in the CB, i.e. u > 0, 
but the situation is generic, and the argument also apphes in the 
VB if one interchanges the polarizations. 



A. Electrons in a Single Relativistic Landau Level 

Quite generally, the origin of strongly-correlated elec- 
tron phases is a quenched kinetic energy, where the par- 
tially filled LL is separated by the cyclotron gap from 
the neighboring ones such that inter-LL excitations con- 
stitute high-energy degrees of freedom. The Coulomb in- 
teraction, which may though be small with respect to the 
cyclotron gap, remains then as the only relevant energy 
scale which dominates the low-energy degrees of free- 
dom if we can neglect disorder effects. This leads to the 
seemingly counter-intuitive finding of strongly-correlated 
phases in weakly-correlated matter. 

In order to quantify the degree of separation between 
the energy scales, one may use a similar argument as 
the one that led us to the definition of the dimension- 
less interaction parameter (|13ip . introduced at the be- 
ginning of Sec. mil Indeed, one needs to compare the 
Coulomb interaction energy Eint = e^ / eRc at the char- 
acteristic length scale Rc = IsV^n + 1 to the LL spacing 
Hujc = heB/mb, where we concentrate on non-relativistic 
electrons first. 



with vp{n^B) = Rcloc- 

(199) 



hevpin, B) 



If one identifies the Fermi wave vector kp ~ \/2n/lB: one 
obtains the same expression as for the zero-field coupling 
constant (|132p . 



= Ts = 



1,-1 



1 



(200) 



This means that the degree of LL mixing is still governed 
by Tg, and the inter-LL excitations are well separated 
from the low-energy intra-LL degrees of freedom unless 
Ts becomes very large. Notice, however, that ^ 1 in 
most 2D electron systems. 

In the case of partially filled relativistic LLs in 
graphene, one is tempted to apply the same argument 
- if the Coulomb interaction e^ jeRc is sufficiently small 
as compared to the LL spacing A„, the relevant degrees 
of freedom are those which couple quantum states in the 
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(a) 



(b) 



LL 

separation 



— O — X — X — O — X — 9—0— 



Figure 23 (a) Completely filled topmost LL. Due to the Pauli 
principle, the only possible excitations are inter-LL transitions, (h) 
Partially filled LL. For sufficiently small values of Tb (or ac), the 
inter-LL excitations constitute high-energy degrees of freedom that 
may be omitted at low energies, where the relevant degrees couple 
states within the same LL. 



same LL, whereas inter-LL excitations may be considered 
as frozen out (see Fig. [^5)1 . Although this seems a rea- 
sonable assumption for the lowest LLs, one is confronted 
with the apparent problem that the LL spacing rapidly 
decreases once the Fermi level resides in higher LLs, 



A, 



(201) 



Notice, however, that this decrease is balanced by the 
l/-\/2n scaling of the characteristic Coulomb interaction, 
such that even in higher LLs the separation between low- 
and high-energy degrees of freedom is governed by the 
dimensionless coupling constant 



e^/elBV2n _ 



(202) 



which coincides with the scale-invariant zero-field cou- 
pling constant (|133p . From the interaction point of view, 
the restriction of the electron dynamics to a single par- 
tially filled LL in the large- rt limit is therefore as justified 
as for the lowest relativistic LLs. Naturally, this state- 
ment only holds true in the absence of disorder that in- 
duces stronger LL mixing for n ^ 1 than in n = or 
±1. 



describe the density fluctuations inside the LL Xn that 
interest us here. Notice that we have dropped the index 
An in the definition of the projected density operators; 
they s atisfy the quantum -mechanical commutation rela- 
tions (|Girvin et all Il986l ) 



[p(q),/5(q')] = 2isin 



q'Aq^l 



p(q + q'), (204) 



where q' A q = (q' x q)^ ~ q'^qy — q^q^ is the 2D vec- 
tor product between q' and q, and these commutation 
relations are independent of the LL index An. The in- 
formation about the LL is indeed waved to the effective 
interaction potential 



27re^ 2 
■i^rUq) = [-7>i(q)] 



(205) 



in terms of the LL form factors [see Eq. (|14ip and their 
exphcit form (jA2p . discussed in Appendix [X] 



(1 ~ Snfi)Ln-l 



2;2 



qH 



g-g^;|/4^ (206) 



independent of the band index A ( Goerbig et all . l2006t 
iNomura and MacDonaldl . I2006D . The Hamiltonian re- 
sulting from Eq. (|15ip 



(207) 



therefore defines, together with the commutation rela- 
tion (|204p the model of strongly-correlated electrons re- 
stricted to a single relativistic LL. The model respects the 
SU(4) spin- valley symmetry, and naturally, there is no 
kinetic energy scale because all processes involve states 
within the same LL. 



1. SU(4)-sym metric model 

Formally, the above-mentioned separation into high- 
and low-energy degrees of freedom may be realized with 
the help of the reduced density operators (|140p . For the 
moment, we only consider the case where ^ = i.e. 
we concentrate on the valley-symmetric model, in which 
case the reduced (intra-valley) density operators (|153p 
fall into two distinct classes: for ?i 7^ n' or A 7^ A', the 
operators PAri,A'n'(q) describe density fluctuations corre- 
sponding to inter-LL transitions of an energy equal to or 
larger than the LL separation A„, whereas the projected 
density operators 

P(q) = PAri,Ari(q) (203) 

= II ("^h""'''^|™')4„,™;C,<xCA„,,n';C,a 

5=± o'=t,i m.m' 



a. Algebraic properties. The SU(4) spin- valley symmetry 
is formally described with the help of the spin and valley- 
pseudospin operators 

5^(q) = (5^® l)®p(q) and /^(q) = (1 ® /'') ® p(q), 

(208) 

respectively, that are tensor products between the pro- 
jected density operators (|203p and the operators S"^ 
and which are (up to a factor 1/2) Pauli matri- 
ces and that describe the spin and valley degrees of 
freedom, respectively. The operators (S"^ (g) 1) and 
(1 (g) are the generators of the SU(2)xSU(2) sub- 
group of SU(4). However, once combined in a ten- 
sor product with the projected density operators /o(q), 
the SU(2)xSU(2)-extcnded magnetic translation group 
is no longer closed due to the non-commutativity of 
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the Fourier components of the projected density opera- 
tors. The commutators [S^{q),I'^{c^)] yield the remain- 
ing generators of the SU(4 )- exten d ed magnetic trans- 
lation group (|Doucot et al. I, l2008t 1 zawa and Haseb j . 
I2OO2I: lEzawa et all \200A 

Physically, the operators introduced in Eq. ()208p play 



the role of projected spin and valley-pseudospin densi- 
ties, where the LL projection is induced by the projected 
charge-density operator p(q). Their non-commutativity 
with the projected charge density, [5'^(q), p(q')] 7^ and 
[/^(q), p(q')] 7^ 0, which are due to the commutation re- 
lation (|204p . is at the origin of the (pseudo-)spin-charge 
entanglement in quantum Hall systems: as we discuss 
in more detail in Sec. IV.B.2| this entanglement yields 
(pseudo-)spin-texture states that carry an electric in ad- 
dition to their topological charge. 



b. Validity of the model. With the help of the Hamilto- 
nian (j207p . we may render more transparent the model 
assumption of electrons restricted to a single relativistic 
LL. We need to show that the energy scale that governs 
the model (|207p and its resulting phases is indeed given 
by e^/eRc and not /eIb- As an upper bound for the 
energy scale, one may use the energy of a completely filled 



the mean- 



LL described by (4«,m;C,crCAn,m';?,a-) = 5„ 

field energy (i?n) /N of which is simply the exchange en- 
ergy,28 



2e 



dg[J■„(g)]^ (209) 



In order to estimate the i ntegral in the large-n lim i t, one 
may use the scaling forni ( Abramowitz and Stegunl . ll970l : 
iGradshtevn and Rvzhikl . 120001 ) of the Laguerre polvno- 
mials 

L„ f ^ j e-«''-/^ ^ Jo(g/i3V2^^), (210) 

in terms of the Bessel function Jo(a;), such that one 
obtains by a simple change of the integration variable 
!^dq[Fn{q)f ^ {IbV2^)-^ !^ dx[Mx)Y = c/^sv^, 
where c is a numerical factor of order one. The exchange 
energy of a completely filled LL n therefore scales with 
n 1 as 



El 



sRc 



(211) 



in agreement with the model assumption of a separation 
between high- and low-energy degrees of freedom and the 
definition (|202[) of the coupling constant Uq. 



The direct energy is compe nsated by the positively charged back- 



1 ne airect energy is compe nsatea Py tne pi 
ground ("jeUium model") l lMahanl . Il993l ) 



2. Symmetry-breaking long-range terms 

When decomposing the Coulomb interaction in the 
two-spinor basis (Sec. IIII.A[) . we have seen that the 
SU(4)-symmetric model yields the leading energy scale, 
whereas the only relevant symmetry-breaking term is as- 
sociated with backscattering processes at an energy scale 
roughly a/Is times smaller than the leading one. When 
restricted to a single relativistic LL An, these backscat- 
tering terms yield a contribution 

= \ E E<'(q)'^^'"^(-q)^"^'^(q)' (212) 

?=± q 

in terms of the effective backscattering potential 



<''(q) 



J": 



An, An 



(q) 



eq 

2^e2 (1 - <5o,„) 



eq 



2n 



Kyfll 



(213) 



2;2 



where we have made use of Eq. (|142p and of the explicit 
expressions for the intervening matrix elements (|A4p . 

The effect of this symmetry-breaking term will be dis- 
cussed in more detail in Sec. IV. Bl in the context of the 
SU(4) quantum Hall ferromagnetism. The term (|212|1 is 
only relevant in relativistic LLs n 7^ as a consequence 
of the factor (1 — (5„.o) in the expres sion (I213P for the 
backscattering potential ([Goerbig et al.. . ,2006. ) . This is a 
consequence of t he ch iral symmetry of the zero-energy LL 
(|Arikawa et all l2008^ where the sublattice index is con- 
funded with the valley pseudospin, as may be seen from 
the expression ([55]) for the associated wave functions. 
Notice, however, that there may occur other symmetry- 
breaking terms in n as a co nsequence of short - range 
interactions on the lattice scale (lAhcea and Fisheil . [200a 
iDoretto and Morals Smithl . 120071: iHerbutl . 1200751 



3. Qualitative expectations for correlated electron phases 

The model of interacting electrons in a single rela- 
tivistic LL has the same structure as the one for non- 
relativistic LLs ~ in both cases, one has an interaction 
Hamiltonian that is quadratic in the projected density 
operators (|203|1 which satisfy the commutation relations 
P04|) . This is a noticeable result in the sense that, 
whereas non-relativistic 2D electron systems are gov- 
erned by Galileian invariance, the electrons in graphene 
are embedded in a Lorentz-invariant "space-time" . How- 
ever, once restricted to a single LL, the electrons forget 
about their original spatial symmetry properties and are 
governed by the magnetic translation algebra, which is at 
the origin of the commutation relations (I204p . As a con- 
sequence and in contrast to the IQHE, the differences 
between strongly-correlated electrons in graphene and 
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(b) 

r= E 
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n=0 (relativislic and non-relativistic) 
^ n=l (relativistic) 

- n=I (non-relativistic) 
n=5 (relativistic) 




n=0 

n=1 (rel.) 
n=1 (non-rel.) 




Figure 24 (a) From Ref. (iGoerbig et all . 120061 '): comparaison 
between the relativistic (black curves) and non-relativistic (grey 
curves) potentials for the LLs n = 0, 1, and 5 in real space. The 
dashed line shows the potential in n = 0, which is the same in 
both the relativistic and the non-relativistic case, (b) Pseudopo- 
tentials for n = (black circles), n = 1 relativistic (blue) and 
n = 1 non-relativistic (green). The lines are a guide to the eye. 
The open circles represent pseudopotentials with even relative pair 
angular momentum that are irrelevant in the case of completely 
spin-valley polarized electronic states. The energies are given in 
units of e^/elB- 



non-relativistic 2D electrons do not stem from their re- 
spective space-time properties, as one wo uld expect frorn 
a mean-field Chern-Si mons approach ( Khveshchenkol 
120071 : iPeres et al\ . \200(ii . 

The differences between graphcnc and non-relativistic 
2D electrons are rather to be scekcd in the larger inter- 
nal symmetry - instead of an SU(2) spin symmetry, one 
has an SU(4) spin-valley symmetry if one neglects the 
small symmetry-breaking term (|212[) in the interactions. 
Another difference arises from the different effective in- 
teraction potential (|205|) . instead of 



,,non— rcl, 



q) 



eq 



q 



(214) 



for the usual 2D electron gas. As one may see from the 
graphcnc form factors (|206p . the effective interaction po- 
tential in graphene for n ^ is the average of the non- 
relativistic ones in the adjacent LLs n and n—1, whereas 
for n — there is no difference between the relativistic 
and the non-relativistic case (see Fig. [24|) . as a conse- 
quence of the above-mentioned chiral properties. 

One notices, furthermore, that the difference between 
the relativistic and non-relativistic effective interaction 
potentials become less prominent in the large-n limit [see 
n = 5 in Fig. I24f a)]. This may be understood from the 
approximate expression (|210p of the form factors, which 
yields J-„(q) ~ [JoiqlsV^^^) + JoiqlBV2^^)]/2 ~ 



J{qlBv2^) + 0{\/n). This result agrees indeed to lead- 
ing order in 1/n with the scaling expression of the form 
factors (|154p for usual non-relativistic 2D electrons. 

The strongest difference in the interaction potentials is 
thus found for n = 1 , which in graphene is quite reminis- 
cent to that in ?? = 0, apart from a reduced repulsion at 
very short distances, whereas for non-relativistic 2D elec- 
trons it has an additional structure [see Fig. [Mf a)]. The 
behavior of the effective interaction potential may also 
be analyzed wi th the help of Haldane's pseudopotentials 
(jHaldand . [1983') 



^^w«(q)i£(<z'^|)e 



(215) 



which represent the interaction between pairs of elec- 
trons in a magnetic field, in a relative angular momentum 
state with quantum number I. This quantum number is 
related to the average distance /bv27+1 between the 
two particles constituting the pair and is a good quan- 
tum number for any two-particle interaction potential 
v{vi — Vj). The pseudopotentials for graphene are shown 
in Fig. ^h) forn = and I. 

Haldane's pseudopotentials are an extremely helpful 
quantity in the understanding of the possible FQHE 
states which one may expect in 2D electron systems. One 
notices first that as a consequence of the anti-symmetry 
of a two-particle wave function under fcrmion exchange, 
the relative angular-momentum quantum number £ must 
be an odd integer, i.e. only the pseudopotentials with 
odd values of £ play a physical role in the description 
of two interacting electrons of the same type (spin or 
valley). Even-.^ pseudopotentials become relevant if the 
SU(4) spin- valley pseudospin is not completely polarized, 
in the treatment of two electrons with different internal 
quantum number a or ^. One then notices that the n = 1 
pseudopotentials, apart from a the difference in Vi^o-, are 
much more reminiscent of those in ri = than of those 
for non-relativistic 2D electrons in the same LL n = 1 
[see Fig. [MTb)]. If one considers polarized electrons, one 
therefore expects essentially the same strongly-correlated 
electronic p hases in graphene for n = I as for n = 
( Goerbig et a i. .200^, as als o corroborated by numerical 
studi es for FQHE states (lApalkov and Chakrabortv , 



2009: 



2006 : Goerbig and Regnaultl. 120071: iPapic et al. 



Toke and Ja inl. l2007t iToke et all 



2006) 



2009; 



and electron-solid phase s (IPoplavskw et all 
IZhang and Joglekail I2007L l2008l) . Because the pseu- 
dopotentials (|215p are systematically larger in ti = 1 
than in 71 = (apart from the short-range component for 
£ ~ Q), the gaps of the FQHE states in n = 1 are larger 
than the corresponding ones i n rt = 0, as one may also see 
from numerical calculations ( Apalkov and Chakrabortvl 
[2006; Toke et ai, 2002,). 

As much as we have emphasized the similarity between 
the n = and n = 1 LLs in graphene, we need to stress 
the difference between the n = 1 LL in graphene as 
compared to n = 1 in non-relativistic 2D electron sys- 
tems. Remember that the physical phase diagram in the 
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non-relativistic n = 1 LL is extremely rich; an intrigu- 
ing e ven-denomina t or FQ HE has been observed at v = 
5/2 (jWillett et all Il987t ) and pro bably possesse s non- 
Abelian quasiparticle excitations ( Greiter et all Il99ll : 
Im oore and Readl . Il99l[) . Furthermore, a particular com- 
petition between FQHE states and electron-solid phases, 
which is chara cteristic of the non-relativistic n = 1 LL 
(|Goerbig et al., . 2003 , 2004) , is at the origin of the reen- 



trance phenomena observed in tran s port m easurements 
(|Eisenstein et al\ . 120021 : iLewis et all boOSh . These phe- 



nomena are absent in the n = LL, and the similarity 
between the n = and the relativistic n = \ LL thus 
leads to the expectation that FQHE states corresponding 
to the 5/2 state in non-relativistic quantum Hall systems 
and the above-mentioned reentrance phenomena are ab- 
sent in the n — 1 LL in graphene. This expectation 
has recently been corroborated in ex act-diagona l izatiq n 
studies on the non- Abelian 5/2 state (|Wois et adl2Q10| ). 



is more involved because there is no physical field that 
couples directly to the valley pseudospin, as suggested by 
the otherwise intuitive form (|217p . There have however 
been proposals that such an effect may be achieved with 
the help of strain- induced disordered gauge fields that 
mimic large-scale ripples ( Mever et g/J . 120071) and that 



yield a n easy-plane anisotropy in n = (jAbanin et al\ . 
l2007al) . similarly to the backscattering term (|213p in 
higher LLs. Quite generally, a valley-degeneracy lifting 
may be achieved indirectly in the zero-energy LL n = 
via fields that couple to the sublattice index. This is due 
to the fact that the components -0^ ^ and ip'^, ^ vanish 
as a consequence of the chiral properties, which identify 
the sublattice and the valley pseudospins in rt = 0, as we 
have discussed in Sec. IH.A.2I 

In order to illustrate this indire ct lifting of the valley 
degeneracy, we consider the term ( Haldang Il988l) 



Mh = ^/f^ (1 



valley ® '''AB ' 



0*, 



(218) 



4. External spin-valley symmetry breaking terms 

Before we consider the different phases due to electron- 
electron interactions, we start with an analysis of the 
different external effects, which are capable of lifting 
the four-fold spin- valley degeneracy. 

The probably most familiar external symmetry- 
breaking term is the Zeeman effect, which lifts the two- 
fold spin degeneracy while maintaining the SU(2) symme- 
try associated with the valley pseudospin. The size of the 
Zeeman splitting is given by the energy Az ~ 1.2i?[T] K, 
for a q -factor that has bee n experimentally determined as 
5 ~ 2 (|Zhang et al ]. 120061 ). If we adopt a compact eight- 
spinor notation to take into account the four different 
spin- valley components, in addition to the two sublattice 
components, the Zeeman term has the form 



A spin 



^'^ (1 valley ® IaB ® <pin) 

+ {A^B), 



(216) 



where the tensor product consists of the valley pseu- 
dospin (represented by the Pauli matrices T^aiicy ^^'^ 
lvaiicy)i the sublattice pseudospin {r^g and Iab), and 
the true spin {rl^p^^ and Ispin)- For a better understand- 
ing, we have given the explicit expression in terms of 
spinor components in the second line of Eq. pi6|) . 

A possible valley-degeneracy lifting, that one could de- 
scribe with the help of a "valley Zeeman effect" similarly 
to Eq. (I2T6)1 . 

AT"' ~ (Galley ® ^AB ® Ispin) (217) 



which breaks the lattice inversion symmetry and opens a 
mass gap at the Dirac point in the absence of a magnetic 
field. In the presence of a S-ficld, the LL spectrum 
is modified by the term (|218p and reads 



e\.n;i - AW A/2 



n. 



(219) 



for n 7^ 0, independent of the valley index ^, whereas the 
fate of the n = LL depends explicitly on ^, 



(220) 



such that the valley degeneracy is effectively lifted. No- 
tice, however, that due to the vanishing components ■0^ ^ 

and ^, the mass term (|218p is now indistinguishable 
(in n = 0) from the above-mentioned valley Zeeman term 

(EHi, 



(1 valley ® l.pin)^' 

«alley ® l^B «> Ispin) (221) 



whereas this is not the case for the LLs n 7^ 0, where 
the valley degeneracy is only lifted by an explicit val- 
ley Zeeman effect. A mass term of the form (|218|) 
typically arises in the presence of a frozen out-of-plane 
phonon that yields a cru mbling of the graphene sheet 
(|Fuchs and Ledereii [2007t) . 

More recent studies have concentrated on a sponta- 
neous deformation of the graphene sheet due to frozen 
i n-plane phonons that yield a Kekule -type distortion 
(|Hou et all \20ld iNomura et ail . |2009[) . This distor- 
tion, which is associated with a characteristic wave vec- 
tor 2K and which therefore couples the two valleys, di- 
rectly breaks the valley degeneracy, via a term Mk — 



My, with 



By external effects we mean tfiose tfiat are not caused by the 
mutual Coulomb repulsion between the electrons. 



M,,y = ^vl/t (^-^^^ ^ ^ VI/. (222) 
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Such a term yields the same energy spectrum (|219p and 
(|220|) as the mass term piSp if one replaces M by Afee/c/2, 

with the characteristic energy scale A^efc = 



A2 + A2 



2B[T] K, (lAiiki and Andd . Il995t IHou et all l2010l) . No- 
tice that this energy scale is slightly larger than, but 
roughly on the same order as, the Zeeman energy scale. 

Finally, we mention another class of terms that break 
the spin- valley degeneracy and that have received recent 
interest in the framework of res earch on topological in- 
sulators [for recent r eviews see ( Hasan and Kand . I2OIOI : 
iQi and Zhan"g| . l201lh ]. In an original work, Haldane ar- 
gued that a time-reversal-symmctry breaking term with 
an inhomogeneous flux distribution inside each hexagon 
opens a gap in a hon eycomb lattice with zero magnetic 
field ( Haldan j . Il988l ). Most saliently, he showed that 
one may thus achieve a quantum Hall effect without an 
external magnetic field, a system that is now often re- 
ferred to as the "quan t um anomalous Hall insulator" 
(iHasan and Kand . I2OIOI: iQi and ZhanS . I2OIII) . A simi- 
lar situation arises when spin-orbit interactions are taken 
into account, which are of the form 



Uso - 



A 



so 



valley 



' AB 



spiny 



(223) 



and which provide again the same LL spectrum (1219 1 
and ([220]) if one replaces M by Aso/2 (|Kane and Mele . 
[2005). In spite of the conceptually appealing prospect 
of the quantum spin Hall effect, which is revealed by 
this model because the spin orientation is locked to a 
particular valley index via the term (|223p . the associ- 
ated energy scale A50 ~ 10 mK turns out to be van- 



ishingly small in graphene, whereas an extrinsic Rashba- 
type sp in-orbit coupling in graphene can be on the order 
of 1 K (|Min et flll . l2006h . 



5. Hierarchy of relevant energy scales 

These energy scales associated with external fields need 
to be compared to the characteristic (bare) interaction 
energy e^/eZs ~ 625{^/B\T]/e) K, which is, for experi- 
mentally accessible magnetic fields, much larger than Az 
or Afcefe. As we have discussed in Sec. IIII.B.41 inter-band 
LL excitations screen the bare Coulomb interaction and 
yield a contribution to the dielectric constant. In the ab- 
sence of a quantizing magnetic field, we have seen that 
this dielectric constant is given by [see Eq. (|176p ] 



TT 

1 + -- 



2 hevp 



(224) 



where we remember that e is the extrinsic dielectric con- 
stant of the surrounding medium. As one may notice 
from Fig. [THl the vacuum contribution H'""^(q) (thick 
dashed lines) is only marginally modified by the magnetic 
field, such that one may use /sEooIb as an approxima- 
tion for the interaction-energy scale for graphene taking 
into account inter-band screening. 

The relevant energy scales are summarized in the table 
below for different values of the magnetic field, in com- 
parison with the interaction energy scales, taking into ac- 
count the effective dielectric constants for several widely 
used substrates from Tab. I. 



energy 


value for arbitrary B 


for B = 6 T 


for B = 25 T 




1.2B[T] K 


7 K 


30 K 




2B[T] K 


12 K 


50 K 


/eIb (bare) 


625(VB[Tl/e) K 


(1550/e) K 


(3125/e) K 


/eEocIb (vacuum) 


139a/-B[T] K 


344 K 


694 K 


/eEooIb (on Si02) 


10A^JB[T] K 


258 K 


521 K 


/eEooIb (on h-BN) 


109a/B[T] K 


270 K 


543 K 


/eEooIb (on SiC) 


7lVB[T] K 


176 K 


355 K 


^sb < [e^ /ElB){a/lB) 


< 1B[T] K 


< 6 K 


< 25 K 



Table 11 Energy scales at different magnetic fields. The first two lines show the energy scales associated with the major external symmetry- 
breaking fields (Zeeman and Kekule-type lattice distortion, and Afcefc, respectively), which scale linearly in B. Below are shown the 
interaction-energy scales (oc y/B), the bare one with an unspecified dielectric constant and the ones for different substrates taking into 
account inter-band screening via the term £00 [Eq. I|224| l]. The last line yields the interaction-energy scale associated with the intrinsic 
symmetry breaking due to inter-valley coupling, discussed in Sec. IV. A. 21 



r 



In view of the above discussion, one may conclude that yields the leading energy scale in the problem of clec- 
the SU(4)-symmetric part of the Coulomb interaction trons in partially filled lower LLs, whereas external terms. 
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such as the Zeeman effect or spontaneous lattice distor- 
tions, play a subordinate role. The remainder of this 
section is therefore concerned with a detailed discussion 
of strongly-correlated electron phases that are formed to 
minimize the Coulomb interaction. 



B. SU(4) Quantum Hall Ferromagnetism in Graphene 

A prominent example of the above-mentioned strongly- 
correlated phases is the generalized quantum Hall fcrro- 
magnet. It arises in systems with a discrete internal de- 
gree of freedom described by an SU(A/') symmetry, such 
that each single-particle quantum state ipn,m occurs in 
M copies. Prominent examples are the non-relativistic 
quantum Hall systems when taking into account the elec- 
tronic spin a =t,i (A/" = 2) or bilaycr quantum Hall sys- 
tems that consist of two parallel 2D electron gases, where 
the layer index may be viewed as a "spin" 1/2 [N = 2, 
or A/" = 4 if one also takes into account the ph ysical spin 
(|Ezawa and Hasebd . [200^ : lEzawa am2003l) ].^° In this 
sense, graphene may be viewed as an SU(4) quantum 
Hall system as a consequence of its four-fold spin-vallcy 
degeneracy. 



1. Ferromagnetic ground state and Goldstone modes 

Quite generally, quantum Hall ferromagnetism arises 
when the filling factor, defined from the bottom of 
the LL,"^^ is an int eger that is not a multiple of N 
( Arovas et al\ . Il999| ). From the point of view of the ki- 
netic Hamiltonian, one is thus confronted with a macro- 
scopic ground-state degeneracy. Even if one has an inte- 
ger filling factor, the situation is thus much more remi- 
niscent of the FQHE, i.e. the relevant energy scale is the 
Coulomb interaction, and the system may be described in 
the framework of the model (j207p of interacting electrons 
in a single (relativistic) LL. For the moment, we consider 
that there are no symmetry-breaking terms, such as the 
backscattering term (|212p or Zeeman-type terms that are 
discussed below in Sec. IV. A. 41 

Qualitatively, one may understand the formation of a 
ferromagnetic ground state as a consequence of the repul- 
sive Coulomb interaction. In order to minimize this inter- 
action, the electrons prefer to form a state described by 
a maximally anti-symmetric orbital wave function that 
must then be accompanied by a fully symmetric SU(A/') 
spin wave function to satisfy an overall fermionic (anti- 
symmetric) wave function. In a usual metal with a fi- 
nite band dispersion, this ferromagnetic ordering (e.g. all 



For a review on non-relativistic mul ti-component systems see 
Ref. llEzawal . I2OO0I : iMoon et aLl . ll995l) . 

Remember that the filling factor in graphene is defined with re- 
spect to the center of the n = LL. There is thus a shift of 2 in 
the filling factor as compared to the non-relativistic case. 



electrons in the spin-f states) is accompanied by a cost 
in kinetic energy - indeed, the Fermi energy for spin- 
t electrons is increased whereas that of spin- J, electrons 
is lowered. The competition between the gain in inter- 
action and the cost in kinetic energy defines the degree 
of polarization, i.e. how ferromagnetic the electrons ef- 
fectively are. In the quantum Hall effect, however, we 
arc confronted with a highly degenerate LL that may be 
viewed as an infinitely flat band, such that the kinetic- 
energy cost for complete spin polarization is zero. 

As an example of an SU(A/') quantum Hall ferromag- 
net, one may consider the state 

k Nb-1 

i™)=n n (225) 

i—l m—0 

which consists of fc < A/" arbitrarily chosen completely 
filled subbranchcs [i e {(A', t), (A'', t), (A', 1), (A'', 1)}, 
for the SU(A/' = 4) symmetry in graphene LLs], where we 
have omitted the LL index An at the fermion operators 
to simplify the notation. The arbirariness in the choice 
of the SU(A/') spin subbranchcs may be viewed as a spon- 
taneous symmetry breaking that accompanies the ferro- 
magnetism. Indeed, the state ()225p is no longer invariant 
under an SU(A/') rotation, but only under a rotation de- 
scribed by the subgroup SU(fc)xSU(A/' — fc), where the 
first factor indicates a symmetry transformation in the 
fully occupied subbranchcs i = 1, fc and the second fac- 
tor one in the empty subbranchcs i = fc + 1, ...,Af. There- 
fore the quantum Hall ferromagnct (|225[) is associated 
with an order parameter with a spontaneous symmetry 
breaking described by the coset space SU(A/')/SU(fc) x 
SV{Af-k) X U(l) ~ U(AA)/U(fc) X U(A^- fc), where the 
additional U(l) is due to the phase differen ce between the 
occup ied and the unoccupied subbranchcs (jArovas et al\ . 
ll999l : lYang et aLl . l2006D . 

The coset space, with its 7V^ - fc^ - (A/" - fc)^ = 
2k{N — fc) complex generators, defines also the Gold- 
stone modes, which are nothing other than the k{N — fc) 
spin-wave excitations of the ferromagnetic ground state 
(|225p .'^^ The number of the spin- wave modes may also 
have been obtained from a simple inspection into the LL- 
subbranch spectrum. Indeed, a spin wave can be de- 
scribed with the help of the components of the projected 
density operators (|203p . 

P^M) = E ("^ ^'''■''1 ™') (226) 

which represent coherent superpositions at wave vector 
q of excitations from the occupied subbranch j to the 
empty subbranch i. One has then fc possibilities for the 
choice of the initial subbranch j and Af — k for the final 



The complex generators come in by pairs of conjugate operators, 
and each pair corresponds to one mode. 
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one, and one obtains therefore k{N — k) different spin- 
wave excitations, in agreement with the above group- 
theoretical analysis. 

Notice that all spin-wave excitations have the 
same dispersion, which ma y be calculated within 

n (lAlicea and Fisher 



a mean-field 



approximati oi 
Doretto and Morais SrnitlJ, I2n07t 



2006; 



Kallin and Halperin 



MM; Yang et al, 200C 

(FM|p„-(-q)i7„p,,(q) - i7„|FM) 
2 / qAk^l^ 



2^i;„(k) sin^ 

k 



which saturates at large values oi q= |q|. 



E„ 



2E 



X 



(227) 



(228) 



i.e. at twice the value of the exchange energy p09p . This 
result is not astonishing insofar as the large-;/ limit cor- 
responds, as we have discussed in Sec. IIILB.2I [see Eq. 
(|167p ]. to an electron- hole pair where the electron is situ- 
ated far away from the hole. The energy (j228p is therefore 
nothing other than the cost in exchange energy to create 
a spin-flip excitation, i.e. an electron with reversed spin 
and a hole in the ferromagnetic ground state. Because of 
the large distance between the electron and the hole in 
such an excitation and the resulting decoupled dynamics, 
one may be tempted to view this energy as the activa- 
tion gap of the quantum-Hall state at v = k, but we see 
below in Sec. IV. B. 21 that there exist elementary charged 
excitations (skyrmions) that have, in some LLs, a lower 
energy than these spin-hip excitations. 

In the opposite limit of small wave vectors {qIb ^ 1), 
one may not understand the excitation in terms of decou- 
pled holes and electrons, and the excitation can there- 
fore not contribute to the charge transport. A Taylor 
expansion of the sine in the spin-wave dispersion p27p 
yields the usual dispersion of the spin- wave Goldstone 
modes, 



E. 



in terms of the spin stiffness 

^" = iE-"(k)|k| 



2j2 



(229) 



(230) 



One notices that the above results for the excitation en- 
ergies do not depend on the size of the internal symmetry 
group, but they can be derived withi n the SU(2) model 
of the quantum Ha ll ferromagnetism ( Moon et aLl . ll99"5t 
ISondhi et al\ . Il993() - the enhanced internal symmetry of 
graphene (or of a general A^-component system) affects 
only the degeneracies of the different modes. 



2. Skyrmions and entanglement 

In addition to the above-mentioned spin-wave modes, 
the SU(A/') ferromagnetic ground state (|225p is charac- 
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Figure 25 Excitations of the SU(2) ferromagnetic state, (a) Spin 
waves. Such an excitation can be continuously deformed into the 
ferromagnetic ground state [spin represented on the Bloch sphere 
(right) by fat arrow] - the grey curve can be shrinked into a sin- 
gle p oint (b ) Skyrmion with non-zero topological charge [from Ref. 
l lGirvinl . 119991) ]. The excitation consists of a a reversed spin at 
the origin z = 0, and the ferromagnetic state is recovered at large 
distances [z/CI c«. Contrary to spin-wave excitations, the spin 
explores the whole surface of the Bloch sphere and cannot be trans- 
formed by a continuous deformation into the majority spin (fat 
arrow) . 



terized by a particular elementary excitation that con- 
sists of a topo l ogical spin texture, the so-called skyrmion 
(jSondhi et all Il993[) . Similarly to a spin-wave in the 
limit qlg 1, the variation of the spin texture in a 
skyrmion excitation is small on the scale of the magnetic 
length I B , such that its energy is determined by the spin 
stiffness (|230p in the small-g l imit. Indeed, o ne may show 
that its energy is given by ( Ezawal l2000t iMoon et all 
ll995l : ISondhi et aLl . ll993l) 



top 1 1 



(231) 



in terms of the topological charge Qtop, which may be 
viewed as the number of times the Bloch sphere is covered 
by the spin texture [see Fig. HSfb)] and which we discuss 
in more detail below. Skyrmions are the relevant elemen- 
tary excitations of the quantum Hall ferromagnetism if 
the energy (|23ip is lower than that of an added electron 
(or hole) with reversed spin, which is nothing other than 
the exchange energy (|209p . i.e. if Esk < E^. Whereas 
this condition is fulfilled in non-relativistic LLs only in 
the lowest one n = 0, skyrmions are the lowest-energy 
el ementary excitations in the graphen e LLs n = 0, 1 and 
2 (|Toke et ad 120061: lYang et ai.l . l2006l ). as a consequence 
of the difference in the form factors. 

As in the case of the spin waves discussed above, the 
skyrmion energy is independent of the size of the internal 
symmetry group, and we will first illustrate the skyrmion 
texture in an effective SU(2) model, where the texture is 
formed only from states within the last occupied (fc) and 
the first unoccupied (fc-|-l) LL subbranch. The skyrmion 
may then be described with the help of the wave function. 
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in terms of the complex coordinate z = (x — iy)/lB, 



k,k+l) — 



+ z 



(232) 

where | f'^ {z)) corresponds to states in the subbranch k 
and I i''^^ {z)) to those in fc + 1, at the position z. One 
notices that at the origin z = the "spin" associated with 
these two components is ^ because the first component 
of Eq. ()232p vanishes, whereas the spins are t at jz/CI — > 
oo (see Fig. [^5]) . where the ferromagnetic ground state 
is then recovered. The parameter ( plays the role of 
the skyrmion size, measured in units of Ib ^ indeed, for 
\z\ = Id, both components are of the same weight and 
the spin is therefore oriented in the xy-plane. 

The skyrmion excitation (|232p can also be illustrated 
on the so-called Bloch sphere on the surface of which the 
(normalized) spin moves (see Fig. [25|). The angles {9 
for the azimuthal and (f) for the polar angle) of the spin 
orientation on the Bloch sphere correspond to the SU(2) 
parametrization lip) = cos(6'/2)| f) +sin(0/2) exp(i0)| \.), 
and the spin orientation at the circle \z\ = |C| in the com- 
plex plane describes the equator of the Bloch sphere. The 
topology of the skyrmion excitation becomes apparent 
by the number of full circles the spin draws when going 
around the origin of the xy-plane on the circle \z\ = \(\. 
More precisely, the topological charge Qtop is not defined 
in terms of such closed paths, but it is the number of full 
coverings of the Bloch sphere in a skyrmion excitation 
[Qtop = 1 in the example (|232p ]. Notice that a spin- wave 



excitation has a topological charge Q 



top 



and corre- 



sponds to an excursion of the spin on the Bloch sphere 
that is not fully covered and that can then be reduced 
continuously to a single point describing the ferromag- 
netic ground state [Fig. I^ST a)]. 

The above considerations may be generalized to sys- 
tems with larger internal symmetries, i.e. to SU(A/') 
quantum Hall ferromagnets. The state ([232]) is invariant 
under the SU(7V) subgroup SU(fc - 1) x SU(A/' -k-1), 
where the first factor describes a rotation of the occupied 
subbranches that do not take part in the skyrmion exci- 
tation and the second factor is associated with a symme- 
try transformation of the corresponding unoccupied sub- 
branches k + 2, ...,Af. A similar group-theoretical analy- 
sis as the one presented in Sec. IV.B.llyields the number 
of residual symmetry transformations ( Yang et aLl . l200"6h 
2k{N — k) + 2{J\f ~ 1), where the first term describes the 
Goldstone modes of the ferromagnetic ground state, and 
the second one corresponds to the A/" — 1 internal modes 
of the skyrmion excitation. 

In addition to the topological charge, skyrmions in 
quantum Hall systems carry an electric charge that co- 
incides, for V = k with the topological charge. Indeed, 
the skyrmion state (|232p describes an electron that is ex- 
pelled from the origin z = in the f-component, and its 
net electric charge is therefore that of a hole. This means 
that skyrmions are excited when sweeping the filling fac- 
tor away from v = k, and the net topological charge is 



given by Qtot = \v~k\NB- The number of internal modes 
is then Qtot{J^ ~ l)i in addition to one mode per charge 
that correspon ds to a simple trans lation z — >■ z + a of 
the excitation (jPoucot et al\ . |2008| ). As a consequence 
of the Coulomb repulsion, it is energetically favorable to 
form a state in which Qtot skyrmions of charge 1 are ho- 
mogeneously distributed over the 2D plane than a single 
defect with charge Qtot (|Moon et al\ . Il995l ) . A natural 
(semi-classical) candidate for the gro und state of Qto t 
skyrmions is then a skyrmion crystal ( Brev e^"all . Il995l ) 
that has recently been revisit ed in the framework of the 
SU(4) symmetry in graphene (jCote et al. I. l2007ll2008[) . In 
this case, the Qtot{-N' — 1) internal modes, which are dis- 
persionless zero-energy modes in the absence of electronic 
interactions or Zeeman-type symmetry-breaking terms, 
are expected to yield M — 1 Bloch bands of Goldstone- 
typc, in addition to the Qtot translation modes that form 
a magnetic-field phonon mode of the skyr mion crystal 
with a characteristic uj oc g^/^ dispersion ( Fukuvainal . 



a. Skyrmions and activation gaps in graphene. Quite gen- 
erally, the activation gap in quantum Hall states is the 
energy required to create a quasi-particle-quasi-holc pair, 
in which the two partners are sufficiently well separated 
to contribute independently to the charge transport. In 
the framework of the quantum Hall ferromagnet, the 
activation gap may be viewed as the energy to create 
a skyrmion of topological charge Q = \ and an anti- 
skyrmion of charge Q = — 1 that are well-separated from 
each other such that one may neglect their residual in- 
teraction. The energy of such a skyrmion-anti-skyrmion 
pair is then given, in the absence of symmetry-breaking 
terms, by twice the energy in Eq. (|23ip . 



(233) 



For graphene, the energies of the theoretical activation 
gaps for n — Q and n = 1 are shown in the table below. 





activation gap 


arbitrary value of B 


n = 


1 pK <•■' 

2 V 2 cEooiH 


400(^5 [T]/eeoc) K 


n = 1 


7 /IT 

32 V 2 EEoois 


175(^5 [T]/eeoc) K 



Table III Theoretical estimates for the activation gaps in the 
n = and 1 graphene LLs due to well-separated skyrmion-anti- 
skyrmion pairs. 



For further illustration, we consider the scenario in 
which the Zceman effect is the only SU(4)-symmetry 
breaking term.'^'^ Due to the Zeeman effect, spin- J, elec- 



The energetic argument remains valid in the case where the dom- 
inant term is a valley-Zeeman effect if one interchanges the role 
of the spin and the valley pseudospin. 
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trons are energetically favored. If only one spin- valley 
branch of a particular LL is filled (fc = 1), the spin mag- 
netization of the spin- valley ferromagnet is preferentially 
oriented in this direction whereas the valley polarization 
may point in any direction. The activation gap would 
then be dominated by valley (anti-)skyrmions with no 
reversed physical spin such that one would not expect 
any dependence of the gap on the in-plane component of 
the mag netic field, in agreem ent with the experimental 
findings (| Jiang et adl2007b[ ). 

The situation is different when both valley branches of 
the spin- J, branch are occupied; an excitation of the SU(4) 
ferromagnet with a full spin polarization would then nec- 
essarily comprise reversed spins, and the corresponding 
Zeeman energy must be taken into account in the energy 
of the (spin) skyrmion-anti-skyrmion pair p33|l . 



(234) 



where Nrs ~ |CP is the number of reversed spins in a sin- 
gle (anti-)skyrmion. Notice that this number depends on 
the competition between the Zeeman effect itself, which 
tries to reduce the skyrmion size and the cost in ex- 
chang e energy due to the strong variation in small tex- 
tures (jMoon ad Il995l : ISondhi ad Il993[ ).^'^ The en- 
ergy of a skyrmion-anti-skyrmion pair in the spin channel 
(with two completely filled valley sublevels) is therefore 
larger than that p33p of a pair in the valley channel 
when only one valley subbranch of the LL is completely 
filled. Notice that this energy increase may even be sig- 
nificant for large skyrmions because of the larger num- 
ber of reversed spins. As a thumbrule, the stability of a 
quantum Hall state is proportional to the activation gap, 
which has in the present case been identified with the 
skyrmion-anti-skyrmion energy and which is dominated 
by the Coulomb interaction energy. Additional external 
symmetry-breaking terms, such as those discussed in Sec. 
IV. A. 41 may enhance this stability although they provide 
only a small correction to the activation energy. 




Figure 26 From Ref. l lDoucot et all l2008t) : Bloch spheres for 
entangled spin-pseudospin systems. Bloch sphere for the spin 
(a), pseudospin (b), and a third type of spin representing the 
entanglement (c). In the case of spin-pseudospin entanglement 
(|coso| 1), the {pseudo)spin-magnetizations explore the interior 
of their spheres, respectively (black arrows). 



four-spinor 



\^{z)) ^ cos -\^Ps)\4'i) 



sin-e*'^|xs)|x/> 
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where a and /? are functions of the complex position z, 
and the local two-component spinors I'ips), \xs), IV'/)) 
and \xi) arc constructed according to 



sin |e"^ 



and 



\x) = 



The angles 6 and (f) define the usual unit vector 
n{6, (f)) = (sin 9 cos cj), sin 9 sin (f>, cos 9) , 



(236) 



(237) 



b. Spin-valley entanglement in graphene. In an experimen- 
tal measurement, one typically has not direct access to 
the full SU(4) spin that describes the internal degrees 
of freedom in graphene LLs, but only to the SU(2) part 
associated with the physical spin, e.g. in a magnetiza- 
tion measurement. It is therefore useful to parametrise 
the SU(4) spin in a manner such as to keep track of the 
two SU(2) copies associated with the physical spin and 
with the valley pseudospin, respectively. This may be 
achieved with the so-called Schmidt decomposition of the 



which explores the surface of the Bloch sphere depicted 
in Fig. [251 Notice that one has two Bloch spheres, one for 
the unit vector n{Os, '/>s) associated with the spin angles 
6s and <f)s and a second one for n{9j, (pi) for the valley- 
pseudospin angles 9i and 0/ (see Fig. [26]) . In addition, 
one may associate a third Bloch sphere with the angles 
a and (3 that describe the degree of factorizability of the 
wave functions and thus the degree of ent anglement be- 
tween the spin and the valley pseudospin ( Doucot et al\ . 
I2008D . 

With the help of the Schmidt decomposition (|235p . one 
obtains immediately the reduced density matrices for the 
spin and the valley-pseudospin sectors 



This energy cost may be evaluated from a gradient expansion of PS Tr/ (|^')(5'|) - 
the energy in the magnetization fields. At leading order, one ob- 
tains, however, a non-linear sigma model that is scale-invariant, pj ~ Tr5 (|^')('I'|) 
such that the ener gy cost must be calculated at higher orders 
llMoon et aLl.ll995l) . 
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respectively, and the local spin and valley-pseudospin 
densities are simply 

= Tr(ps^") = cosa(^5|^1^s) = cos an'' {9 s,4's) 

(239) 

and 

= Tr(p//^) = cosa(i/'/|/''|V'/) = cos a 71^(6*/, 0/), 

(240) 

where S*" and represent the components of the spin 
and valley-pseudospin operators, respectively [see Eq. 
(|208p ]. One notices from these expressions that, for the 
case a 7^ or TT (i.e. cos^ a < 1), the local (pseudo)spin 
densities are no longer normalized, but they are of length 
|ms/7p = cos^ a. Thus, in a semiclassical picture, the 
(pseudo)spin dynamics is no longer restricted to the sur- 
face of the Bloch sphere, but explo res the entire volum e 
enclosed by the sphere (Fig. ^ (jPoucot et all |2008( ). 
This result indicates that one may be confronted, in the 
case of full entanglement (e.g. a = tt/2), with an SU(4) 
quantum Hall fcrromagnct the (spin) magnetization of 
which completely vanishes, as one would naively expect 
for an unpolarized state. 



3. Comparison with magnetic catalysis 

An alternative scenario proposed for the deg eneracy 
lifting in n = is that of t he magnetic ca t alysis (lEzawa , 



2007 



2007b 



Gorbar et all . 120081: iGusvnin et "oil . l2006t iHerbut 
20081 ). which was discussed even before the dis- 



covery of graphene ( Gorbar et all . 120021 : iKhveshclicnkoi 
l200l[ ). According to this scheme, the Coulomb interac- 
tion spontaneously generates a mass term for the (origi- 
nally masslcss) 2D electrons once the magnetic field in- 
creases the density of states at zero energy by the for- 
mation of the highly degenerate n = LL. As a con- 
sequence of this mass generation, the particles condense 
in a state of coherent particle-hole pairs (excitonic con- 
densation). The effect is at first sight reminiscent of the 
excitonic condensation at = 1 in non-relativistic bi- 
laycr quantu m Hall sy stems ("Ezaw a and Iwazakil Il993l : 
iFcrtig , 1989 : Wen and Zee . , 1992a) . Its superfiuid behav- 
ior gives rise to a zero-bias anomal y in the tunneling con - 
ductance between the two layers ( Spielman et aLl . l2000f ) 



as well as to a simultaneous suppression of the longitudi- 
nal and the Hall res i stance in a counte rflow experiment 
(|Kellogg et adl2004l : lTiriuc et a/.l . l2004l) . The bilayer ex- 
citonic condensate may be described as an ea sy-plane 
quantum Hall ferromagnet ( Moon et all Il995l ). where 
the spin mimics the layer index. The origin of this easy- 
plane anisotropy stems from the difference in the inter- 
actions between electrons in the same layer as compared 
to the weaker one for electron pairs in different layers. 

This comparison with non-relativistic 2D electrons in 
bilayer systems indicates that there may exist a close re- 
lation between the quantum Hall ferromagnetism and the 
scenario of the magnetic catalysis also in graphene in a 
strong magnetic field. Notice, however, that the excitonic 



state in graphene is not in the same universality class as 
that of the quantum-Hall bilayer - in the latter case the 
symmetry of the (interaction) Hamiltonian is U(l) as a 
consequence of the easy-plane anisotropy, and the sym- 
metry breaking is associated with a superfiuid mode that 
disperses linearly with the wave vector, uj q. In con- 
trast to this system, the interaction Hamiltonian (|207p 
has the full SU(4) symmetry, and even for a sufficiently 
strong Zeeman effect, the symmetry is quite large with 
SU(2)^ X SU(2)^, i.e. each spin projection f and | is gov- 
erned by the residual SU(2) valley symmetry and has the 
characteristic a; cx pseudospin-wave modes. 

The connection between the two scenarios becomes 
transparent within a mean-field treatment of the 
Coulomb interaction Hamiltonian. The quantum Hall 
ferromagnetic states discussed in the previous subsec- 
tions may be described equivalently with the help of the 
mean-field order parameters 



valley 



spin ' 



(241) 



where denotes the same eight-spinor as in Sec. IV. A. 41 
Remember that a pure spin quantum Hall ferromagnet 
is obtained for T^^aHoy = 1 valley; whereas a pure valley- 
pseudospin ferromagnet is described by an order parame- 



ter ((24T|) with ■ 



spin 



The remaining order param- 



eters describe states with a certain degree of spin-valley 
entanglement, as discussed above. 

Notice, however, that the choice of order parameters 
is not restricted to those in Eq. (|241[) . Indeed, one may 
also opt for a mean-field calculation of t he interaction 
Ham i ltonian with the ord er parameters ( Gorbar et al\ . 
l2008l : IGusviAT" et al. I l2006l) 



and 



Ms = (' 



valley 



valley ' 



'AB ' 



• AB ' 
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which describe mass gaps. Indeed, we already encoun- 
tered a term of the form ()243p in Sec. IV. A. 41 and showed 
that it lifts the valley degeneracy of the n = LL. 
Whereas such a term arises naturally in the context of 
an out-of-plane distortion of the graphene lattice, here, 
it is generated dynamically via the repulsive electron- 
electron interaction. The difference between the two 
mass terms Ms and Mt stems from the residual sym- 
metry of the SU(2)ct groups. The term (|242p . which may 
be viewed as a singlet mass term explicitly breaks this 
symm etry, whereas the term (l243l) has been coin ed triplet 
mass ( Gorbar et al I l2008HGusvnin et aLl . [2006h . 

In Sec. IV.A.4I we have argued that mass terms of the 
above form only lift the valley degeneracy in the zero- 
energy LL n = 0, whereas they simply renormalize the 
LL energy for n ^ 0. Furthermore we have seen that as 
a consequence of the vanishing spinor components ip"^ ^ 



and Tp 



the mass term Mt is indistinguishable, in 
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0, from a valley-pseudospin ferromagnetic state, 



(Ivalloy ® T^B » Ispin) *) 



(Galley ® l^B ® Ispin) *) , (244) 



whereas the singlet mass term simply renormalizes the 
overall chemical potential, 



M 



n=0 



(vl/t 



^valley ® '''AB ' 



0*) 



(Ivalloy IaB 



*) . (245) 



These arguments lead to the conclusion that the mag- 
netic catalysis in n = 0, i.e. the spontaneous genera- 
tion of a mass gap due to electron-electron interactions, 
may be fully described in the framework of the SU(4) 
quantum Hall ferromagnetism. Furthermore, the recent 
observation of a ful ly lifted spin- valle y degeneracy in the 
n = 1 graphene LL ( Dean et all . \201)h is naturally under- 
stood in the framework of quantum Hall ferromagnetism, 
whereas the mass terms ()242p and (j243p . obtained from 
magnetic catalysis, would not provide a fully lifted spin- 
valley degeneracy. 



4. The quantum Hall effect at u = ±1 and = 

Before discussing the experimental results on the quan- 
tum Hall effect, a clarification on the filling factor is re- 
quired. In the preceding parts of this section, which were 
concerned with general aspects of the quantum Hall fer- 
romagnet in LLs with internal degrees of freedom, the 
filling factor v = k has been defined with respect to the 
bottom of the partially filled LL. However, in graphene, 
this is at odds with the natural definition of the filling 
factor (|10ip in terms of the electronic density measured 
from the charge neutrality point in undoped graphene - a 
zero filling factor therefore corresponds to two completely 
filled spin- valley subbranches {k ~ 2) oi the n — LL. In 
the remainder of this section, we therefore make a clear 
distinction between the two filling factors, and v denotes 
the filling of the n = LL measured from the bottom of 
the level, whereas the natural filling factor poip is from 
now on denoted by vg- Explicitly, the relation between 
the two filling factors reads 



(246) 



Early transport measurements in exfoliated graphene 
on a Si02 have revealed broken s pin- valley-symmetry 
states at at yg = 0,±1 and ±4 ( Jiang etd] . l2007bl : 
IZhang et all |2006| ) , where the latter corresponds to the 
LLs ±1. More recent experiments on exfoliated graphene 
on a h-BN substrate have furthermor e revealed quantum 
Hall states at i/g = ±3 (jPean et at . 2011), thus com- 
pleting the full resolution of the spin- valley quartet, not 
only in n = 0, but also in ±1. 

The observed states may generally be understood in 
the framework of the quantum Hall ferromagnetism, but 
the understanding of the situation at vg = requires 



an additional consideration of the subleading external 
symmetry-breaking terms discussed in Sec. IV. A. 41 The 
two-stage picture, which we adopt here based on the 
above discussions, may be summarized as follows, (a) 
The quantum Hall ferromagnetic states are formed to 
minimize the leading energy given by the Coulomb in- 
teraction. However, because of the (approximate) SU(4) 
symmetry of the interaction, the orientation of the quan- 
tum Hall ferromagnets is not fixed - a polarization in 
the spin channel is as probable as one in the valley chan- 
nel, and this yields the high degeneracy of the Goldstone 
modes described in Sec. IV.B.ll (b) Therefore, in spite 
of the small energy scale of the external fields, the latter 
are relevant for the orientation of the ferromagnets and 
for the degeneracy lifting of the Goldstone modes. 



a. The quantum Hall effect at vq = ±1. For vq = — 1, 
only one spin-valley branch is completely filled by elec- 
trons.'^^ The Zeeman effect would give a small ener- 
getic advantage to spin-^ electrons, such that the two 
spin Goldstone modes associated with collective excita- 
tions to the spin-'f acquire a g = gap, given by A^. 
In contrast to the spin excitations, the Goldstone mode, 
which couples the two valleys in the spin- J, branch of 
n = 0, remains gapless, and the ground state may thus 
be viewed as a valley-pseudospin ferromagnet in the spin- 
\. branch. The activation gap would be given by Eq. 
PSSP for pseudospin skyrmion-anti-skyrmion pairs, and 
its associated scaling Jeln oc ^fB has ind eed been ob- 
served experimentally ( Jiang et alV l2007b() . The resid- 
ual valley SU(2) symmetry may be broken by the lat- 
tice distortions, which we have discussed in Sec. IV. A. 41 
Whereas an out-of-plane lattice distortion would yield a 
gapped valley-pseudospin-wave mode, a Kekule-type in- 
plane distortion orients the pseudospin ferromagnet in 
the X F-plane, associated w ith a gapless U(l) superfluid 
mode ( Nomura et al\ . [2009() . Notice that the lattice dis- 
tortion, characterized by the energy scale A^efc is not 
in competition, at vq = ±1, with the Zeeman effect, 
such that the resulting ferromagnetic state is the same 
for > Afcefe as for A^ < ^kek- In the remainder 
of this section we restrict the discussion of the valley- 
pseudospin degeneracy lifting to in-plane distortions that 
seem to be energetically more relevant than out-of-plane 
distortions, but the overall picture remains unchanged if 
the latter are more relevant. 



b. The quantum Hall effect at vg = 0. The situation is 
more subtle at = 0, where it is not possible to fully po- 
larize both the spin and the valley pseudospin and where 
the Zeeman effect is in competition with a lattice distor- 



For VQ = "1-1 J the same arguments apply in terms of holes due 
to particle-hole symmetry. 



47 




(b) 





Figure 27 Possible scenarios for the lifted spin-valley degeneracy 
at VQ = 0. (a) Az > A^efe in the bulk. When approaching 
the edge, the energy difference between the two valleys increases 
drastically, and two levels {K' and {K,^) cross the Fermi energy 
at the edge depicted by the dashed line (Quantum Hall state), (b) 
Afeefc > A^ in the bulk. The K subbranches are already located 
above the Fermi energy, and those of K' below, such that the energy 
difference is simply increased when approaching the edge with no 
states crossing the Fermi energy (Insulator). 



tion that orients the valley pseudospin. For > A^efc, 
it is favorable to fill both valley sublevels of the spin- 
^ branch and the resulting state is a spin ferromagnet 
with gapped spin- wave excitations. For < A^efe, a 
pseudospin-ferromagnetic ground state is favored with 
both spin sublevels completely filled. The two different 
situations are depicted in Fig. [^Tj Most saliently, the two 
phases reveal drastically different transport properties as 
one may see from their behavior at the sample edges. 

The electronic behavior at the edges may be described 
within a model of electron confinement, in which the 
sample edge is described via a mass confinement term 
M{y)T^g in the Hamiltonian, which has the symmetry 
of the term ()218p or else, in n = 0, that of a valley- 
Zeeman term (|217p . as argued in Sec. IV.A.4I The pa- 
rameter M{y) is zero in the bulk and increases drasti- 
cally at the edge at a certain value of the coordinate 
y.^^ Although the model is a simplification to treat the 
graphcne edges in the continuum description of the Dirac 
equation, a more sophisticated treatment that takes into 
account the geometry of the edges yields, apart from a 
fine structt ire of the levels a t the e dge, qualitatively simi- 
lar results ( Brev and FertisJ . |2006() . The mass term M{y) 
modifies the valley coupling due to the lattice distortion 



and yields a y-dependent term -\/A|^/j + M{y)'^, which 
therefore equally diverges at the sample edge.'^^ 

These preliminary considerations on the gap behav- 
ior at the edges allow us to appreciate the difference in 
the expected electronic transport between a spin ferro- 
magnet and a valley-pseudospin ferromagnet at i>g ~ 0. 
Indeed, for Az > Afcefc, one obtains a quantum Hall 



^® For the present argument, we consider translation invariance in 
the x-direction. 

^"^ In the case of an out-of-plane distortion, the term M(y) simply 
adds up to the energy scale A™""^^ [see Eq. l|217| l] . but the 
physical picture remains unaltered. 



state at = that is characterized by a bulk gap as- 
sociated with two counter-propagating edge states [Fig. 
W\si)]. In the bulk, where M[y) = 0, both valley sub- 
levels of the spin- J, branch are filled (spin ferromagnet). 
When approaching the edge, however, the energy term 
\J A^,g^, + AI{y)'^ is enhanced by the rapidly increasing 
contibution from M{y), and the (i^, t) level eventually 
crosses the (K ' j) one at the edge at the Fermi energy 
( Abanin et all l2007bl ). This situation corresponds to 



a quantum Hall state with a bulk insulator and (two 
counter-propagating) conducting channels. In contrast 
to usual quantum Hall states, the edge states are not chi- 
ral, but the chiralitiy, i.e. the transport direction, of each 
channel is linked to its spin orientation.'^^ The quantum 
Hall state therefore remains stab le unless magneti c impu- 
rities couple the two chiralities ( Shimshoni et al\ . [2009i) . 
One notices furthermore a change in the spin polariza- 
tion at the edge; whereas the spin polarization in the 
bulk is complete, the system becomes spin unpolarized 
at the edge. If one takes into account the exchange inter- 
action, the change in the polarization takes place over a 
certain distance, and the conducting properties may be 
described in terms of spin-carrying one-dimensional edge 
excitations ( Shimshoni et al ]. l2009l) . 

In the opposite limit with Az < Akek in the bulk [Fig. 
[27l b)]. the system at i^g = is already valley-polarized, 
and an increase of y^l^^^T+iW^^/p when approaching 
the edge does not induce a level crossing at the Fermi 
energy. Thus, there are no zero-energy states at the edge, 
and the system would be insulating both in the bulk and 
at the edge. 

From an experimental point of view, it is not fully set- 
tled which of the two phases describes the state at i^g = 
0. Whereas early experiments in exfoliated graphene on 
Si02 samples were di scussed in the f ramew o rk of a dom- 
inant Zeeman effect (I Abanin et all l2007bl : iJiang et all 
l2007bl : IZhang et al. more recent experiments at 

very large magnetic fields ( Checkelskv et am2008l . [20091 ) 
and on suspended grap hene samples with increasing mo- 
bility (|Du et all l2009( ) favor the insulating scenario of 
Fig- I27f b) with a dominant valley degeneracy lifting. 
Especially the high-field measurements hint at an easy- 
plane or XY (valley-pseudospin) ferromagnetic ground 
state because the transition between the metallic and the 
insulating state i s reminiscent of a Koste r litz-T houless 
phase transition ( Kosterlitz and Thoulessl . Il973l ) if one 
replaces the temperature by the magne tic field as the 
parameter driving the transition ([C hecke lskv et all l2008l 
l2009l : lHou et adboiOtlNomura etallhOO^ . However, it 
has also been argued that this effect may be understood 
within the above scenario of a Zeeman-dominated quan- 
tum Hall ferromagnet in the bulk, in the framework of a 
Luttinger-liquid description of the domain wall seperat- 



These helical edges are the signature of a quantu m spin Hall 
effect llHasan and Kand [2OI0I: iQi and Zhanel . l20lih . 
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i ng the polarized f rom t he unpolarized region at the edge 
(|Shimshoni et aZj . |2009| ) . 

One notices that both the Zecman effect and the 
Kekule-type lattice distortion are very close in energy 
(see Tab. II) such that one may speculate that other 
effects, as e.g. impurities, strongly influence the forma- 
tion and the orientation of the quantum Hall ferromag- 
net. Further experimental and theoretical studies there- 
fore seem to be necessary to clearly identify the leading 
symmetry-breaking mechanisms, which need not be uni- 
versal, in the zero-energy LL at = and ±1. 

We finally mention scanning-tunneling spectroscopic 
results for the level splitting at z^g = that were per- 
formed on graphene on a graphite substrate (|Li et 
l2009a[ ). Although a gap has been observed as one may 
expect in the framework of the above scenario, it satu- 
rates as a function of the magnetic field. This is in dis- 
agreement with both the -s/S-behavior of an interaction- 
dominated gap as well as with the linear dependence of 
the Zeeman or lattice-distortion effects. A probable ori- 
gin of this gap is the commensurability of the graphene 
lattice with the graphite substrate that may break the in- 
version symmetry between the two sublattices by a term 
of the type (|218p . The coupling to the substrate being es- 
sentially electrostatic, one would then expect no or only 
a weak magnetic-field dep endence o f the s plitting, as ob- 
served in the experiment ( Li et all . l2009al) . 



C. Fractional Quantum Hall Effect in Graphene 

The most salient aspect of strongly-correlated 2D elec- 
trons in partially filled LLs is certainly the FQHE, which 
is due to the formation of incompressible liquid phases at 
certain magical values of the filling factor. As we have 
already argued in Sec. IV. A. 31 on the basis of the pseu- 
dopotentials, the FQHE is expected to be present in the 
graphene LLs n ~ and n = 1, and the main differ- 
ence with respect to non-relativistic 2D electron systems 
should arise from the internal S U(4) symmetry [for a re- 
cent theoretical review see Ref. (|Papic et a/.l . l2009f )] . 

On the experimental level, recent progress in the fabri- 
cati on of high-mobility samples, e . g. via current anneal- 
ing (|Bolotin aLl . 120081: IDu et adl2008( ). has allowed for 
the observation of several FQHE states in graphene. The 
first observations of a state at z^g = ±1/3 were reported 
in 2009 on current-annealed suspen ded graphene s a mples 
in the two-term inal configuration (iBolotin et al I 120091: 
iDu et aU [20091 ).^^ More recently, in 2010, the FQHE has 
also been observed in the four-terminal geometry, which 
allows for the simultaneous measurement of the longi- 
tudinal and the Hall resistance, in suspended graphene 
(jGhahari et adl201l[) and on graphene on a h-BN sub- 



strate (iDean et a/.l . l201ll) . 

Before commenting in more detail on these first exper- 
imental results (indeed, this part of graphene research 
has just started), we introduce the theoretical four- 
component or SU(4) picture of the FQHE in graphene, 
in terms of generalized Halperin wave functions. These 
wave functions, which may be viewed as multi-component 
generalizations of Laughlin's wave function, provide the 
natural framework for the description of the phenomenon 
in view of the model of electrons restricted to a single rel- 
ativistic LL (Sec. IV.Ap 



] 1. Generalized Halperin wave functions 



The theoretical study of the FQHE is intimitcly related 
to trial TV-particle wave functions. In 19 83, Laughlin pro - 
posed a one-component wave function ( Laughlim Il983() . 



N 



({^4)-n( 



Zk- zi) 



(247) 



k<l 



which allows for an understanding of incompressible 
FQHE states at the filling factors ly = 1/m that are de- 
termined by the exponent m for the particle pairs fc, I in 
Eq. (|247p . The variable Zk = {xk~iyk)/lB is the complex 
position of the fc-th particle, and the form of the Laughlin 
wave function (I247P is dictated by the analyticity condi- 
tion for wave functions in the lowest LL.'**' Furthermore, 
the exponent m must be an odd integer as a consequence 
of the fcrmionic statistics imposed on the electronic wave 
function. Even if Eq. ([247]) describes only a trial wave 
function, one can show that it is the exact ground state 
for a class of model interactions tha t yield, with th e help 
of Eq. ((2T5)) . the pseudopotentials (jHaldanel . [l983[) 



Vi>0 for e <m and Ve = for £ > m. (248) 

Although the Coulomb interaction does not fulfill such 
strong conditions, the pseudopotentials decrease as 
l/\/rn for large values of m. Because the incompressible 
ground state is protected by a gap that is on the order 
of Vi, one may view the pseudopotentials Ve>m as an 
irrelevant perturbation that does not change the nature 
of the ground state. Indeed, exact-diagonalization calcu- 
lations have shown that, for the most prominent FQHE 
sX V = 1/3, the overlap between the true ground state 
and the Laughlin state (l247l) is extremely larg e (> 99%) 
(|Fano et a?.l . ll986l : iHaldane and RezaviL Il985l ). 



Soon after Laughlin's original proposal, Halperin gen- 
eralized the wave function (|247p to the SU(2) case of 



^® There are also some weak indications for FQHE states at other 
fiUing factors than vq = ±1/3 in these samples. 



The lowest-LL condition of analytic wave functions may seem a 
very strong restriction when discussing FQHE states in higher 
LLs. However, the model II207II indicates that all LLs can be 
treated as the lowest one, n = 0, if the interaction potential is 
accordingly modified. We adopt this point of view here. 
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electrons with s pin, in the absence of a Zeeman effect 
( Halperiiil . Il983( ) - one has then two classes of particles, 
spin-f and spin- J, particles, which are described 

ft") (1) 

by the complex positions zjj^ and zjj^ , respectively. In 
the (theoretical) absence of interactions between elec- 
trons with different spin orientation, the most natural 
ground-state candidate would then be a simple product 
of two Laughlin wave functions (|247p . 

(249) 

one for each spin component with the exponents m-f and 
m^, respectively, that need not necessarily be identical. 
Inter-component correlations may be taken into account 
by an additional factor 



(T) 



(250) 



where the exponent n can now also be an even integer be- 
cause the fcrmionic anti-symmetry condition is concerned 
only with electrons in the same spin state. 

Halperin's idea is easily generalized to the case of more 
than two components, and the corresponding trial wave 
function for an SU(A^) quantum Hall system w ith Af com- 
ponents reads (|Goerbig and Regnaultl . I200l1) 

(251) 

in terms of the product 



''m-i,...,mu 



(252) 



i=l kj<lj 

of M Laughlin wave functions and the term 



nnn 

i<j ki kj 



(253) 



which describes inter-component correlations. As in 
the case of Ha lperin's two-component wave function 
(iHalperinl . I1983D . ;he exponents mj must be odd inte- 



gers for fermionic particles whereas the exponents Uij 
may also be even integers. These exponents define a 
symmetric JV x J\f matrix M = riij, where the diago- 
nal elements are nu = nii. This exponent matrix encodes 
the statistical properties of the quasi-particle excitations, 
such as their (fractional) charge and their statistical an- 
gle (IWen and Zeel . ll992allbh . 

Moreover, the exponent matrix M determines the com- 
ponent densities pj - or, equivalently, the component fill- 
ing factors Vj = pj/ns, 



(254) 
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where u = vi + . . . + is the total filling factor mea- 
sured from the bottom of the lowest LL. Notice that Eq. 
(|254p is only well-defined if the exponent matrix M is in- 
vertible. In this case, all component filling factors Vj are 
completely determined, whereas otherwise some of the 
component fillings remain unfixed, e.g. vi and V2 for the 
sake of illustration, although the sum of them [vi + V2) 
is fixed. This is nothing other than a consequence of the 
underlying ferromagnetic properties of the FQHE state 
that, similarly to the states ai v = k discussed in Sec. 
IV. Bl are described by subgroups of SU(A/'). 

Finally, we notice that not all SU(A/') wave functions 
(|25ip describe incompressible quantum liquids with a ho- 
mogeneous charge density for all components. A general- 
ization of Laughlin's plasma picture, according to which 
the modulus square of the trial wave function corre- 
sponds to t he Bo ltzmann weight of a classical 2D plasma 
(iLaughlinl . 1X9831) . shows that all eigenvalues of the ex- 
ponent matrix M must be positive (or zero for states 
with ferromagnetic order) . Otherwise, some of the differ- 
ent components phase-separate in the 2D plane because 
the inter-component repulsi on between them ex ceeds the 
intra-component repulsion ( de Gail et ail . [20081 ) . 



2. The use of generalized Halperin wave functions in graphene 

These general considerations allow us to define the 
framework for a basic description of the FQHE in 
graphene where the SU(4) spin- valley symmetry imposes 
M ~ A. Four-component Halperin wave functions are 
therefore expected to play an equally central role in the 
description of the graphene FQHE as Laughlin's in a one- 
component or Halperin's in two-component systems. In 
the remainder of this section, we attribute the four spin- 
valley components as 1 = (f, i^), 2 = (f, K'), 3 = (4-, A'), 
and 4 = {i,K'). 



a. Fractional SU(4) quantum Hall ferromagnet. In a first 
step, we consider a four-component Halperin wave func- 
tion in which all components are equal (odd) integers, 
nij = Hij = m, regardless of whether they describe 
intra- or inter-component correlations. One obtains then 
a completely anti-symmetric orbital wave function that 
is accompanied by a fully symmetric SU(4) spin-valley 
wave function. 

As we have argued in Sec. IV.B.li this situation repre- 
sents precisely a perfect SU(4) quantum Hall ferromagnet 
- indeed, for m = 1, the generalized Halperin wave func- 
tion (|25ip is nothing other than the orbital wave function 
of the state at 1/ = 1, i.e. when one of the subbranches is 
completely filled. The SU(4) symmetry is then sponta- 
neously broken, and the group-theoretical analysis pre- 
sented in Sec. IV.B.ll vields 3 degenerate Goldstone modes 
that are generalized spin waves. 

The situation is exactly the same for any other odd 
exponent m, but the orbital wave function (|25ip is then 
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a Laughlin wave function (|247p in terms of the particle 
positions Zk regardless of their internal index j = 1, 4. 
The ferromagnetic properties of these wave functions may 
be described by the same equations as the spin- wave and 
skyrmion modes derived in Sec. IV.Bl if one takes into ac- 
count a renormalization of the spin stiffness, as it has 
been discussed extensively i n the l iteratur e for SU(2) 
quan t um Hall ferromagne ts ( Ezawai. 2000; Moon et all . 
1 19951 : ISondhi et all . Il993l ). States described by such a 
wave function are ground-state candidates for the filling 
factors ly = 1/m, which correspond to the graphene fill- 
ing factors [see Eq. (|246p ] ug = —2 + 1/m or hole states 
at i^G = 2 — 1/m. 

There are now two different manners to break the in- 
ternal SU(4) symmetry explicitly. The simplest one is 
the same as for the quantum Hall ferromagnetism at 
I'G = or ±1. in terms of external symmetry- breaking 
fields such as those discussed in Sec. IV. A. 41 However, 
one may also change some of the exponents in the gen- 
eralized Halperin wave function (|25ip . in which case one 
also changes the filling factor. One may for instance con- 
sider the [m; m — 1, m] wave function with mj = m for all 
j, ni3 = 7124 = m - 1, and ni2 = riu = ^23 = ^34 = m, 
which correspond to a filling factor^^ 

2 

or i^G = - 



2m - 1 



2m- 1 



(255) 



Indeed, the difference in the inter-component exponents 
explicitly breaks the spin-valley symmetry - electrons in 
different valleys are more weakly correlated (with an ex- 
ponent m — 1) than electrons in the same valley (expo- 
nent m), regardless of their spin orientation. As a con- 
sequence, the filling factors in each of the two valleys, 
i^K ^ vi + and vk' ^1^2 + ^4, respectively, are fixed, 
i^K = i^K' ~ l/(2m — 1), and one may view the wave 
function as a state with ferromagnetic spin ordering, but 
that is valley-pseudospin unpolarized. Alternatively, the 
[m; m— 1, m] wave function may be interpreted as a tensor 
product of an SU(2) Halperin ( r n.m. rn — 1) pseudospin- 
singlet wavefunction (iHalperinl . [l98l and a completely 
symmetric (ferromagnetic) two-spinor that describes the 
physical spin. The relevance of the [m;m — l,m] wave 
function, with m = S {v = 2/5) has been corroborated in 
recent exact-diagonali zation studes, both in the graphene 
LLs n ~ and 71 = 1 ( Papic et all . l2009l : iToke and Jainl . 
I2OO7D . 

The SU(4) spin- valley symmetry is fully broken, e.g., 
in the case of the [m; m — 1, m — 1] wave function with 



all m. 



m and off-diagonal 



function, which describes a state at 



or VG 



4m- 3 



1. This wave 



(256) 



4m — 3 



'^^ We only discuss electronic states here, but the arguments are 
equally valid for the particle-hole symmetric states at = 2 — 
2/(2m - 1). 



may be viewed as an SU(4) singlet where the filling fac- 
tors of all spin-valley components are l/(4m — 3). Exact- 
diagonalization calculations for = 4 and 8 particles 
have shown that the [m; m — 1, m — 1] wave function with 
m = 3 {'At V = 4/9) describes to great accuracy the 
ground state for a Coulomb interaction (|205p . with over- 
laps On=8 = 0.992 in n = and On=r = 0.944 and in 
the n = 1 graphene LL ( Papic et a/.l . l2009l) . 



b. A route to understanding the graphene FQHE at vq = 
±1/3. The discussion of the above-mentioned states was 
based on the understanding acquired from quantum Hall 
systems in semiconductor heterostructures, where the fill- 
ing factor is defined with respect to the bottom of the 
n = LL. First experimental observations, however, in- 
dicated a prominent FQHE at vg — ±1/3, which cor- 
responds to two completely filled spin-valley sublevels of 
the graphene 71 LL, and a third one that is 1/3 filled, 
v = 2 + 1/3. Such a state would naturally arise in a 
system where the SU(4) symmetry is strongly broken, 
e.g. by a strong Zeeman effect. However, as argued in 
Sec. IV.A.4[ these external fields are weak as compared 
to the leading interaction energy scale, and it is therefore 
natural to ask how such a state may arise from the inter- 
action point of view in the framework of four-component 
Halperin wave functions. 

A Halperin wave fu nction that describ es the above- 
mentioned situation is ( Papic et aLl . l2010l ) 



^2+1/3 = n n 

i=K,K' i<] 

X n n 

^=K,K' i<j 



(257) 



or any permutation of the spin-valley components. One 
notices that this wave function implicitly breaks the 
SU(4) spin-vallcy symmetry and, moreover, is not an 
eigenstate of the full SU(4) pseudospin, such that it can- 
not describe the ground state in the total absence of 
an external symmetry-breaking field. However, exact- 
diagonalization calculations have shown that even a 
tiny external Zeeman field is capable to stabelize the 
state (|257p , which becomes th e ground state for — 
O.OleVds (jPapic et all . l2010[ ). Furthermore, the state 
(|257p possesses, in addition to the valley-pseudospin- 
wave Goldstone mode in the spin-f branch, low-lying 
spin-flip excitations for moderately small Zeeman fields, 
even if the charge (activation) gap is the same as for the 
usual 1/3 Laughlin state. This particular interplay be- 
tween the leading Coulomb energy and subordinate ex- 
ternal spin- valley symmetry breaking terms, illustrated 
at the VG = 1/3 example, shows the complexity of the 
graphene FQHE, and further surprises may be expected 
in future experiments. 
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3. Experiments on the graphene FQHE 

We terminate this section on the graphene FQHE with 
a short discussion of experimental observations in the 
Ught of the above-mentioned theoretical four-component 
picture. 



a. Two-terminal measurements. In the first observations 
of the FQHE, the two-terminal configuration was used, 
where the voltage (and thus the resistance) is measured 
between the same two contacts us ed to drive th e elec- 
tric current th rough the sample (iBolotin et al. I, 120091: 
In this two-terminal configuration, it 
is not possible to measure simultaneously the Hall and 
the longitudinal resistance. It is nevertheless possible to 
extract the Hall and the longitudinal conductivities from 
the two-terminal resistance with the help of a conformal 
mapping, as a consequence of the 2 T) nature of the quan- 
tum transport in these s ystems ( Abanin and Levitovl . 
l2008l : I Williams et a/.l . [20091 ) . This technique has been ap- 
plied to obtain insight into the longitudinal conductivity 
the expected activated behavior of which yields a rough 
estimate of the activation gap at vg = 1/3 (A1/3 ^ 4.4 K 
at B = 12 T) (j Abanin et all\201(J ). which is an order of 
magnitude smaller than the theoretically expected v alue 
(|Apalkov and Chakrabortvll2006l: iToke et al. I, [200^.42 



b. Four-terminal measurements. The activation gap of 
the 1/3 FQHE state has also been measured in suspended 
graphene in the four-terminal configuration, in which 
the longitudinal resistance can be measur ed directly and 
indep endently from the Hall resistance ( Ghahari et al\ . 
I2OIII ). In this case, the activation gap has been estimated 
to be Ai/3 ~ 26. ..50 K at B = 14 T, a value that agrees 
reasonably well with the t heoret i cally expec t ed va lue 
(jApalkov and Chakrabortvl l2006t iToke et all |2006[ ) if 
one considers the energy scale e'^/eSoolB, which takes into 
account the RPA dielectric constant Eoo for graphene in 
vacuum (see Sec. IIII.B.4[) . 

Finally, we would mention very recent high-field trans- 
port measurements in the four -terminal configur ation on 
graphene on a h-BN substrate (jPean et al ]. l201lh . These 
experiments allowed for a clear identification of several 
states of the 1/3 family, at lyg = ±1/3, ±2/3, and ±4/3 
corresponding to the zero-energy LL n = 0, as well as 
at = ±7/3, ±8/3, ±10/3, and ±11/3 which reside in 



Notice that the theoretical estimates have been obtained within 
a simplified two-component model, with a completely frozen 
spin degree of freedom. In spite of this simplification, the 
above-mentioned exact-diagonalization calculations with an im- 
plemented SU{4) symmetry have shown that the charge gap, 
which is responsible for the activated behavior, coi ncides indeed 
with that obtained in the two-component model l|PaDic et al\ . 
l201Cf) . 



the n = 1 LL. The estimation of the activation gap at 
i/G — 4/3 agrees again reasonably well with the theo- 
retical expectation for the 1/3 state. The perhaps most 
salient (and unexpected) feature of the transport mea- 
surement is the absence (or extreme weakness) of the 
vg = ±5/3 representative of the 1/3 family, which would 
correspond to the Laughlin state (i^=l/3-r->z^G = — 5/3 
and the corresponding hole state) with a full SU(4) spin- 
valley ferromagnetic order, as argued in Sec. IV.C.2I 

Whereas the absence of this state remains to be un- 
derstood, these findings corroborate the theoretical four- 
component picture of the graphene FQHE. Indeed, it 
clearly shows that the SU(4) symmetry of the n = 
LL is essential because the only correspondence between 
the FQHE states is particle-hole symmetry that maps 
i>G "i^G- If the SU(4) symmetry were broken, e.g. 
by a sufficiently strong Zeeman effect, the only symme- 
try would be the valley-SU(2) symmetry in each spin 
branch of the n = LL, in which case there exist the 
further mappings — 2 + -s-> — in the spin-^ branch and 
1/ O 2 — 1/ in the spin-f branch. However, the (observed) 
±1/3 state would than be mapped on the (unobserved 
or extremely weak) ±5/3 state, and the strong difference 
in the visibility between these two states is therefore dif- 
ficult to understand. This is also the case if the SU(4) 
symmetry is fully broken by strong external spin and 
valley Zeeman fields, such that all spin- valley sublevels 
are completely resolved, and ±5/3 would be mapped on 
±4/3, in the same manner as ±1/3 on ±2/3. 



VI. CONCLUSIONS AND OUTLOOK 

We have reviewed the quantum-mechanical properties 
of relativistic 2D electrons in monolayer graphene ex- 
posed to a strong magnetic field. The main parts of this 
review are concerned with the role of electronic interac- 
tions in graphene LLs. Whereas we have argued that 
these interactions may be treated perturbatively in the 
regime of the relativistic (integer) quantum Hall effect 
(RQHE) , they constitute the relevant energy scale in par- 
tially filled graphene LLs due to the quenching of the ki- 
netic energy. This is reminiscent of partially filled LLs in 
non-relativistic 2D electron systems, and the most promi- 
nent consequence of this quenched kinetic energy and the 
macroscopic LL degeneracy is certainly the FQHE. The 
graphene FQHE is expected to be reminiscent of that of 
non-relativistic 2D electrons but it is governed by a larger 
internal degeneracy described to great accuracy by the 
SU(4) group. The experimental study of the FQHE in 
graphene is still in its infancy, and novel surprises may 
be expected. Only recently have been reported mea- 
surements in the four-terminal geometry which allow for 
an analysis of prominent characteristics of FQHE states, 
such as the activation gaps. In view of the generally 
accepted universality of the quantum Hall effect, it will 
certainly be interesting to make a systematic compari- 
son with the activation gaps of related FQHE states in 
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conventional 2D electron gases with a parabolic band. 

In the perturbative regime of the RQHE, the theoret- 
ical study of electron-electron interactions indicates the 
presence of fascinating novel collective modes, such as 
linear magneto-plasmons, that are particular to graphene 
and do not have a counterpart in non-relativistic 2D elec- 
tron systems in a perpendicular magnetic field. Also the 
upper-hybrid mode, which is the magnetic-field counter- 
part of the usual 2D plasmon, is expected to behave in 
a particular manner in graphene as a consequence of the 
linear disperison relation and the vanishing band mass. 
Whereas these studies are at present only theoretical, 
these collective modes may find an experimental verifica- 
tion in inelastic light-scattering measurements. 

Similarly to the role of electron-electron interactions 
in the RQHE regime, the electron-phonon coupling 
yields exciting resonance phenomena in graphene in a 
strong magnetic field. The electron-phonon interaction 
in graphene LLs has been discussed in the framework of 
a perturbative approach. Indications for the magneto- 
phonon resonance, e.g., have recently been found in Ra- 
man spectroscopy of epitaxial graphene. 

The present review has been limited to monolayer 
graphene, and it is definitely a reasonable research pro- 
gram to ask how the effects described here manifest them- 
selves in bilayer graphene. For example, the particular 
collective excitations described in Sec. IIIII have been at- 
tributed to the lack of equidistant LL spacing and the 
presence of two bands. Whereas bilayer graphene also 
consists of two (particle-hole-symmetric) bands in the 
low-energy regime, the approximate parabolicity there 
yields almost equidistant LLs. The presence of additional 
high-energy bands (in the 300 meV range) certainly also 
affects the plasmonic modes. 
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Appendix A: Matrix Elements of the Density Operators 

The matrix elements that intervene in the expres- 
sion for the density operators (jl39p are of the form 
(n, m|exp(— iq • r)|n',m') and may be calculated with 
the help of the decomposition of the cyclotron variable r] 
and the guiding centre R into the ladder operators a and 
b, respectively [see Eqs. dTS]), (HH) and We further- 

more define the complex wave vectors q = {qx + iqy)lB 
and q~{qx~ iqy)lB,'^^ One finds 

(n,TO|e-'i""|n',m') ^ (m|e-'i-^|m') ® (nle^'i'^ln') 



(m|e 



TO') (Al) 

n'). 



The two matrix elements may be simplified with the 
help of the Baker-Hausdorff formula exp(^) exp(i?) = 
exp(A -I- B) cx p{\A,B]/2), for the case \A [A, B]] = 
[B, [A, B]] = (|Cohen-Tannoudii et a/.l . Il973[) . The sec- 
ond matrix element thus becomes, for n > n' 



(gat+ga)|^^, 



(A2) 



e-l«|-/4(n|e-^*"" e~^*"|n') 
e-l«l^/^5](n|e-7f'^'^*|j)(j|e-7f''>') 




^ (n- j)!(n' - j)!j! 



2 



2 \ " -J 




■ -\ n—7i /I |9 

nW^J ^2 



where we have used 



(n|e-75«-<^'|j-) 



%\ 1 



jl. {n~jy. \ 72 



for j > n 

n-j 

=q for j < n 



in the third line an d the definition of the associate d La- 
guerre polynomials ( Gradshtevn and Rvzhikl . l2000i ). 



i-xY 



{n' — my.{n — n' + to)! to! 



We use this notation solely in the present appendix. Throughout 
the main text, q denotes the modulus of the wave vector q, g = 

|q|- 



53 



In the same manner, one obtains for m> m! 
(m|e-*i-^|m') = (m|e-7!(?^'+«^) |,7i') 

j 77 / . \ m—m' 

= e-l«l'/s/— ( ^) 
V m! Vy2/ 

xi"r"' (^M!^ . (A3) 
With the help of the definition 

one may rewrite the expressions without the conditions 
n> n' and m > m' , 

{n\e-'^-'l\n') ^ [e(n - n')a„,„' (g) (A4) 
+e(n'-n-l)e„^„(-g)]e'l«l'/4 

and 

(m|e-*i-^|m') = [e(m - m')^,„,™' (g) (A5) 
+e(m' - m - e-l-'l'/^. 
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